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FOREWORD 


The  US  Army  Electronic  Proving  Ground  was  responsible  for  the  execution 
of  this  project.  Mr.  C.  M.  Giorgi  was  the  USAEPG  Project  Officer  with  as¬ 
sistance  from  Mr.  L.  Butler  and  1LT  G.  Christman.  Technical  guidance,  as  well 
as  the  idea  for  the  project,  was  provided  by  the  Test  and  Evaluation  Command 
Project  Officer,  Dr.  N.  E.  Pentz. 

Dr.  E.  A.  Milne  and  Dr.  E.  C.  Crittenden  of  the  Naval  Postgraduate 
School,  Monterey,  California,  were  responsible  for  the  overall  technical 
direction  of  the  project,  to  Include  test  planning,  data  acquisition,  and  data 
analysis.  The  material  presented  in  the  body  of  this  report  primarily  re¬ 
flects  the  results  of  their  work  as  documented  by  Dr.  Milne. 

Dr.  J.  Engel  and  Mr.  R.  Ransier,  Bell  Technical  Operations,  Textron, 
Tucson,  Arizona,  provided  support  for  the  operation  of  the  video  data 
acquisition  systems  and  analysis  of  the  video  data.  The  documentation 
covering  the  video  data  acquisition,  processing,  and  the  associated  computer 
programs  was  prepared  by  Dr.  Engel. 

Mr.  E.  Burgess  of  the  Dugway  Proving  Ground,  Utah,  provided  the  smoke 
expertise  for  the  project  and  also  directed  its  deployment  during  the  test. 
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SECTION  1 .  SUMMARY 


1.1  BACKGROUND 

a.  The  spot  patterns  produced  by  laser  designators  suffer  from  a  number 
of  defects.  The  patterns  are  broadened  and  exhibit  wander  and  intensity  fluc¬ 
tuation  because  of  turbulence  in  the  atmosphere.  Designator  spot  patterns 
also  exhibit  similar  defects  because  of  internal  laser  instability  and  motion 
of  the  transmitter  platform.  In  order  to  correct  the  latter  defects,  a  system 
is  needed  that  is  capable  of  separating  the  contributions  of  the  several  ef¬ 
fects. 


b.  Investigative  effort  at  the  Naval  Postgraduate  School  (refs  1,  2,  and 
3,  app  G)  and  elsewhere  (ref  4,  app  G)  previously  demonstrated  that  the  ef¬ 
fects  of  atmospheric  turbulence  on  designator  spot  patterns  can  be  expressed 
in  terms  of  the  Optical  Transfer  Function  (OTF)  of  the  atmosphere.  This 
quantity  is  the  Fourier  transform  of  the  spot  profile  due  to  the  atmosphere 
alone.  If  other  causes  of  spot  broadening  such  as  laser  instability  and  plat¬ 
form  motion  are  present,  the  spot  profile  will  be  the  convolution  of  the_spot 
profiles  due  to  each  of  the  effects.  The  spot  profile  broadening,  due  to 
laser  and  platform  Instability,  can  then  be  separated  from  the  atmospheric  ef¬ 
fects  by  dividing  the  Fourier  transform  of  the  spot  profile  by  the  OTF  of  the 
atmosphere  ( polnt-by- point  as  a  function  of  spatial  frequency).  The  necessary 
OTF  of  the  atmosphere  can  be  obtained  by  means  of  a  slit-scanning  telescope 
system  which  views  a  point  laser  in  the  target  vicinity  from  a  location  near 
the  designator  transmitter.  An  additional  imaging  system,  viewing  the  de¬ 
signator  spot  on  a  target  screen,  from  a  location  near  the  target,  yields  the 
composite  spot  profile.  Data  from  this  Imager,  together  with  that  from  the 
telescope  imager  at  the  transmitter  site,  can  separate  the  various  effects. 
However,  because  designator  lasers  use  short  pulses,  the  previously  developed 
mechanical  slit-scan  techniques  cannot  be  used.  Storage  imaging  systems,  us¬ 
ing  VIDICON,  CCD,  or  CID  techniques,  are  needed  to  sense  the  image  for  the 
very  short  pulse  periods  and  provide  the  equivalent  of  a  slit-scan  during  the 
inter pulse  period. 

c.  A  study  and  experimental  Investigation  to  evaluate  and  demonstrate 
the  feasibility  of  the  video  technique  was  accomplished  as  Phase  I  of  this 
investigation.  The  results  of  the  Phase  I  effort  are  reported  in  references  7 
and  8,  appendix  G.  The  feasibility  of  the  technique  was  demonstrated  and  a 
test  plan  to  evaluate  the  technique  under  field  conditions  was  proposed  as  a 
followup  effort  (See  app  E).  This  report  covers  the  results  of  that  field 
test  or  Phase  II  of  the  investigation. 

1.2  OBJECTIVES 

a.  To  compare  atmospheric  turbulence  measurements  made  with  a  mechanical 
slit  scan  telescope  system  to  those  made  with  a  video  imaging  system. 

b.  To  develop  and  document  video  data  reduction  and  analysis  procedures 
and  software  for  measuring  atmospheric  turbulence. 

c.  Apply  the  measured  field  data  to  account  for  the  effects  of  atmos¬ 
pheric  turbulence  on  the  performance  of  a  ground-based  laser  designator. 


1.3  SUMMARY  OF  PROCEDURES 


a.  Test  Period.  The  field  tests  were  conducted  on  10,  11,  and  12 
September  1984.  With  regard  to  the  test  schedule,  the  conduct  of  the  field 
test  involved  the  application  of  instrumentation  and  personnel  resources  which 
were  also  being  utilized  to  support  the  ongoing  Remotely  Piloted  Vehicle  (RPV) 
test  program  at  USAEPG.  Consequently  scheduling  of  the  lower  priority  method¬ 
ology  project  produced  a  series  of  postponements  during  the  investigative 
period  to  Include  having  to  eventually  schedule  around  the  July- August  monsoon 
season  in  southern  Arizona.  As  it  turned  out,  in  consonance  with  one  of 
Murphy's  Laws,  the  first  two  days  of  the  test  produced  all-time  reco-d  rain¬ 
fall  for  those  days. 

b.  Test  Site.  The  test  site  was  on  the  Vest  Range,  Fort  Huachuca,  at  a 
location  previously  used  for  a  "Smoke  Week"  test  program.  The  site  was  chosen 
for  its  clear  line-of-sight  path  and  the  availability  of  a  large  roofed 
instrumentation  platform. 

c.  Instrumentation  Configuration.  The  general  configuration  is  de¬ 
scribed  in  the  test  plan  (app  E).  The  specific  video  measurement  setup  and 
equipment  used  in  the  field  test  is  shown  in  figure  1.  The  MTF  measuring 
instrumentation  was  stationed  on  a  large  platform.  However,  the  vibration 
sensitivity  of  the  measurement  systems  required  the  use  of  a  separate  platform 
for  the  scanning  telescope  and  video  camera/telescope  units  which  were  mounted 
on  tripods.  The  extension  platform  which  was  weighted  down  with  sandbags 
served  to-  Isolate  the  vibration  caused  by  personnel  movement  on  the  main  plat¬ 
form. 


d.  Test  Procedures:  The  test  procedures  are  described  in  the  test  plan. 
Preliminary  checkout  of  all  test  Instrumentation  was  accompllshd  on  the  morn¬ 
ing  of  10  September.  Comparison  data  was  acquired  on  the  11th  and  smoke  test¬ 
ing  completed  September  12th.  A  malfunction  of  the  laser-designator  simula¬ 
tor,  which  had  just  undergone  a  factory  checkout,  precluded  acquiring  data  for 
the  designator  spot  analysis  application  part  of  the  field  test.  The  video 
data  was  recorded  and  later  processed  using  programs  originally  written  by  Dr. 
E.  Milne  of  NPS  for  an  IBM  system  and  adapted  for  the  HP  1000  by 

Dr.  J.  Engel  of  Bell  Technical  Operations.  The  data  from  the  NPS  system  was 
processed  on-site  in  near  real  time  using  a  dedicated  microprocessor.  The 
test  data  acquisition  and  analysis  processes  are  further  described  in  Section 
2  and  in  appendices  A  through  D.  The  appendices  Include  a  complete  descrip¬ 
tion  of  the  computer  programs  used  for  the  data  reduction  and  analysis  and 
provide  a  user-type  manual  for  applying  the  technique  to  support  electro- 
optical  field  measurements. 

e.  Test  Conditions.  As  noted  above,  the  weather  did  not  cooperate  dur¬ 
ing  the  brief  opening  in  the  availability  of  the  test  support  resources.  How¬ 
ever,  on  the  positive  side,  the  three  days  did  include  a  range  of  clear  to 
very  cloudy  conditions  and  a  larger  than  normal  range  of  ambient  temperatures 
due  to  the  uncharacteristic  rain  squalls. 

1.4  SUMMARY  OF  RESULTS 

a.  although  the  total  amount  of  test  data  was  restricted  by  the  weather 
conditions,  the  limited  results  demonstrated  "hat  a  video  system,  viewing  a 
HeNe  laser,  can  be  used  to  measure  Cn2  if  care  is  taken  to  avoid  satura¬ 
tion  of  the  video  signal. 


Figure  1.  Data  acquisition  setup. 


(2)  A  Cn2  data  analysis  program  for  the  video  system  which  runs 
on  the  HP  1000  has  been  developed  and  documented.  This  capability  will  allow 
transfer  of  the  technology  to  a  number  of  potential  users. 

c.  No  experimental  results  of  applying  the  data  for  laser  designator 
performance  analysis  were  obtained  due  to  the  malfunction  of  the  designator 
simulator. 

1.5  ANALYSIS 

a.  Comparison  of  Results.  Although  there  were  differences  in  the  two 
measurements,  the  Cq.  measured  by  the  NPS  slit-scan  system  and  the  Cn  meas¬ 
ured  by  the  video  system  had  a  correlation  coefficient  of  0.81.  This  is  not 
as  good  as  one  would  like  but  it  indicates  that  the  two  systems  are  measuring 
essentially  the  same  thing.  The  lack  of  better  correlation  can  be  attributed 
to  several  problems;  one,  the  video  signal  is  subject  to  saturation  which  will 
give  unpredictable  results;  two,  the  focus  of  the  two  systems  will  give  larger 
Cn  values  if  out  of  adjustment;  three,  the  vibrations  of  the  platform  will 
give  higher  results  for  the  wander  and  the  long-term  measurements. 

b.  Problems  with  the  TV  System 


(1)  Some  observations  with  the  video  system  showed  problems  that 
should  be  corrected  if  the  system  is  to  be  used  effectively.  First,  the 
mounting  of  the  camera  and  the  telescope  should  allow  easy  access  to  the 
adjustments  of  the  polaroids;  second,  a  measurement  of  the  OTF  of  the  video 
optical  system  should  be  carefully  made  in  the  laboratory  before  any  calcula¬ 
tions  are  made.  The  diffraction  limit  OTF  was  used  in  these  calculations; 
third,  the  vibrations  of  the  platform  due  to  the  wind  were  evident  because  the 
Cq,  as  measured  from  the  wander,  was  consistently  high  in  both  the  slit-scan 
system  and  the  video  system.  The  platform  was  braced  with  diagonal  2  x  4s  but 
this  did  not  completely  stop  the  vibration  by  the  wind. 

(2)  The  video  system  needed  a  spotting  scope  on  the  telescope  for 
alignment  and  a  way  to  adjust  the  polaroids  without  taking  the  system  apart. 
The  ability  to  adjust  the  polaroids  is  very  Important  because  the  dynamic 
range  of  the  video  system  is  rather  small.  This  often  leads  to  saturation  of 
the  signal  from  the  video  system  if  the  intensity  of  the  laser  spot  is  too 
high;  small  changes  in  the  extinction  or  Cn2  will  require  that  the  setting 

be  changed  to  prevent  wash-out  of  the  signal  if  the  intensity  is  too  low  or 
saturation  of  the  signal  if  the  intensity  is  too  high. 

c.  Video  System  Viewing  a  Distant  Lrfser 


(1)  The  tests  show  that  a  video  system  viewing  a  distant  laser  could 
be  used  for  making  measurements  on  the  atmospheric  modulation  transfer  func¬ 
tion  (MTF).  However,  there  are  some  details  that  must  be  resolved  before  such 
a  system  can  work  reliably  in  the  field. 
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(2)  The  first  of  these  is  the  focus  of  the  telescope,  which  must  be 
accurately  adjusted  for  the  best  focus.  If  the  system  is  not  in  focus,  the 
measured  values  of  OTF  will  be  high,  especially  during  periods  of  low  turbul¬ 
ence.  This  was  not  much  of  a  problem  for  this  test  since  the  laser  source  and 
the  telescope  and  video  system  were  at  fixed  positions  and  sufficient  time  to 
adjust  for  the  best  focus.  However,  if  either  the  laser  or  the  video  system 
is  on  a  moving  platform,  the  focus  must  be  continually  adjusted.  Associated 
with  this  is  the  need  to  know  the  range,  which  is  necessary  for  the  calcula¬ 
tion  of  the  OTF. 

(3)  Another  detail  is  the  fact  that  the  video  system  has  a  very  low 
dynamic  range  between  the  noise  level  and  the  maximum  intensity  before  satura¬ 
tion.  A  normal  video  system  uses  different  f-stops  to  adjust  the  light  level, 
but  this  would  require  a  different  system  OTF  for  every  f-stop.  Crossed 
polaroids  were  used  for  this  test  which  seemd  to  work  well  except  as  mentioned 
above;  they  were  not  accessible  for  easy  adjustment.  As  a  result,  during 
periods  of  low  turbulence  the  light  level  of  the  laser  spot  on  the  video 
signal  became  too  high  and  needed  to  be  adjusted.  This  was  overlooked  on  part 
of  the  test  and  gave  anomalous  results  such  as  a  negative  OTF. 

(4)  The  problem  of  a  moving  platform  can  be  overcome  by  using  only 
the  short-term  calculations.  Even  with  the  telescopes  on  a  separate  platform, 
there  was  enough  motion  of  the  platform  to  be  noticeable  in  the  wander  meas¬ 
urements  when  the  wind  was  blowing.  If  a  platform  is  being  used,  such  as  a 
helicopter,  a  tracking  system  may  be  necessary  to  keep  the  telescope  pointed 
at  the  laser. 

(5)  With  regard  to  data  turnaround,  the  test  demonstrated  the  need 
for  near  real  time  processing  of  the  video  data.  With  the  currently  available 
powerful  and  Inexpensive  microcomputers,  the  video  data  could  easily  be  pro¬ 
cessed  on  a  near  real  time  basis  with  a  dedicated  unit. 

d.  Comments  on  a  Video  System  Observing  a  Laser  Designator 

(1)  This  part  of  the  test  was  not  run  because  the  laser  designator 
simulator  did  not  operate  during  the  test.  However,  a  comment  can  be  made 
about  the  probable  results  of  such  a  system.  The  divergence  of  the  laser 
designator  beam  must  not  be  greater  than  twice  the  divergence  due  to  the 
atmosphere  on  a  turbulent  day  which  is  on  the  order  of  50  microradians,  in 
order  to  be  able  to  make  measurements  of  Cn^.  It  may  be  possible  to  infer 
the  C,^  from  the  spatial  distribution  of  the  Intensity  about  the  center  of 
the  designator  spot  on  the  target. 

(2)  Because  of  the  failure  of  the  laser  designator  simulator  to  ope¬ 
rate,  no  experimental  results  were  obtained  with  the  designator.  Theoretical 
calculations  indicate  that  a  video  system  viewing  a  laser  designator  spot 
should  work  to  separate  out  the  effects  of  turbulence,  internal  laser  insta¬ 
bility,  and  platform  motion,  provided  that  the  beam  divergence  of  the  designa¬ 
tor  is  of  the  order  of  100  microradians  or  less.  That  is,  the  spot  on  a 
target  1  kilometer  away  should  be  of  the  order  of  10  centimeters  or  less  in 
diameter.  Also  the  video  system  viewing  the  target  should  be  close  enough  so 
that  Che  Image  of  the  spot  is  at  least  15  to  20  pixels  in  diameter. 
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1 .6  CONCLUSIONS 


a.  The  comparison  tests  of  the  NPS  split-scan  system  with  the  video  sys¬ 
tem  demonstrated  that  a  video  system  viewing  a  HeNe  laser  at  0.5  to  2  kilome¬ 
ters  distance  can  he  used  to  measure  the  OTF  of  the  atmosphere  if  care  is 
taken  to  avoid  saturation  in  the  vidicon  tube  and  that  the  best  focus  is 
maintained. 

b.  A  video  system  for  viewing  laser  designator  spots  can  be  used  to 
separate  the  effects  on  laser  designators  of  atmospheric  turbulence,  Internal 
laser  designator  instability,  and  platform  motion,  provided  the  divergence  of 
the  designator  beam  is  less  than  about  100  microradians  when  there  is  no 
atmospheric  turbulence.  The  effects  of  the  turbulent  atmosphere  will  be  less 
than  about  50  aicroradians  which  amounts  to  about  5  cm  at  the  target. 1  km 
away.  This  is  likely  to  be  much  smaller  than  the  effects  of  motion  of  the 
platform. 

1.7  RECOMMENDATION 

If  a  video  system  is  built  for  measuring  the  OTF  of  the  atmosphere,  it  is 
recommended  that  it  include  a  system  for  obtaining  the  best  focus,  a  feedback 
system  for  preventing  saturation  of  the  vidicon  tube,  and  an  on-line  micro¬ 
computer  system  for  immediate  data  reduction. 


SECTION  2.  DETAILS  OF  INVESTIGATION  (FIELD  TEST) 


2.1  EQUIPMENT  DESCRIPTION 

a.  CW  Lager  Source.  A  1 -milliwatt  helium-neon  laser  source  (  «  0.6328 

micrometers)  was  positioned  1147  meters  downrange  from  the  platform  as  measured 
by  a  laser  ranging  system.  The  laser  was  mounted  on  a  tripod  and  pointed 
directly  at  the  camera  on  the  observation  platform.  In  front  of  the  laser  was 
an  8-foot-square  board  painted  white  with  a  small  hole  in  the  center  as  a 
passthrough  for  the  laser  beam.  A  diverging  lens  was  used  to  reduce  the  inten¬ 
sity  of  the  laser  and  to  provide  a  uniform  intensity  cross  section  at  the 
obervatlon  point.  To  ensure  that  the  center  of  the  laser  beam  fell  upon  the 
camera,  the  outer  edges  of  the  beam  were  visually  located  at  the  observation 
area,  and  the  laser  was  moved  so  that  these  points  were  equidistant  from  the 
camera.  The  beam's  divergence  was  about  20  milliradlans. 

b.  Video  System.  A  Cohu  7210B  camera,  incorporating  an  extended  range 
silicon  Jiode  array  vidlcon,  was  used  to  observe  the  laser  beam.  The  front  end 
optics  Included  a  Celestron  2000-millimeter  f/10  Cassegrain  lens  and  an  adapter 
for  joining  the  lens  and  the  camera.  The  adapter,  manufactured  by  the  NFS, 
consists  of  a  1.55X  Barlow  lens,  a  0. 6328-micrometer  Interference  filter,  and  a 
pair  of  adjustable  polarizers  for  setting  the  illuminance  of  the  vidlcon.  This 
camera  assembly  was  mounted  on  a  tripod  and  positioned  on  a  small  platform  in 
front  of,  and  separate  from,  the  Smoke  Week  platform  to  ensure  vibration  iso¬ 
lation.  The  camera  control  unit  was  used  to  adjust  the  level  of  the  video 
signal  from  the  camera.  Care  had  to  be  taken  that  this  level  was  not  above  the 
clipping  level  of  the  camera  system.  Thus,  the  video  levels  were  continuously 
monitored  with  an  oscilloscope.  The  video  signal  was  recorded  on  a  3/4-inch 
format  Panasonic  video  cassette  recorder.  Time  and  date  annotations  were 
Inserted  on  the  video  signal  prior  to  recording. 

c.  NPS  Slit  Scanning  System.  The  system  consists  of  a  Celestron  8-lnch 
diameter  Schmidt  telescope  with  an  equivalent  focal  length  of  160  Inches  with  a 
times  2  Barlow  lens.  It  has  a  central  obscuration  ratio  of  0.344.  The  tele¬ 
scope  is  equipped  with  up-down  and  right-left  tracking,  and  llnar  scan  across  a 
slit,  for  generating  the  line-spread  function.  The  field  lenses  pass  the  light 
through  a  narrowband  Interference  filter  before  falling  on  the  detector,  which 
is  a  GE  silicon  avalanche  multiplier.  The  auxiliary  equipment  Includes  the 
scan  servos,  amplifiers,  and  on-line  data  processing  equipment.  The  total 
package  occupies  a  volume  of  about  that  of  an  office  desk. 

• 

d.  Equipment  Checkout.  The  preliminary  setup  and  equipment  check-out 
occurred  in  the  morning  of  10  September  1984.  The  measurements  made  on  the  un- 
callbrated  system  were  satisfactory.  The  calibration  was  deferred  until  the 
following  morning  when  the  scintillation  was  likely  to  be  low.  A  multisensor 
meteorological  station  was  located  at  the  test  site  to  provide  meteorological 
data.  A  severe  storm  developed  in  the  early  afternoon  which  terminated  further 
checkout  and  required  the  disassembly  of  the  equipment.  It  had  been  planned  to 
obtain  photo  coverage  of  the  test  site  and  the  lnstrumention.  Unfortunately 
the  adverse  weather  and  the  photographer  showed  up  at  the  same  time  each  day. 


2.2  MEASUREMENTS  OF  11  SEPTEMBER  1984 

On  11  September,  test  setup  was  initiated  at  0600.  The  systems  were 
calibrated  and  the  first  measurements  were  made  at  0748.  The  weather  started 
out  clear  but  by  0800  It  was  partly  cloudy.  It  continued  to  cloud  up  until  It 
was  mostly  cloudy  by  0900.  It  started  to  rain  by  noon  and  the  storm  lasted  all 
afternoon  with  very  heavy  rains.  No  data  were  taken  in  the  afternoon.  A  list¬ 
ing  of  the  data  is  shown  In  Table  1. 


TABLE  I.  DATA  TAKEN  ON  11  SEPTEMBER  1984 


Time 

Conditions 

Cn  (LT) 

Cn  (ST) 

Cn  (Wan) 

0747 

partly  cloudy 

3.036e-07 

— 

— 

0750 

partly  cloudy 

— 

3.795e-07 

— 

0803 

cloudy 

1.733e-07 

— 

— 

0805 

cloudy 

— 

2.152e-07 

1.049e-07 

0824 

cloudy 

2.174e-07 

— 

— 

0827 

cloudy' 

2.880e-07 

— 

— 

0846 

full  sun 

4.488e-07 

— 

— 

0847 

full  sun 

— 

5.125e-07 

2.814e-07 

0859 

partly  cloudy 

3.346e-07 

— 

— 

0900 

partly  cloudy 

— 

4.221e-07 

2.929e-07 

0916 

cloudy 

3.394e-07 

— 

— 

0918 

part  sun 

— 

4.121e-07 

3.316e-07 

1011 

part  sun 

4.121e-07 

— 

- — 

1013 

windy 

— 

4.085e-07 

8.053e-07 

(wind?) 

1044 

cloudy 

2.57le-07 

— 

— 

- 

1045 

cloudy  and  wind 

— 

2.045e-07 

3.535e-07 

1102 

cloudy  and  wind 

1.199e-07 

— 

— 

NO 

1984. 

data  for  the  video 

system  were 

reduced  for  data 

taken  on  11 

September 

2.3.  MEASUREMENTS  OF  12  SEPTEMBER  1984 

INCLUDING  SMOKE  TEST 

On  12  September,  test 

setup  began 

at  0615.  Except 

for  a  delay 

—  the 

generator  at  the  laser  end  of  the  range  ran  out  of  gas  —  everything  was  ready 
to  take  data  by  0730.  The  first  measurements  were  made  just  before  0900.  The 
weather  was  sunny  during  the  morning,  then  became  partly  cloudy  in  the  early 
afternoon,  and  finally,  cloudy  later  in  the  afternoon.  A  crosswind  blew  most 
of  the  afternoon.  The  wind  was  often  strong  enough  to  cause  a  vibration  of  the 
jdatfora  holding  the  telescopes.  This  showed  up  in  the  measurements  of  long 
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term  and  wander.  No  long-term  measurements  were  made  during  the  rest  of  the 
afternoon.  Since  the  wander  measurements  are  very  sensitive  to  vibration,  they 
give  an  indication  of  the  wind.  The  separate  platform  did  Isolate  the  tele-, 
scopes  from  the  vibrations  due  to  people  moving  about  on  the  main  platform,  but 
the  wind  caused  vibrations  which  affected  the  wander  and  the  long-term  measure¬ 
ments.  A  little  rain  fell  between  1245  and  1330,  then  it  started  to  rain  again 
after  1630.  The  smoke  tests  were  made  after  1130.  There  was  no  apparent  ef¬ 
fect  of  the  smoke  from  the  grenades  except  for  increasing  the  extinction. 

After  1600,  black  smoke  from  burning  rubber  was  tested.  This  smoke  definitely 
increased  the  effective  Cq  that  was  measured.  This  was  also  observed  quali¬ 
tatively  with  the  video  system.  Video  recordings  were  made  throughout  the  day 
on  September  12  at  roughly  half-hour  intervals.  Twice  during  the  day,  in  per¬ 
iods  of  low  turbulence,  a  diffraction  grating  was  held  in  front  of  the  camera 
optics  and  recordings  were  made  of  the  laser  spot  diffraction  pattern.  The 
data  provided  scaling  Information  used  in  the  analysis  performed  later.  A 
total  of  about  2-1/2  hours  of  video  data  was  obtained  during  the  2  days  of  data 
collection.  The  data  taken  on  12  September  is  shown  in  Table  II. 

Efforts  were  made  to  operate  the  laser  designator  but  it  behaved  errati¬ 
cally,  hence,  measurements  using  the  laser  designator  were  not  accomplished. 


TABLE  II.  TESTS  PERFORMED  ON  12  SEPTEMBER  1984. 


Time 

Conditions 

Cn  O.T) 

Cn  (ST) 

Cn  (Wan) 

Smoke 

0855 

Clear 

2.481e-07 

— 

— — 

No 

0859 

Clear 

3.346e-07 

— 

— 

No 

0900* 

Clear 

2.828e-07 

2.205e-07 

2.989 

No 

0900 

Sunny 

— 

2.955e-07 

— 

No 

0902 

Sunny 

— 

3.128e-07 

3.465e-07 

No 

0930* 

Sunny 

3.447e-07 

2.976e-07 

8.595e-07 

No 

0931 

Sunny 

3.406e-07 

3.840e-07 

3.326e-07 

No 

0957 

Sunny 

3.981e-07 

3.934e-07 

3.297e-07 

No 

1001 

Sunny 

4.199e-07 

4.300e-07 

3.304e-07 

No 

1001* 

Sunny 

4.127e-07 

2.855e-07 

1.824e-06 

No 

1039 

Sunny 

4.752e-07 

4.553e-07 

4.970e-07 

No 

1039* 

Sunny 

3.801e-07 

3.122e-07 

1.317e-06 

No 

1100 

Some  grenade  set  off  15  m 

upwind  from  beam 

and  114  meters 

down 

range  from  telescopes. 


Indicates  video  data. 


TABLE  II.  TESTS  PERFORMED  ON  12  SEPTEMBER  1984  (CONT) 


Time 

Conditions 

Cn  (IT) 

Cn  (ST) 

Cn  (Wan) 

Smoke 

1106 

4.293e-07 

— 

— 

Yes 

1108 

3.842e-07 

— 

— 

Yes 

1110 

4.458e-07 

— 

— 

Yes 

1110* 

Prtly  Cloudy 

3.494e-07 

3.068e-07 

1.236e-07 

Yes 

1112 

Prtly  Cloudy 

4.458e-07 

— 

— 

Yes 

1113 

Prtly  Cloudy 

4.443e-07 

— 

— 

Yes 

1130 

Cloudy 

1.686e-07 

— 

— 

Yes 

1130* 

Cloudy 

2.518e-07 

2.240e-07 

5.485e-07 

Yes 

1131 

Clearing 

1.936e-07 

— 

— 

Yes 

1132 

Prtly  Cloudy 

2.998e-07 

— 

— 

Yes 

1133 

Sunny 

3.698e-07 

— 

— 

Yes 

1134 

Sunny 

3.590e-07 

— 

— 

Yes 

1135 

Sunny 

4.574e-07 

— 

— 

Yes 

1136 

Prtly  Cloudy 

4.468e-07 

— 

— 

Yes 

1137 

Mstly  Cloudy 

3.983e-07 

— 

— 

Yes 

1138 

Cloudy 

3.835e-07 

— 

— 

Yes 

1139 

Cloudy 

3.007e-07 

— 

— 

Yes 

1140 

Cloudy 

2.401e-07 

— 

— 

Yes 

1141* 

Cloudy 

1.837e-07 

1.441e-07 

6.734e-07 

Yes 

1141 

Cloudy 

1.634e-07 

— 

— 

Yes 

1224 

Cloudy 

2.701e-07 

— 

— 

No 

1225 

Cloudy 

— 

1.297e-07 

1.302e-07 

Jo 

1228 

Cloudy 

9.426e-08 

— 

— 

No 

1229 

Windy  Cloudy 

— 

l.255e-07 

5.060e-07 

No 

There  was  rain  for  45  minutes 

from  1245 

to  1315. 

1340 

Cloudy 

— 

l.547e-07 

2.877e-07 

No 

1348 

Prtly  Cloudy 

3.130e-07 

— 

— 

No 

1349 

Cloudy 

2.988e-07 

— 

— 

No 

Smoke  test  begins  at  114  meters  down  range. 
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TABLE  II.  TESTS  PERFORMED  ON  12  SEPTEMBER  1984  (CONT) 


Time 

Conditions 

Cn  (LT) 

Cn  (ST) 

Cn  (Wan) 

Smoke 

1405 

Cloudy 

2.701e-07 

— 

- — 

Yes 

1406* 

Cloudy 

2.937e-07 

1.063e-07 

9.861e-07 

Yes 

1407 

Cloudy 

— 

3.382e-07 

2.237e-07 

Yes 

1408 

Cloudy 

2.099e-07 

— 

— 

Yes 

1410 

Cloudy 

2.562e-07 

— 

— 

Yes 

1411* 

Cloudy  1.542  e-07 

Telescope  out-of-focus  since 

on  the  NPS  slit-scan  system. 

1.234e-07 

the  rain  - 

4.769e-07 

refocused  the 

Yes 

telescope 

1414 

Cloudy 

9.654e-08 

— 

— 

No 

1416 

Cloudy 

9.737e-08 

— 

— 

Yes 

1417 

Cloudy 

9.547e-08 

— 

— 

Yes 

1418 

Cloudy 

9.454e-08 

— 

— 

Yes 

1419 

Cloudy 

9.643e-08 

— 

— 

Yes 

1421 

Cloudy 

1.040e-08 

— 

— 

Yes 

1422 

Cloudy 

9.498e-08 

— 

— 

No 

1424 

Cloudy 

1.060e-07 

— 

— 

Yes 

1425 

Cloudy 

1.077e-07 

— 

— 

No 

1427 

Cloudy 

9.423e-07 

— 

— 

Yes 

1428 

Cloudy 

— 

i.334e-07 

— 

Yes 

1429 

Cloudy 

— 

1.367e-07 

9.181e-08 

Yes 

1430 

Cloudy 

Smoke  moved  to 

1.605e-07 

800  meters  down  range. 

Yes 

1439 

Overcast 

8.583e-08 

— 

— 

No 

1440 

Overcast 

— 

8.431e-08 

8.916e-08 

No 

1447 

Overcast 

8.950e-08 

— 

— 

No 

1448 

Overcast 

— 

9.442e-08 

1.745e-07 

No 

1454 

Overcast 

— 

1.196e-07 

1.156e-07 

No 

1500 

Overcast 

1.4l2e-07 

— 

— 

Yes 

1500* 

Overcast 

6.902e-08 

4.219e-07 

3.306e-07 

Yes 

1501 

Overcast 

1 .643e-07 

— 

— 

Yes 

1502 

Overcast 

— 

1 . 545e-07 

2.931e-07 

Yes 
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TABLE  II.  TESTS  PERFORMED  ON  12  SEPTEMBER  1984  (CONT) 


Time 

Conditions  Cn  (LT) 

Cn  (ST) 

Cn  (Wan) 

Smoke 

1504 

Overcast 

— 

1.660e-07 

1.397e-07 

Yes 

1505 

Overcast 

— 

1.807e-07 

1.162e-06 

Yes 

There  was  some  vibration 

of  the  platform  during  last  run 

at  1505 

1507 

Overcast 

— - 

1.712e-07 

4.934e-07 

Yes 

1508 

Overcast 

— 

1.704e-07 

2.89le-07 

Yes 

1510 

Overcast 

— 

1.036e-07 

2.733e-07 

Yes 

1512 

Overcast 

— 

1.050e-07 

2.681e-07 

Yes 

1513 

Overcast 

— 

1.041e-07 

1.678e-07 

Yes 

1515 

Thick  Smoke 

— 

1.026e-07 

1.847e-07 

Yes 

1516 

Medium  Smoke 

— 

1.048e-07 

1.913e-07 

Yes 

1519 

Thin  Smoke 

— 

1.012e-07 

1.116e-07 

Yes 

1520 

No  Smoke 

— 

9.361e-08 

1.306e-07 

No 

Moved  smoke  to  498 

meters 

down  range  from 

detectors. 

1529 

Overcast 

— - 

1.037e-07 

3.272e-07 

No 

1531 

Overcast 

- — 

1.060e-07 

1.203e-07 

No 

1541 

Cloudy/Windy 

— 

1.078e-07 

9.332e-07 

Yes 

1546 

Overcast 

— 

1.048e-07 

3.408e-07 

Yes 

1548 

Overcast 

— 

1.066e-07 

4.719e-07 

Yes 

1553 

Overcast 

— 

1.005e-07 

1.953e-06 

Yes 

1555 

Thick  Smoke 

— - 

1 .072e-07 

1.036e-06 

Yes 

1557 

Overcast 

— 

1.169e-07 

2.882e-07 

Yes 

1559 

Overcast 

— 

1.964e-07 

2.059e-07 

Yes 

1600 

Overcast 

- — 

1.060e-07 

1.846e-07 

Yes 

1602 

Thick  Smoke 

— 

1.116e-07 

1.776e-06 

Yes 

1603 

Medium  Smoke 

— 

1.086e-07 

2.307e-07 

Yes 

1604 

Thick  Smoke 

— 

1.146e-07 

1.698e-07 

Yes 

1606 

Light  Smoke 

1.034e-07 

2.261e-07 

Yes 

1607 

Almost  No  Smoke 

— 

1.078e-07 

3.023e-07 

Yes 

1608 

Overcast 

— 

9.577e-08 

9.809e-08 

No 

1610 

Almost  No  Smoke 

— 

1.094e-07 

2.297e-07 

Yes 

1612 

Black  Smoke 

— 

1.043e-07 

8.193e-07 

Yes 

from  Rubber 


TABLE  II.  TESTS  PERFORMED  ON  12  SEPTEMBER  1984  (CONT) 


Tiae  Conditions  Cn  (LT)  Cn  (ST)  Cn  (Wan)  Smoke 


1613* 

Black  Smoke 

5.908e-07 

1614 

Med  Black  Smoke 

— 

1614* 

Some  Smoke 

9.526e-08 

1615 

Very  Light 

— 

Black  Smoke 

1617 

Almost  No  Smoke 

— 

1619 

Threatening 

— 

Rain 


Negative 

2.476e-07 

Yes 

1.721e-07 

2.728e-07 

Yes 

Negative 

5.678e-08 

Yes 

1.076e-07 

8.331e-07 

Yes 

1.138e-07 

8.414e-08 

Yes 

1.038e-07 

1.47le-07 

No 

The  effect  of  the  red  phosphorus  smoke  was  small  and  noticeable  only  when 
the  sky  was  overcast  and  very  little  turbulence.  The  change  in  the  Cn?  was 
only  a  3  to  10  percent  Increase.  The  black  smoke  from  the  burning  rubber 
caused  about  60  percent  increase  in  the  Cn2.  The  Increase  in  the  Cn2 
is  due  to  the  heat  of  combustion  causing  the  region  of  the  smoke  to  have  a 
larger  temperature  variation  than  the  nearby  air  without  smoke.  The  burning 
rubber  gave  off  considerably  more  heat  than  the  smoke  grenade,  hence,  there  was 
much  more  effect  due  to  the  smoke  from  the  rubber  than  the  phosphorus.  If  it 
had  not  been  overcast  this  effect  probably  would  have  been  hidden  by  the  higher 
turbulence  when  the  sun  was  shining. 

2.4  VIDEO  DATA  ANALYSIS 

a.  The  video  data  obtained  in  the  field  was  analyzed  with  the  equipment 
setup  shown  in  figure  2.  The  video  cassette  recordings  were  first  dubbed  onto 
an  Ampex  VPR-80  1-inch  reel-to-reel  video  tape  recorder  to  make  use  of  some  of 
the  Ampex  recorder's  special  features,  such  as  the  ability  to  view  single  video 
fields.  The  time  base  corrector  was  an  Ampex  TBC-80,  designed  to  enable  the 
Ampex  recorder  to  perform  special  video  effects.  The  time  base  corrector  also 
allowed  for  adjustment  of  the  video  signal  levels  and  for  the  video  black 
level. 


b.  Each  field  of  video  to  be  analyzed  was  digitized  into  the  memory  of  a 
Quantex  DS-30  video  processor.  The  field  of  video  is  stored  as  512  x  512  array 
of  picture  elements  (pixels)  represented  by  a  gray  scale  value  between  0 
(black)  and  255  (white).  These  values  could  be  randomly  accessed  and  transfer¬ 
red  via  the  Hewlett-Packard  Interface  bus  (HPIB)  to  an  HP  1000  minicomputer. 

c.  Before  a  processing  run  was  made,  a  typical  field  of  video  was  digi¬ 
tized.  The  digitized  values  were  read  into  the  computer  using  the  program 
WINDOW,  developed  by  the  Electromagnetic  Environmental  Test  Facility  (EMETF), 
which  helped  locate  the  position  of  the  laser  image  in  the  field.  Only  data 
without  a  256-  by  256-pixel  window,  centered  on  the  laser  image,  were  sub¬ 
sequently  processed.  WINDOW  also  enabled  the  scaling  parameter,  SCALE,  to  be 
determined.  SCALE  was  determined  by  first  digitizing  a  field  of  video  acquired 
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when  a  diffraction  grating  was  in  front  of  the  camera  optics.  WINDOW  was  used 
to  count  the  number  of  pixels  between  the  central  and  first  maxima  of  the  dif¬ 
fraction  pattern.  The  calculated  angular  separation  in  microradians  of  the 
maxima  was  divided  by  their  pixel  separation  to  derive  the  value  of  SCALE. 


d.  The  actual  processing  of  the  digitized  video  data  was  performed  by  a 
FORTRAN  Program  called  MTFLTV ,  which  was  supplied  by  NPS.  This  program  has  as 
inputs  the  digitized  video  data  as  well  as  test  parameters  such  as  the  distance 
from  the  laser  to  the  camera  and  specifications  of  the  front-end  optics.  It 
outputs  the  optics  MTF,  the  atmospheric  MTF,  and  the  turbulence  structure. 

Ca2- 


2.5  VIDEO  DATA  RESULTS 


a.  Data  from  12  different  time  periods  on  September  12  were  processed. 
The  Cn2  values  calculated  for  long-term,  short-term,  and  wander  data  are 
shown  separately  in  table  III.  This  table  also  indicates  the  start  time  of 
each  processing  run  and  its  duration  in  seconds.  For  each  run,  60  separate 
video  fields  were  processed.  These  fields  were  evenly  distributed  within  the 
run's  duration  time;  so  several  fields  were  not  processed. 


TABLE  III.  CALCULATED  ATMOSPHERIC  STRUCTURE  PARAMETERS 


Run 

No .  Time 


Duration 

(3) 


C^2  (Short-Term)  Cn2  (Long-Term)  Cn2  (Wander) 


1 

9:00:00 

6 

4.862 

X 

10“14 

8.000 

X 

10-14 

8.936 

X 

10-13 

2 

9:30:50 

12 

8.8S8 

X 

10-14 

1.188 

X 

10-13 

7.387 

X 

10-13 

3 

10:00:50 

12 

8.153 

X 

10-14 

1.703 

X 

10-13 

3.328 

X 

10-12 

4 

10:39:00 

12 

9.747 

X 

10-14 

1.445 

X 

10-13 

1.735 

X 

10-12 

5 

11:10:00 

12 

9.412 

X 

10-14 

1.221 

X 

10-13 

1.528 

X 

10-12 

6 

11:30:00 

12 

5.018 

X 

10-14 

6.339 

X 

10-14 

3.008 

X 

10-13 

7 

11:40:48 

12 

2.077 

X 

10-14 

3.733 

X 

10-14 

4.535 

X 

10-13 

8 

14:06:00 

12 

1.130 

X 

10-14 

4.148 

X 

10-14 

9.724 

X 

10-13 

9 

14:11:00 

12 

1.522 

X 

10-14 

2.378 

X 

10-14 

2.274 

X 

10-13 

10 

15:00:00 

12 

1.780 

X 

10-15 

4.764 

X 

10-15 

1.028 

X 

10-13 

11 

16:12:53.30 

6 

-3.320 

X 

10-16 

3.490 

X 

10-15 

6.133 

X 

10-14 

12 

16:13:35 

6 

■  -3.711 

X 

10-16 

9.074 

X 

10-15 

3.224 

X 

10-13 

( 


b.  The  short-term  Cn2  for  runs  11  and  12  turned  out  to  be  m.’ative, 
which  is,  of  course,  an  Invalid  result.  The  specific  reasons  for  this  result 
are  unclear:  however,  In  both  of  these  runs,  the  spot  size  was  small  and  smoke 
was  present  in  the  optical  path. 

c.  It  was  discovered  that  the  Quantex  video  digitizer  tended  to  clip  high 
level  video  signals;  therefore,  the  laser  image  had  an  intensity  platea.  in  its 
center.  When  this  plateau  occurred,  the  processing  program  produced  negative 
values  of  Cn2  for  both  the  short  term  and  long  term.  The  problem  was  cor¬ 
rected  by  lowering  the  video  output  level  from  the  time  base  corrector  until  it 
was  just  below  the  clipping  point*  It  was  also  discovered  that  the  program 
produced  negative  Cn2  values  if  the  background  video  level  was  high.  This 
problem  was  corrected  by  adjusting  the  black  level  on  the  time  base  corrector 
until  all  the  background  was  at  or  below  the  black  video  level. 


2-10 


(THIS  PAGE  INTENTIONALLY  LEFT  BLANK) 


APPENDIX  A.  VIDEO  TEST  DATA  ACQUISITION  FOR  ATMOSPHERIC  MTF  MEASUREMENT 

a.  This  appendix  outlines  the  general  procedures  for  obtaining  the  video 
images  used  as  inputs  to  programs  MTFLTV  and  MTFLD.  The  procedures  discussed 
here  are  not  complete  and  do  not  require  strict  adherence. 

b.  The  equipment  essential  to  obtaining  the  images  consists  of: 

(1)  a  laser 

(2)  a  video  camera 

(3)  a  telescopic  lens  adapted  to  the  camera 

(4)  a  filter  which  passes  light  only  at  the  laser  frequency 

(5)  a  video  recorder 

(6)  a  time  annotator 

(7)  a  video  monitor 

(8)  an  oscilloscope 

(9)  a  large  white  target  board. 

The  camera  system  and  the  laser  system  must  be  free  from  vibrations. 

c.  To  obtain  the  images  needed  by  MTFLTV,  the  laser  (e.g.,  a  low-power 
helium-neon  laser)  is  set  up  a  short  distance  behind  the  target  board.  The 
camera  and  its  recording  equipment  are  set  up  at  least  one  kilometer  in  front 
of  the  target  board.  The  camera  is  aimed  at  the  laser  which  radiates  through 

a  small  hole  in  the  center  of  the  board.  The  camera's  field  of  view  is  limited 
to  the  surface  of  the  target  board  by  the  telescopic  lens.  Thus,  the  desired 
image  recorded  from  the  camera  is  a  bright  spot  (the  laser  beam)  in  the  center 
of  a  dark  background  (the  target  board).  The  images  arc  recorded  by  a  video 
recorder  such  as  a  VCR  for  later  processing  by  MTFLTV. 

d.  Care  must  be  taken  to  insure  the  laser  beam  does  not  overdrive  the 
camera's  amplifiers,  causing  the  spot  image  to  be  clipped.  Clipping  may  be 
detected  by  connecting  the  oscilloscope  to  the  camera  video  outputs.  To  re¬ 
duce  the  intensity  of  the  beam  at  the  camera,  a  diverging  lens  may  be  placed 
in  front  of  the  laser. 

e.  Time  annotation  should  be  placed  on  the  video  so  that  the  results  of 
data  processing  can  he  related  to  the  time  of  day.  A  log  should  be  kept  of 
the  times  when  measurements  are  taken  and  of  the  current  weather  conditions. 

f .  At  some  point  of  low  atmospheric  turbulence,  a  large  diffraction 
grating  should  be  placed  in  front  of  the  camera  optics  and  a  video  recording 
made  of  the  resulting  diffraction  pattern.  The  orientation  of  the  diffractidn 
grating  must  be  parallel  to  that  of  the  video  scan  direction.  This  recording 
through  the  grating  is  necessary  to  determine  the  scaling  parameter  (SCALE) 
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which  is  input  to  MTFLTV.  Also,  the  camera- to- laser  distance  must  be  measured 
to  provide  the  ranging  parameter  (RANGE) .  If  the  telescopic  lens  is  of  the 
Cassegrainian  type,  then  its  obscuration  ratio  (OBS)  must  be  determined  by  tak¬ 
ing  the  ratio  of  the  diameters  of  the  central  obscuration  ratio  and  the  whole 
lens.  Finally,  the  diameter  of  the  lens  must  be  measured  to  provide  the  optics 
diameter  (DOPT). 

g.  Because  MTFLD  was  designed  for  use  with  laser  designator  beams,  its 
equipment  setup  differs  somewhat  from  that  of  MTFLTV.  Laser  designator  beams 
are  generated  by  high-power  lasers,  which  usually  cannot  be  viewed  directly 
through  cameras.  Thus,  for  MTFLD,  the  camera  and  the  laser  both  face  the 
target  board  from  the  same  side.  The  camera  is  set  up  close  to  the  board  so 
that  the  board  fills  its  field  of  view  without  the  telescopic  lens.  The  laser 
is  set  up  at  least  one  kilometer  away  from  the  board.  Recordings  are  then  made 
of  the  laser  spot  that  appears  on  the  target  board.  The  rest  of  the  setup  and 
measurements  are  the  same  as  for  MTFLTV. 


AFPENDIX  B  -  PART  1.  HARDWARE  AND  SOFTWARE  SYSTEMS  FOR  VIDEO 
“  DATA  ANALYSIS 


1 .  INTRODUCTION 

1.1  Test  Background 

USAEPG  and  the  Naval  Postgraduate  School  (NPS)  have  jointly  developed  a 
video  data  technique  for  making  a  path-weighted  measurement  of  a  modulation 
transfer  function  (MTF) .  (See  part  2  of  this  appendix  for  a  general  discussion 
of  MTF  calculation  theory).  The  feasibility  of  the  technique  has  been  demonstr¬ 
ated  in  a  laboratory  environment  by  comparing  the  results  from  the  video  system 
with  those  from  the  NPS  slit-scanning-telescope  system.  To  validate  the 
laboratory  comparison  of  the  two  techniques,  and  to  validate  the  laboratory 
tests  themselves  as  valid  measures  of  field  performance,  a  joint  field  test  was 
run  in  September  1984  at  Fort  Huachuca,  Arizona.  This  appendix  describes  hard¬ 
ware  and  software  systems  used  to  analyze  the  video  data  taken  during  the  field 
test. 


1.2  Software  Considerations 

The  NPS  developed  software  consisting  of  three  main  programs  (MTFLTV , 
MTFLD,  and  PREDICT)  and  several  subroutine  and  function  modules.  The  develop¬ 
ment  was  not  done  on  the  present  hardware;  however,  prior  to  delivery,  the  NPS 
modified  the  software  for  use  on  an  Hewlett-Packard  (HP)  HP-1000  minicomputer. 
Upon  delivery,  the  modified  software  had  not  been  tested  or  verified  on  an  HP- 
1000.  (At  the  time  of  the  writing  of  this  document,  MTFLD  and  PREDICT  remain 
untested  and  unverified.)  In  addition,  the  software  delivered  was  originally 
written  for  a  much  larger  and  faster  computer,  where  efficiency  was  not  so  im¬ 
portant  as  on  the  HP-1000;  processing  which  took  a  few  minutes  on  the  original 
computer  sometimes  requires  over  an  hour  on  the  HP-1000.  Those  inefficiencies 
which  could  be  remedied  by  additional  code  modifications  are  addressed  in  this 
document . 

1.3  Required  Peripherals 

a.  Two  pieces  of  equipment  other  than  the  HP-1000  computer  are  essential 
for  running  MTFLTV  and  MTFLD.  The  first  is  a  video  digitizer  with  computer- 
accessible  memory.  The  second  is  a  video  recorder  capable  of  stepping  through 
individual  video  fields. 

b.  The  video  digitizer  used  as  the  basis  for  the  present  versions  of 
MTFLTV  and  MTFLD  is  the  Quantex  DS-30  Video  Processor.  An  effort  has  been 
made  to  isolate  the  code  relating  to  the  Quantex  so  that  the  programs  can  be 
easily  modified  for  use  with  another  digitizer.  The  modules  affected  are 
VIDEODATA  and  DIGITIZE. 

c.  Any  video  recorder  which  can  be  stepped  through  individual  fields  may 
be  used,  e.g.,  an  Eigen  videodisc  unit  or  an  Ampex  VPR  80  video  tape  recorder. 
The  programs  currently  do  not  depend  on  the  model  of  video  recorder  used  be¬ 
cause  the  recorder  is  assumed  to  1 _  manually  operated  and  is  linked  to  the 
digitizer  only  through  a  video  cable. 


1.4  Device  Efficiency 

a.  To  increase  efficiency,  the  Quantex  DS-30  bypasses  the  FORTRAN  for¬ 
matter  and  interfaces  the  program  through  EXEC  calls  (operating  system  level 
commands).  The  EXEC  calls  are  probably  the  most  system-dependent  part  of  the 
programs  and  should  be  kept  in  mind  when  transporting  the  code  to  other  operat¬ 
ing  environments. 

b.  The 'manual  operation  of  the  video  recorder  is  very  time-consuming  for 
the  program  user.  Ideally,  the  code  should  be  modified  so  that  the  programs 
control  both  the  video  recorder  and  the  digitizer.  These  modifications  would 
be  made  in  the  D INPUT  module. 

2.  GENERAL  PROGRAM  DESCRIPTION 

a.  The  MIF  software  contains  three  main  programs:  MTFLTV,  MTFLD,  and 
PREDICT.  MTFLTV  and  MTFLD  are  very  similar  and  use  many  of  the  same  program 
modules.  Both  use  as  main  Inputs  a  series  of  digitized  video  images  which, 
in  combination  with  other  inputs  such  as  scale  factors  and  the  distance  from 
the  laser  to  the  target,  are  used  to  calculate  the  atmospheric  MTF.  The  MTF 
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is  then  used  to  calculate  C  ,  a  measure  of  the  atmospheric  turbulence  affect¬ 
ing  image  quality. 

b.  The  major  difference  between  MTFLTV  and  MTFLD  field  test  methods  is 
MTFLTV  uses  images  of  a  helium-neon  laser  beam  directed  at  a  camera,  whereas 
MTFLD  uses  images  of  the  reflection  of  a  laser  designator  beam  off  a  white 
target  board.  Both  tests  place  the  laser  approximately  one  kilometer  away 
from  the  target. 
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c.  The  third  program,  PREDICT,  uses  a  previously  calculated  value  of  C 

n 

to  reconstruct  the  atmospheric  MTF.  PREDICT  then  uses  the  atmospheric,  source, 
and  optics  MTF's  to  predict  the  image  quality  of  a  given  optical  system. 

d.  Figures  1-3  show  the  logical  Interconnections  between  the  modules  in 
MTFLTV,  MTFLD,  and  PREDICT.  Details  of  the  separate  modules,  including  their 
functions,  inputs/ outputs,  and  source  code  listings,  are  discussed  in  appendix 
C.  The  programs  are  structured  vertically,  l.e.  most  of  the  modules  have  a 
direct  connection  only  with  the  main  routine  and  not  with  each  other.  Thus, 
the  program  flow  is  fairly  simple. and  straightforward.  However,  many  of  the 
modules  do  interface  through  common  data  blocks,  and  this  fact  must  be  kept  in 
mind  when  program  modifications  are  made. 

e.  For  MTFLTV  and  MTFLD,  all  the  input  data  is  provided  in  modules  GINPUT 
and  DINPUT.  GINPUT  allows  for  the  input  of  test-specific  parameters  via  an  in¬ 
teractive  terminal  or  a  data  file.  DINPUT  preprocesses  the  digitized  video 
data  and  stores  the  results  in  arrays.  The  data  in  these  arrays  is  then  proc¬ 
essed  further  by  other  modules  in  MTFLTV  and  MTFLD.  Each  time  the  data  is  proc¬ 
essed,  the  new  results  overwrite  the  previous  data.  Intermediate  results  may 
be  output  to  data  files  using  the  OUTPUT  module.  The  ultimate  values,  the 
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Cn  parameters,  are  displayed  on  the  operator's  terminal  as  well  as  being  out¬ 
put  to  a  data  file. 


(Text  continued  ontpage  B-6) 


Figure  1.  MTFLTV  Module  Interconnections 


Figure  3.  PREDICT  Module  Interconnections 
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f.  For  PREDICT,  Che  module  GINPUT  performs  all  of  Che  inpuC.  Again, 
many  of  Che  lncermediace  calculacions  as  well  as  Che  final  resulcs  may  be  out- 
puc  co  daca  files.  This  ouCpuC  is  performed  by  che  POUTPUT  module. 

3.  RUNNING  THE  PROGRAMS 

3.1  Hardware 

a.  The  programs  can  be  run.  on  HP-1000  series  compucers  under  Che  RTE-6/ 

VM  operating  system.  They  have  also  been  run  on  Che  HP  A-600  computer,  under 
the  RTE-A  operating  system.  All  program  modules  are  coded  in  HP's  FORTRAN  7X. 
The  programs  currently  require  that  a  Quantex  DS-30  Video  Processor  be  linked 
Co  the  computer  via  the  HP  Interface  Bus  (HPIB). 

b.  The  programs  MTFLTV  and  MTFLD  are  especially  large  and  barely  fic  into 
a  64  kilobyte  memory  partition.  When  modifications  are  made  to  che  program,  it 
may  be  necessary  to  use  a  larger  partition,  reduce  che  sizes  of  che  arrays,  or 
put  the  arrays  into  che  Extended  Memory  Areas  (EMA) . 

3.2  Program  Inputs 

The  general  inpuc  parameters  which  encer  a  program  via  module  GINPUT  may 
be  input  interactively  at  the  operator's  terminal  or  from  the  file  INPUT.DAT. 

In  either  case,  these  parameters  must  be  determined  by  the  operator  prior*  to 
running  a  program.  Below  is  a  brief  description  of  each  input  parameter  and 
how  its  value  is  determined.  [NOTE:  LTV  denotes  use  in  MTFLTV;  LD  in  MTFLD; 

P  in  PREDICT.] 

(1)  A  (real).  The  maximum  intensity  of  the  laser  source  in  watts/ 
steradian.  (P) 

(2)  DATE  (character *80) .  Time  and  date  of  the  run;  placed  in  the 
header  of  the  output  files.  (LTV,  LD) 

(3)  DOPT  (real) .  The  diameter  of  the  objective  of  the  optics  in 
meters;  determined  by  measuring  across  che  objective.  (LTV,  LD,  P) 

(4)  EXT  (real) .  The  extinction  coefficient  of  the  atmosphere  in 
meters  ^  (P) . 

(5)  HEADER1  (character*80) .  Contents  about  the  run;  placed  in  the 
header  of  the  output  files  (LTV,  LD).  Information  about  the  laser  source.  (P) 

(6)  BEADER2  (character*80) .  Comments  about  the  weather;  placed  in 
the  header  of  the  output  files  (LTV,  LD).  Information  about  the  optics.  (P) 


(7)  10  (integer) .  Parameter  specifying  whether  or  not  all  inter¬ 
mediate  calculations  are  to  be  output.  If  0,  then  no  output  occurs;  if  1, 
then  all  intermediate  calculations  are  to  be  output.  (LTV,  LD,  P) 

(8)  LOPTICS  (integer).  Parameter  specifying  the  type  of  optics 
diffraction  function  to  be  used  by  the  program.  If  1,  then  the  function  is 
computed;  if  2,  it  is  input  from  a  file.  (P) 


(9)  LSOURCE  (integer) .  Parameter  specifying  the  type  of  laser 
source  function  to  be  used  by  the  program.  If  1,  then  a  Gaussian  source  func¬ 
tion  is  computed;  if  2,  then  the  source  function  is  input  from  a  file;  if  3, 
then  the  source-optics  function  combination  is  input  from  a  file.  (P) 

(10)  NFLDS  (integer) .  The  number  of  fields  of  video  data  to  be 
processed.  Note  that  the  fields  may  be  separated  by  any  fixed  time  interval. 
(LTV,  LD) 

(11)  OBS  (real).  The  ratio  of  the  diameter  of  the  central  obscura¬ 
tion  of  the  objective  to  the  diameter  of  the  objective;  determined  by  measur¬ 
ing  across  the  obscuration  and  dividing  by  DOPT.  (LTV,  LD,  P) 

(12)  RANGE  (real).  The  distance  between  the  optics  and  the  target; 
determined  by  a  surveyor's  range  finder.  (LTV,  LD,  P) 

(13)  SCALE  (real) .  Angular  scale  of  a  digitized  video  image  as 
specified  by  the  angle  subtended  by  one  pixel;  in  microradians  per  pixel.  De¬ 
termined  by  digitizing  a  video  field  acquired  through  a  diffraction  grating, 
calculating  the  angular  separation  of  the  central  and  first  maxima  of  the  laser 
image  diffraction  pattern  (a  function  of  the  diffraction  grating),  counting  the 
number  of  pixels  between  these  maxima  in  the  .digitized  image,  and  finally  tak¬ 
ing  the  ratio  of  the  calculated  angular  separation  of  the  maxima  and  the  number 
of  pixels  between  the  maxima.  (LTV,  LD,  P) 
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(14)  SQCN  (real).  The  previously-determined  value  of  C  for  the 

-2/3  n 

atmosphere  in  units  of  meters  .  (P) 


(15)  SSIGMA  (real).  The  divergence  of  the  laser  source,  i.e.  the 
angle  between  the  beam  center  and  the  line  at  which  the  laser's  energy  falls 
to  1/e  of  its  peak  energy;  in  radians.  This  value  is  a  characteristic  of  the 
laser  used;  usually  determined  in  the  laboratory.  (LD,  P) 

(16)  WAV  (real).  Wavelength  of  the  laser  in  micrometers;  deter¬ 
mined  by  the  type  of  laser  used,  e.g.,  a  helium/ neon  laser  radiates  at 
0.6328  pm.  (LTV,  LD,  P) 

(17)  XW1,  XW2,  YW2,  YW1  (integer).  The  left,  right,  top,  and  bot¬ 
tom  locations  of  the  digitized  video  window,  which  is  centered  on  the  laser 
spot  image.  These  locations  are  specified  in  pixel  units  of  the  digitizer. 
Care  must  be  taken  to  insure  that  the  window  is  large  enough  to  completely 
enclose  the  laser  spot  images;  however,  making  it  unnecessarily  large  signi¬ 
ficantly  increases  processing  time.  Currently,  the  top  and  bottom  of  the 
window  must  be  separated  by  either  256  or  512  pixels.  (LTV,  LD) 

3.3  Compiling  and  Linking  the  Programs 

a.  The  programs  MTFLTV,  MTFLD,  and  PREDICT  consist  of  several  modules, 
each  of  which  exists  in  a  separate  file.  A  command  file  has  been  created 
which  automatically  compiles  and  links  all  of  these  modules.  The  file,  named 
COMP&LINK.CMD,  is  listed  in  appendix  D.  To  utilize  the  file,  one  simply 
enters  "TR  COMP&LINK.CMD”  at  a  terminal,  and  all  of  the  command  lines  in  the 
file  are  executed  sequentially. 


b.  COMP &L INK . CMD  does  the  following: 

(1)  The  operator  must  verify  that  the  list  of  source  code  file  names 
in  MERGE.CMD  is  up-to-date 

(2)  All  source  code  files  (format:  NAME.FTN)  listed  in  MERGE.CMD 
are  compiled 

(3)  The  extensions  of  the  compiled  versions  are  changed  to  .REL  (for 
“relocatable") 

(4)  MERGE  concatenates  all  .REL  files  into  one  file,  MTF.REL 

(5)  LINDX  creates  MTF.LIB,  an  indexed  version  of  MTF.REL 

(6)  To  save  disk  space,  all  .REL  files  are  purged 

(7)  LINK  links  MTFLTV ,  MTFLD,  and  PREDICT  into  the  system. 

c.  Control  is  then  returned  to  the  operator.  The  programs  are  ready  to 
be  run.  Note  that  MTF.LIB  is  not  purged  and  may  be  accessed  to  link  with  other 
programs  using  MTF  modules. 

3.4  MTFLTV  Run  Procedure 

At  present  MTFLTV  and  MTFLD  require  a  significant  amount  of  operator  ef¬ 
fort  to  run  successfully.  The  following  set  of  step-by-step  procedures  is 
presented  to  help  the  operator  properly  setup  and  run  MTFLTV  and  MTFLD. 

(1)  Set  up  the  processing  equipment  as  shown  in  figure  4.  Note 
Chat  all  connections  are  made  using  75S2  video  cable  except  the  connection  to 
the  computer,  which  uses  an  HPIB  cable. 

(2)  Make  sure  that  the  logical  unit  numbers  of  the  HPIB  card  and 
video  digitizer  correspond  to  those  in  variable  LUHPIB  and  LUDS30  of  modules 
VIDEODATA  and  DIGITIZE.  If  they  do  not,  then  change  the  values  of  LUHPIB  and 
LUDS30.  Programs  MTFLTV  and  MTFLD  must  then  be  recompiled  and  relinked. 

(3)  Make  sure  Chat  the  program  is  ready  to  run  on  the  system  (see 
para  3.3).  . 

(4)  Set  the  working  directory  to  /MTF/. 

(5)  Make  sure  that  there  is  enough  room  in  directory  /MTF/  to  store 
all  the  program  output  if  Che  output  is  to  be  stored.  If  data  from  previous 
runs  is  to  be  saved,  make  sure  the  data  files  have  been  renamed. 

(6)  With  the  stop-field  video  recorder  (in  this  case,  an  Ampex  VPR- 
80),  locate  the  start  field  of  the  video  data  to  be  processed.  For  the  time  in¬ 
terval  over  which  the  video  data  is  to  be  processed,  visually  check  the  picture 
for  gross  problems,  such  as  picture  break-up,  which  might  affect  the  results. 

If  the  video  signal  is  bad,  then  the  data  may  not  be  reducible.  Also,  monitor 


the  oscilloscope  for  any  signs  of  signal  clipping,  which  causes  the  spot  inten¬ 
sity  profile  on  the  oscilloscope  to  have  a  plateau.  Clipping  can  usually  be 
alleviated  by  lowering  the  video  signal  level  out  of  the  time  base  corrector. 

If  the  signal  is  clipped  on  the  recording,  then  that  video  data  is  not  reduc¬ 
ible. 


(7)  Digitize  one  field  of  video  with  the  video  digitizer.  Make  sure 
that  the  digitized  video  appears  the  same  as  the  original  video  and  has  not 
been  amplified  or  attenuated  in  any  way. 

(8)  Determine  the  values  of  the  general  input  parameters  for  the 
program  (see  para  3.2).  If  these  are  to  be  input  from  a  data  file  by  the 
program,  then  enter  the  values  into  file  INPUT.DAT. 

(9)  Again,  locate  the  start  field  of  the  video  data  to  be  processed. 

(10)  Run  the  program. 

(11)  The  program  asks  if  the  general  input  parameters  are  to  be 
entered  at  the  operator's  terminal.  Entering  ''l”  at  the  terminal  initiates  the 
interactive  input  routine.  Entering  any  other  number  causes  the  parameters  to 
be  read  from  the  input  file.  The  Interactive  input  routine  asks  the  operator 
to  enter  the  parameters  one-by-one  at  the  terminal. 

(12) '  Once  all  of  the  general  input  parameters  are  entered,  the 
program  sounds  the  beep  on  the  operator's  terminal  and  displays  the  message: 

"Advance  to  field  #1  out  of  [NFLDS],'.' 

where  [NFLDS]  is  the  number  of  fields  of  video  to  be  processed.  Be  sure  the 
stop-field  recorder  is  displaying  the  correct  field.  Once  it  is,  enter  "1" 
at  the  terminal.  The  program  inputs  the  video  data  from  the  digitizer  and  pre- 
processes  it.  At  this  point  the  program  is  operating  in  module  DINPUT.  Once 
the  current  video  field  is  preprocessed,  the  program  loops  back  to  the  message 
which  tells  the  operator  to  advance  to  .the  next  field.  This  loop  continues 
until  all  the  fields  have  beetr  digitized  and  the  data  preprocessed.  Note 
that  it  is  not  necessary  to  process  every  field  of  video  in  the  time  interval 
of  Interest.  However,  the  time  separation  between  fields  must  remain  fixed. 

For  example,  suppose  ten  seconds  of  video  data  is  to  be  processed.  This  ten 
seconds  contains  600  fields.  If  the  operator  wants  to  only  process  60  fields 
in  total,  then  every  tenth  field  should  be  processed. 

(13)  After  the  video  data  has  been  digitized  and  preprocessed,  the 
program  processes  the  data  further  without  operator  intervention.  The  values 
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of  the  long  term,  short  term,  and  wander  Cn  are  output  both  to  the  terminal 

and  to  the  file  0UTPUTCN.DAT.  If  the  full  output  option  was  chosen  (10  set  to 
one),  then  the  intermediate  calculations  of  the  FSF,  LSF,  atmospheric  MTF,  etc., 
are  output  to  files  OUTPUT [nn] .DAT,  where  [nn]  is  a  different  number  for  each 
type  of  output. 

(14)  The  program  terminates.  Each  run's  output  files  must  be  renamed 
if  they  are  to  be  protected  from  overwriting. 


3.5  PREDICT  Run  Procedure 


The  procedures  for  running  PREDICT  are  straightforward. 


3.3). 


(1)  Make  sure  the  program  is  ready  to  run  on  the  system  (see  para 


(2)  Set  the  working  directory  to  /MTF/. 


(3)  Make  sure  that  directory  /MTF/  has  enough  room  to  store  the  pro¬ 
gram  outputs.  If  previous  data  is  to  be  saved,  make  sure  the  data  files  have 
been  renamed  so  it  is  not  overwritten. 


(4)  Determine  the  values  of  the  general  input  parameters  for  the 
program  (see  para  4.2).  If  these  are  to  be  input  from  a  data  file  by  the  pro¬ 
gram,  then  enter  the  values  into  file  INPUT.DAT. 

(5)  Run  the  program. 

(6)  The  program  asks  if  the  general  input  parameters  are  to  be 
entered  at  the  operator's  terminal.  Entering  “1”  at  the  terminal  initiates 
the  interactive  input  routine.  Entering? any  other  number  causes  the  parameters 
to  be  read  from  the  input  file.  The  interactive  input  routine  asks  the  opera¬ 
tor  to  enter  the  parameters  one-by-one  at  the  terminal. 

(7)  The  program  then  runs  without  further  operator  intervention. 

The  program  outputs  are  written  into  files  OUTPUT [nn] .DAT,  where  [nn]  is  a 
different  number  for  each  type  of  output. 

(8)  The  program  terminates.  Each  run's  output  files  must  be  renamed 
if  they  are  to  be  protected  from  overwriting. 


APPENDIX  B,  PART  2.  THEORY  OF  MTF  CALCULATIONS 


a.  This  appendix  presents  some  of  the  background  theory  behind  the  cal¬ 
culations  performed  by  the  subject  programs.  This  discussion  is  not  meant  to 
be  rigorous;  rather,  the  reader  is  advised  to  consult  the  references  listed 
in  appendix  E  for  more  complete  expositions. 

b.  Anything  which  lies  between  object  (source)  and  image  (response)  af¬ 
fects  the  linage  quality  of  an  optical  system.  Optical  systems  themselves  de¬ 
grade  images  from  light  sources  due  to  effects  such  as  diffraction  and  aber¬ 
rations.  The  planar  description  of  this  image  degradation  is  the  point  spread 
function  (PSF).  The  one-dimensional  (single  variable)  analog  of  the  PSF  is  the 
line  spread  function  (LSF),  the  measure  of  an  optical  system's  response  to  a 
slit  source  of  light.  The  LSF  can  be  generated  from  the  PSF  by  integrating 
over  one  of  the  two  dimensions  of  a  point  source  image.  The  MTF  is  the  real 
part  of  the  Fourier  transform  of  the  PSF  or  of  the  LSF.  MTF's  are  often  used 
when  analyzing  image  quality,  rather  than  PSF's  or  LSF's,  because  the  overall 
MTF  (and  thus  the  image  quality)  of  a  complex  optical  system  can  be  simply  de¬ 
termined  by  multiplying  the  MTF's  of  each  component  of  the  system. 

c.  An  MTF  can  be  determined  for  any  source  of  image  degradation  which 
can  be  isolated,  including  lenses,  cameras,  imaging  sensors,  and  the  atmos¬ 
phere.  The  image  degradations  caused  by  the  atmosphere  are  a  function  of  den¬ 
sity  fluctuations,  and  can  be  described  by  the  index  of  refraction  structure 
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parameter,  C  .  C  may  be  determined  statistically  by  relating  distance  to 
ii  n  ^ 

the  fluctuation  amplitude  of  the  index  of  refraction.  Alternatively,  may 

be  approximated  using  Fried's  formula  by  performing  a  least  squares  fit  of  the 
measured  atmospheric  MTF: 

"Cn2* 

MTF(x)  -  Ae 

The  second  method  is  used  by  MTFLTV  and  MTFLD.  On  the  other  hand,  PREDICT 

2 

uses  Fried's  formulas  and  a  value  of  C  (either  from  the  programs  above  or 

n 

from  manual  input)  to  calculate  the  atmospheric  MTF. 
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APPENDIX  C.  MODULE  DESCRIPTIONS  AND  SAMPLE  DATA  PRINTOUTS  FOR  THE 

MTF  ANALYSIS  COMPUTER  PROGRAM 

Part  1  of  this  appendix  comprises  module  information  sheets,  each  containing 
a  brief  description  of  a  module,  a  list  of  which  modules  reference  it,  a  list  a 
data  interfaces  with  other  modules,  and  notes  which  a  programmer  might  find  useful 
when  maintaining  and  modifying  the  code.  Note  that  the  following  compiler  instruc¬ 
tions  are  specific  to  present  compiler:  "FTN7X,L,Q"  at  the  beginning  of  each  module 
listing,  and  ?FIL£S(0,2)"  in  the  main  program  listings. 

Part  2  of  this  appendix  contains  sample  data  printouts  for  a  representative  run 
(#3)  processed  by  the  program  MTFLTV.  The  following  sets  of  data  are  included: 

1.  Long-term  raw  data 

2.  Short-term  raw  data 

3.  Long-term  LSFT  adjusted  data 

4.  Short-term  LSFT  adjusted  data 

5.  FFT  of  long-term  data 

6.  FFT  of  short-term  data 

7.  Long-term  MTF  of  atmosphere 

8.  Short-term  MTF  of  atmosphere 

9 .  Wander 

10.  Long-term  MTF  of  atmosphere  calculated  from  FIT 

11.  Short-term  MTF  of  atmosphere  calculated  from  FIT 

In  addition,  the  following  sample  data  sets  are  included: 

1.  Optics  diffraction  pattern 

2.  LSF  of  optics  diffraction 

3.  MTF  of  optics 

4.  MTF  of  TV 

5.  Typical  input  data 


MODULE  TYPE: 

USAGE: 

FUNCTION: 

REFERENCED  BY: 
INTERFACES : 


APPENDIX  C,  PART  1.  MODULE  DESCRIPTIONS 

ABEL 

FORTRAN  7X  subroutine 
CALL  ABEL  (NPTS,  FN) 

Converts  a  line  spread  function  into  a  point  spread  func¬ 
tion 

PREDICT 

External  References:  ABS,  FLOAT,  SQRT 
Common  Data:  None 
Farnal  Arguments : 

(1)  NPTS  (integer) .  number  of  points  in  array  FN  to 
be  processed 

(2)  FN  (real  array  of  size  NPTS).  items  in  the  array 


are  functional  values;  on  input,  FN  is  'the  LSF; 
on  output,  it  is  the  PSF 


MODULE  TYPE: 


USAGE: 


FUNCTION: 


REFERENCED  BY: 


INTERFACES: 


FORTRAN  7X  REAL  function  subprogram 
Y  -  AIRY  (X) 

Calculates  the  value  of  the  Airy  function  of  X. 

AIRY(X)  •  2  *  J1(X)/X,  where  J1  is  the  Bessel  function  of 
order  one. 

OPTICS 

External  References:  ABS,  FLOAT,  MAXI 
Common  Data:  None 
Formal  Arguments : 


X  (real),  value  of  the  independent  variable  of  the 
Airy  function 


Inputs :  X 


Outputs :  AIRY 


MODULE  TYPE: 

USAGE: 

FUNCTION: 

REFERENCED  BY: 
INTERFACES : 


FORTRAN  7X  subroutine 

CALL  AMTF  (LTERM,  SQCN,  EX,  NPTS,  MTF) 
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Calculates  the  atmospheric  MTF  from  given  values  of  CQ  , 

range,  wavelength,  atmospheric  extinction,  and  for  short 
term  data,  the  optics  diameter. 

MTFLTV,  MTFLD,  PREDICT 

External  References:  COS,  EXP 

Common  Data:  /PARAMS/  WAV,  DOPT,  OBS,  SCALE,  RANGE,  EXT, 
A,  S SIGMA 

Formal  Arguments : 

(1)  LTERM  (integer) .  if  set  to  zero  on  input,  then 
the  long  term  atmospheric  MTF  is  calculated;  if 
set  to  one,  then  the  short  term  calculation 

(2)  SQCN  (real).  the  value  of  CQ2 

(3)  EX  (real),  the  atmospheric  extinction  factor 

(4)  NPTS  (Integer) .  the  number  of  elements  in  array 
FN 

(5)  MTF  (real  array  of  size  NPTS) .  the  calculated 
atmospheric  MTF 


Inputs :  LTERM,  SQCN,  EX,  NPTS,  WAV,  DOPT,  SCALE,  RANGE 


i 


MODULE  TYPE: 

USAGE: 

FUNCTION: 

REFERENCED  BY: 
INTERFACES : 


BITINV 

FORTRAN  7X  subroutine 
CALL  BITINV  (K,  L,  M) 

Finds  the  bit  inverse  of  an  integer,  i.e.,  its  base  2 
mirror  image,  of  a  given  bit  length. 

FFT 

External  References:  BTEST,  IBSET 
Common  Data:  None 
Formal  Arguments: 


(1) 

K  (integer), 
to  be  found 

the  integer  whose  bit 

inverse  is  i 

(2) 

L  (integer) . 

the  bit  inverse  of  K 

(3) 

M  (integer). 

the  number  of  bits  in 

the  integers  3 

K  and  L 
Inputs :  K,  M 
Outputs :  L 


DIGITIZE 


MODULE  TYPE; 

USAGE: 

FUNCTION: 

REFERENCED  BY: 
INTERFACES: 


NOTE: 


FORTRAN  7X  subroutine 
CALL  DIGITIZE 

Digitizes  one  frame  of  video  data  using  the  Quantex  DS-30 
Video  Processor.  Control  is  via  the  HPIB. 

DINPUT 

External  References:  EXEC 


Common  Data:  None 
Formal  Arguments :  None 
Inputs :  None 
Outputs :  None 

Currently  the  logical  unit  number  of  the  Quantex  DS-30 
(LUDS30)  is  set  to  24  and  that  of  the  HPIB  controller 
card  (LUHPIB)  is  set  to  20.  These  should  be  set  for  the 
particular  system  configuration  in  use. 


DIGITS 


MODULE  TYPE: 

USAGE: 

FUNCTION: 

REFERENCED  BY: 
INTERFACES: 


FORTRAN  7X  CHARACTER *16  function  subprogram 
CHAR  -  DIGITS  (N,  BASE) 

Converts  an  Integer  to  a  character  string  of  digits  which 
represents  the  integer  in  a  given  base. 

OUTPUT 

External  References:  IABS ,  MOD 
Common  Data:  None 
Formal  Arguments: 

(1)  N  (integer) .  the  Integer  to  be  converted 

(2)  BASE  (integer) .  the  number  base  in  which  the 
input  integer  is  to  be  represented 


D INPUT 


MODULE  TYPE: 

USAGE: 

FUNCTION: 

REFERENCED  BY: 
INTERFACES: 


FORTRAN  7X  subroutine 

CALL  DINPUT  (NPTS,  DLTA,  DSTA,  DWD) 

Obtains  digitized  video  data,  calculates  the  long  and 
short  term  average  data,  and  determines  the  wander. 

MTFLTV,  MTFLD 

External  References:  DIGITIZE,  VIDEODATA,  FLOAT 

Common  Data:  /WINDOW/  XW1,  XW2,  YW1,  YW2,  NFLDS; 

/PARAMS/  WAV,  DOPT,  OBS,  SCALE,  RANGE,  EXT, 

A,  SSIGMA 

Formal  Arguments: 

(1)  NPTS  (integer) »  number  of  points  in  arrays  DLTA, 
DSTA,  DLTI,  DSTI 

(2)  DLTA  (real  array  of. size  NPTS).  long  term  aver¬ 
age  of  the  LSF 

(3)  DSTA  (real  array  of  size  NPTS) ♦  short  term  aver¬ 
age  of  the  LSF 

(4)  DWD  (real  array  of  siae  NFLDS).  wander  of  the 
center  of  the  LSF 

Inputs :  NPTS,  NFLDS,  XW1,  XW2,  YW1,  YW2,  SCALE 
Outputs:  DLTA,  DSTA,  DWD 


DIVIDE 


MODULE  TYPE: 

USAGE; 

FUNCTION: 


REFERENCED  BY: 
INTERFACES: 


FORTRAN  7X  subroutine 

CALL  DIVIDE  (NPTS,  ARRAY1,  ARRAY2) 

Divides  each  element  of  first  array  (ARRAY1)  by  the  cor' 
responding  element  of  the  second  array  (ARRAY2).  If  a' 
divisor  is  zero,  then  the  quotient  is  set  to  zero. 

MTFLTV,  MTFLD 

External  References:  None 

Common  Data:  None 

Formal  Arguments: 

(1)  NPTS  (integer) .  number  of  elements  in  ARRAY1 
and  ARRAY2 

(2)  ARRAY!  (real  array  of  size  NPTS) .  dividend 

(3)  ARRAY2  (real  array  of  size  NPTS) .  divisor 
Inputs :  NPTS 

Outputs:  ARRAY1,  ARRAY2 


At. 


*r 


FFT 

FORTRAN  7X  subroutine 


CALL  FFT  (LINV,  NPTS,  FN) 

Calculates  the  Fourier  transform  (or  inverse  Fourier  trans¬ 
form)  of  a  real-valued  symmetric-  function  using  the  Fast 
Fourier  Transform  (FFT)  algorithm.  The  number  of  func¬ 
tional  values  must  be  an  integral  power  of  two. 

MTFLTV ,  MTFLD,  PREDICT 

External  References:  BITINV,  CEXP,  CMPLX,  CONJG,  FLOAT 
Common  Data;  None 
Formal  Arguments: 

(1)  LINV  (Integer) .  specifies  whether  forward 
(LINV*0)  or  inverse  (LINV»1)  Fourier  transform 
is  to  be  calculated 

(2)  NPTS  (integer) .  on  input,  NPTS  is  the  number  of 
elements  in  array  FN;  on  output,  it  is  the  largest 
non-zero  term  of  FN 

(3)  FN  (real  array  of  size  NPTS).  on  input,  FN  is 
an  array  of  functional  values;  on  output,  it  is 
the  forward  (or  Inverse)  Fourier  transform  of 
these  functional  values. 

Inputs:  LINV,  NPTS,  FN 

Outputs:  NPTS,  FN 

NOTE:  Complex  arithmetic  is  used  in  this  routine. 


MODULE  TYPE: 

USAGE: 

FUNCTION: 

REFERENCED  BY: 
INTERFACES: 
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MODULE  TYPE: 

USAGE: 

FUNCTION: 

REFERENCED  BY: 
INTERFACES : 


FORTRAN  7X  subroutine 

CALL  FIT  (CNSQ,  COEFF,  NPTS,  MFT) 

Performs  a  least  squares  fit  of  the  measured  MTF  data  to 
Fried's  model  of  the  atmospheric  MT F.  This  operation  de- 
2 

termines  the  C  parameter  for  the  atmosphere, 
n 

MTFLTV,  MTFLD 

External  References:  ALOG,  EXP,  FLOAT 

Comnon  Data:  /PARAMS/  WAV,  DOPT,  OBS,  SCALE,  RANGE,  EXT, 

A,  SSIGMA 

Formal  Arguments: 

2 

(1)  CNSQ  (real).  Cq  ,  the  atmospheric  turbulence 
struc ture  parameter 

(2)  COEFF  (real) ♦  the  coefficient  factor  in  the 
MTF  formula 

(3)  NPTS  (Integer) .  number  of  points  in  the  MTF  array 

(4)  MTF  (real  array  of  size  NPTS).  MTF  functional 
values 


Inputs :  NPTS,  MTF,  WAV,  DOPT,  SCALE,  RANGE 


GINPUT 


MODULE  TYPE: 

USAGE; 

FUNCTION: 

REFERENCED  BY: 
INTERFACES: 


FORTRAN  7X  subroutine 

CALL  GINPUT  (LPROG,  10) 

Obtains  general  parameter  inputs  for  programs  MTFLTV, 

MTFLD,  and  PREDICT.  The  parameters  may  be  input  interac¬ 
tively  or  via  data  file. 

MTFLTV,  MTFLD,  PREDICT 

External  References:  None 

Common  Data:  /PARAMS/  WAV,  DOPT,  OBS,  SCALE,  RANGE,  EXT, 

A,  S SIGMA 

/WINDOW/  XW1,  XW2,  YW1,  YW2,  NFLDS 
/HEADER/  HEADER1,  HEADER2,  DATE 
/PPARMS/  LSOURCE,  LOPTICS,  SQCN 

Formal  Arguments: 

(1)  LPROG  (integer),  parameter  specifying  which  pro¬ 
gram  is  being  run.  If  1,  then  MTFLTV  is  being 
run;  if  2,  then  MTFLD;  if  3,  then  PREDICT. 

(2)  10  (integer) .  parameter  specifying  whether  or 
not  all  data  is  to  be  output  to  files.  If  0, 
then  data  is  not  output;  if  1,  then  it  is. 

Inputs :  LPROG 

Outputs :  10,  WAV,  DOPT,  OBS,  SCALE,  RANGE,  EXT,  A,  SSIGMA, 

XW1,  XW2,  YW1,  YW2,  NFLDS,  HEADER1 ,  HEADER2,  DATE, 
LSOURCE,  LOPTICS,  SQCN 


MODULE  TYPE: 

FORTRAN  7X  subroutine 

USAGE: 

CALL  INT  (LNORM,  NPTS,  FN) 

FUNCTION: 

Integrates  the  point  spread  function  to  find  the  power 
inside  a  circle  as  a  function  of  radius.  Will  return  the 
power  function  normalized  to  total  power  if  LNORM  «  1. 

REFERENCED  BY: 

PREDICT 

INTERFACES : 

External  References:  FLOAT 

Common  Data:  /PARAMS/  WAV,  DOPT,  OBS,  SCALE,  RANGE,  EXT, 
A,  S SIGMA 


Formal  Arguments: 

(1)  LNORM  (integer 1 


if  0,  returns  unnormalized 
power  function; 

if  1,  returns  normalized  power 
function 


(2)  NPTS  (integer) .  number  of  points  in  array  FN 

(3)  FN  (real  array  of  size  NPTS) .  items  in  the  ar¬ 
ray  are  functional  values;  on  input  FN  is  the  PSF; 
on  output  it  is  the  power  integral  of  the  PSF  as  a 
function  of  radius 

Inputs ;  LNORM,  NPTS,  FN,  SCALE,  RANGE 


Outputs :  FN 


MODULE  TYPE: 


FORTRAN  7X  subroutine 


USAGE: 

FUNCTION: 

REFERENCED  BY: 
INTERFACES : 


CALL  LSF  (NPTS,  FN) 

Converts  a  radially  symmetric  point  spread  function  into  a 
line  spread  function 

MTFLTV,  MTFLD,  PREDICT 

External  References:  FLOAT,  IFIX,  SQRT 


Common  Data:  None 
Formal  Arguments: 

(1)  NPTS  (integer) .  number  of  points  in  the  array  FN 

(2)  FN  (real  array  of  size  NPTS) .  items  of  the  array 
are  functional  values;  on  input,  FN  is  the  PSF; 
on  output,  it  is  the  LSF 


MODULE  TYPE: 

USAGE: 

FUNCTION: 

REFERENCED  BY: 
INTERFACES : 


FORTRAN  7X  subroutine 
CALL  LSFS  (NPTS ,  FN) 

Shifts  a  line  spread  function  so  that  it  is  centered  at 
the  origin 

MTFLTV,  MTFLD 

External  References :  None 


Common  Data:  None 
Formal  Arguments: 

(1)  NPTS  (Integer) .  number  of  points  in  array  FN 

(2)  FN  (real  array  of  3ize  NPTS).  items  in  the  array 
are  functional  values;  on  input,  FN  is  the  non- 
centered  LSF;  on  output,  it  is  the  centered  LSF 

Inputs:  NPTS ,  FN 


Outputs:  FN 


I 

i 

I 

I 


MODULE  TYPE: 

USAGE: 

FUNCTION: 


REFERENCED  BY: 
INTERFACES : 


FORTRAN  7X  subroutine 

CALL  MUL  (NPTS,  ARRAY1,  ARRAY2) 

Performs  an  element  by  element  multiplication  of  ARRAY1  by 
ARRAY2  and  returns  the  result  in  ARRAY1 

MTFLTV,  MTFLD,  PREDICT 

External  References:  None 


Common  Data:  None 
Formal  Arguments: 

(1)  NPTS  (integer) .  number  of  elements  in  ARRAY1 
and  ARRAY2 

(2)  ARRAY1  (real  array  of  size  NPTS) .  on  input,  ARRAY1 
is  the  first  factor  array;  on  output,  it  is  the 
resultant  product  array 

(3)  ARRAY2  (real  array  of  size  NPTS) .  the  second 
factor  array 

Inputs:  NPTS,  ARRAY1,  ARRAY2 


Outputs:  ARRAY1 


OPTICS 

FORTRAN  7X  subroutine 
CALL  OPTICS  (NPTS,  FN) 

Calculates  the  point  spread  function  (diffraction  pattern) 
of  the  telescope  optics 

MTFLTV,  MTFLD,  PREDICT 

External  References:  AIRY,  FLOAT 

Common  Data;  /PARAMS/  WAV,  DOPT,  OBS,  SCALE,  RANGE,  EXT, 

A,  S SIGMA 

•Formal  Arguments: 

(1)  NPTS  (integer) .  number  of  points  in  array  FN 

(2)  FN  (real  array  of  size  NPTS) .  items  in  the  array 
are  functional  values;  the  computed  optics  PSF 

Inputs;  NPTS,  DOPT,  OBS,  SCALE,  WAV 

Outputs :  FN 

■ 

i 


MODULE  TYPE: 

USAGE: 

FUNCTION: 

REFERENCED  BY: 
INTERFACES: 


C-17 


OUTPUT 


MODULE  TYPE:  FORTRAN  7X  subroutine 

USAGE:  CALL  OUTPUT  (LHEAD,  NPTS,  ARRAY) 

FUNCTION :  Outputs  calculations  of  the  PSF,  LSF,  MTF,  etc.  to  data 

files. 


REFERENCED  BY: 


MTFLTV,  MTFLD 


INTERFACES-: 


External  References:  DIGITS 


Common  Data:  /HEADER/  HEADER1,  HEADER2,  DATE 
Formal  Arguments : 

(1)  LHEAD  (integer),  parameter  specifying  which 
title  is  to  be  printed  at  top  of  the  data  file 

(2)  NPTS  (integer) .  the  number  of  elements  in  array 
ARRAY 

(3)  ARRAY  (real  array  of  size  NPTS) .  array  of  values 
which  are  to  be  output  to  the  data  file 


Inputs :  LHEAD,  NPTS,  ARRAY,  HEADER1,  HEADER2,  DATE 


Outputs:  None 


MODULE  TYPE: 


FORTRAN  7X  subroutine 


USAGE: 


FUNCTION: 


REFERENCED  BY: 


INTERFACES : 


CALL  POUTPUT  (LHEAD,  NPTS,  ARRAY) 


Outputs  calculations  of  the  PSF,  LSF,  MTF,  etc.  to  data 
files . 


PREDICT 

External  References:  DIGITS 

Common  Data:  /HEADER/  HEADER1,  HEADER2,  DATE 

Formal  Arguments: 


(1)  LHEAD  (integer) .  parameter  specifying  which 
title  is  to  be  printed  at  top  of  the  data  file 


(2)  NPTS  (integer) .  the  number  of  elements  in  array 
ARRAY 


(3)  ARRAY  (real  array  of  size  NPTS).  array  of  values 
which  are  to  be  output  to  the  data  file 


Inputs :  LHEAD,  NPTS,  ARRAY,  HEADER1,  HEADER2 
Outputs :  None 
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cl 


i 
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SOURCE 


MODULE  TYPE: 

USAGE: 

FUNCTION: 

REFERENCED  BY: 
INTERFACES: 


FORTRAN  7X  subroutine 
CALL  SOURCE  (NPTS,  FN) 

Calculates  the  values  of  a  Gaussian  source  distribution 
function 

MTFLD* 

External  References:  EXP,  FLOAT 

Common  Data:  /PARAMS/  WAV,  DOPT,  OBS,  SCALE,  RANGE,  EXT, 

A,  SSIGMA 

Formal  Arguments: 

(1)  NPTS  (integer),  number  of  points  in  array  FN 

(2)  FN  (real  array  of  size  NPTS)-  items  in  the  array 
are  functional  values;  the  computed  Gaussian 
source  distribution  function 


Inputs:  NPTS,  A,  SCALE,  SSIGMA 


MODULE  TYPE: 


FORTRAN  7X  subroutine 


CALL  TVMTF  (NPTS,  FN) 

Obtains  the  modulation  transfer  function  of  the  TV  system 
MTFLTV,  MTFLD 
External  References:  None 
Common  Data:  None 
Formal  Arguments: 

(1) _  NFTS  (integer),  number  of  elements  in  array  FN 

(2)  FN  (real  array  of  size  NPTS).  items  in  the  array 
are  functional  values;  the  TV  system  MTF 

Inputs :  NPTS 

Outputs:  FN 

NOTE:  Currently,  this  module  calculates  a  simple  step  function 

for  the  TV  system  MTF.  It  would  have  to  be  modified  to 
allow  for  a  more  complicated  MTF. 


USAGE: 
FUNCTION: 
REFERENCED  BY: 
INTERFACES: 
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MODULE  TYPE: 


USAGE: 


FUNCTION: 


REFERENCED  BY: 


INTERFACES: 


NOTE: 


VIDEODATA 

FORTRAN  7X  subroutine 
CALL  VIDEODATA  (ROW,  ROWSUM) 

Finds  the  sum  of  pixel  values  along  a  given  row.  The  row 
must  lie  between  the  top  and  bottom  window  limits,  YW2  and 
YW1,  and  the  length  of  the  row  is  determined  by  the  left 
and  right  window  limits,  XW1  and  XW2 . 

DINPUT 

External  References:  ABRT,  EXEC,  LAND,  ISHFT 
Common  Data:  /WINDOW/  XW1,  XW2,  YW1,  YW2,  NFLDS 
Formal  Arguments: 

(1)  ROW  (Integer) .  the  digitized  video  field  row 
which  is  to  be  summed 

(2)  ROWSUM  (real) .  the  sum  of  pixel  values  along 
the  row  of  pixels  given  by  ROW 

Inputs:  ROW,  XW1,  XW2,  YW1,  YW2 

Outputs :  ROWSUM 

Currently  the  logical  unit  number  of  the  Quantex  DS-30 
(LUDS30)  is  set  to  24  and  that  of  the  HPIB  controller 
card  (LUHPIB)  is  set  to  20.  These  should  be  set  for  the 
particular  system  configuration- in  use. 


m 
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WANDER 


MODULE  TYPE: 

USAGE: 

FUNCTION: 

REFERENCED  BY: 
INTERFACES: 


FORTRAN  7X  subroutine 

CALL  WANDER  (CNSQ,  NPTS,  FN) 

2 

Calculates  the  value  o£  CQ  from  the  wander  of  the  laser 
spot  Image 

MTFLTV,  MTFLD 

External  References:  FLOAT 

Common  Data:  /PARAMS/  WAV,  DOPT,  OBS,  SCALE,  RANGE,  EXT, 

A,  SSIGMA 

Formal  Arguments.: 

2 

(1)  CNSQ  (real),  value  of  C  ,  calculated  from  wander 

data  n 

(2)  NPTS  (integer) .  number  of  elements  in  array  FN 

(4)  FN  (real  array  of  size  NPTS) .  items  in  the  array 
are  functional  values;  locations  of  the  laser 
spot- image  center  as  a  function  of  time 

Inputs:  NPTS,  FN,  DOPT,  RANGE,  SCALE 


Outputs :  CNSQ 


I 
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RUN  3 

LONG  TERM  RAW  DATA 


..063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

.938 

..063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

..063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

.938 

1.938 

-.063 

-.063 

..063 

-.063 

-.063 

-.063 

..063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

1.938 

..063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

..063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

..063 

-.063 

-.063 

..063 

-.063 

-.063 

-.063 

-.063 

..063 

..063 

-.063 

-.063 

..063 

-.063 

-.063 

-.063 

..063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

—.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

989  .*9 U 

-.063 

-.063 

.938 

99.938 

89.938 

396.938 

397.938 

996.938 

2021.938 

2007.938 

31*9.938 

3117.938 

*289.938 

*31*. 938 

i??66i5:??? 

6831.938 

6967.938 

831*. 938 

8555.938 

9715.938 

10223.937 

1*590.938 

1*939.937 

16851.937 

17285.938 

2022*. 937 

20260.938 

22257.937 

2222 1.937  25819.937 

2601*. 937 

30*08.937 

30357.937 

35793.937 

35577.937 

38869.937 

38577.937 

**15*. 937 

*3357.937 

*700*. 937 

*6617.937 

52*80.938 

522*9.937 

5*127.937  53851.937 

57050.938  56870.938 

55355.938 

55250.938 

55083.937 

5*973-937 

50722.937 

507**. 937 

*5910.938 

*5677.937 

38287.937 

38180.938 

31911.937 

313*6.937 

23250.938  23058.937 

17328.937 

17162.937 

12**2.937 

12297.937 

87*9.938 

8528.938 

5391.938 

5313.938 

2760.938 

2701.938 

1072.938 

1075.938 

376.938 

3*7.938 

93.938 

92.938 

-.063 

12.938 

7.938 

.938 

.938 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

..063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

..063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

—  063 

..063 

..063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

..063 

-.063 

-.063 

—  063 

.  -.063 

..063 

-.063 

..063 

-.063 

-.063 

-.063 

—  063 
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RUN  3 


SHORT  TERH  RAW  DATA 


-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

262.938 

11449.937 

59802.937 

89402.937 

TO;flr 

7!:!H 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

.938 


-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

505.938 

16540.938 

68249.937 

85151.937 

40605.938 

4239.938 

18.938 

-.063 

-.063 

—.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 


—.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

695.938 

20164.937 

72840.938 

80405.938 

35206.937 

2952.938 
3-938 
-.063 
,  -063 
—.063 
-.063 
-.063 
-.063 
—  063 
-.063 
-.063 
-.063 
-.063 
-.063 


-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

.1.938 

-.063 

-.063 

-.063 

-.063 

-.063 

6.938 

1535.938 

26804.937 

81012.937 

75522.937 

27411.937 

1679.938 

1.938 
-.063 
-.063 
-.063 
-.063 
-.063 
-.063 
-.063 
-.063 
-.063 
-.063 
-.063 
-.063 


-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

1.938 

-.063 

10.938 

2101.938 

32171.937 

84893.937 

71458.937 

23579.937 

1169.938 
-.063 
-.063 
-.063 
-.063 
-.063 
-.063 
-.063 
-.063 
-.063 
-.063 
-.063 
-.  063 
-.063 


-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

39.938 

3936.938 

39221.937 

87399.937 

65740.938 

16666.937 

453.938 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 


.938 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

44.938 

5345.938 

45499.937 

90120.938 

61156.937 

13501.937 

365.938 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 


1.938 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

153-938 

8966.938 

54432.937 

90492.937 

53901.937 

8977.938 

106.938 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

-.063 

—.063 

-.063 
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RUN  3 
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FFT  OF  SHORT  TERM 


[ 

\ 

L 


' 

i 

I 


107520.70 

107411.22 

107083.37 

106538.97 

105780.97 

104813.47 

103641.73 

100711.59 

98968.891 

97053.109 

76836.797 

94974.359 

92743.453 

90371.844 

87871.531 

82534.844 

79724.375 

73885.625 

70884.359 

67846.344 

64784.711 

58641.414  55584.094 

52551.812 

49555.508 

46605.586 

43711.781 

40883.078 

35453-070 

32865.641 

30371.168 

27974.574 

25680.047 

23490.980 

21410.020 

17578.723 

15829.682 

14191.365 

12662.740 

11242.244 

9927.834 

8716.969 

6593.698 

5674.041 

4843.556 

4097.937 

3432.622 

2842.980 

2324.379 

1481.447 

1147.648 

865.982 

631.732 

440.422 

287.711 

169.612 

21.257 

-16.002 

-33.231 

-33-596 

-20.051 

4.760 

38.512 

124.786 

173.864 

224.868 

276.438 

327.372 

376.678 

423.570 

508.070 

545.021 

578.1-06 

607.191 

632.140 

652.917 

669.591 

691.490 

697.330 

700.191 

700.310 

697.906 

693.175 

686.324 

667.524 

656.296 

644.322 

631*871 

619.141 

606.250 

593.309 

567.793 

555.567 

543.955 

533.059 

522.907 

513.455 

504.634 

488.725 

481.710 

475.366 

469.770 

464.859 

460.521 

456.569 

449.303 

422.228 

445.821 

442.420 

439.146 

436.070 

432.935 

429.668 

418.040 

413.522 

408.707 

403.637 

398.286 

392.588 

379.750 

372.503 

364.769 

356.659 

348.284 

339.700 

330.943 

312.724 

303.183 

293.364 

283.384 

273-361 

263.383 

253.521 

234.076 

224.434 

155.788 

214.873 

205.452 

196.288 

187.474 

179.052 

163.291 

148.460 

141.307 

134.389 

127.795 

121.588 

110. 316 

105.118 

100.101 

95.224 

90.494 

85.952 

81.642 

73.758 

70.063 

66.414 

62.779 

59.152 

55.584 

52.157 

45.881 

42.977 

40.124 

37.234 

34.287 

31.302 

28.379 

22.866 

20.398 

18.083 

15.867 

13.715 

11.651 

9.735 

6.648 

5.551 

4.726 

4.104 

3.641 

3- 308 

3. 161 

3-774 

4.659 

5.927 

7.523 

9.353 

11.352 

13.507 

18.435 

21.315 

24.504 

27.951 

31.599 

35.369 

39.236 

47.365 

51.724 

56.308 

61.058 

65.901 

70.718 

75.445 

84.582 

89.095 

93-632 

98.187 

102.707 

107.110 

111.280 

118.834 

122.239 

125.423 

128.436 

131.238 

133.756 

135.925 

139.172 

140.14.5 

141.033 

141.527 

141.967 

141.993 

141.765 

140.066 

138.633 

136.893 

134.936 

132.827 

130.594 

128.197 

122.826 

119.882 

116.828 

113.773 

110.792 

107.900 

105.077 

99.606 

96.940 

94.386 

92.051 

90.015 

88.340 

87.002 

85.130 

85.683 

84.468 

83-970 

83.652 

83.581 

83.769 

84.219 

86.563 

87.476 

88.407 

89.400 

90.449 

91.538 

93.560 

92.178 

94.299 

94.760 

94.923 

94.814 

94.474 

93-928 

90.852 

89.141 

87.021 

84.513 

81.660 

78.519 

71.515 

67.662 

63.565 

59.235 

54.737 

50.164 

45.613 

36.779 

32.574 

28.480 

24.487 

20.621 

16.936 

13.511 

7.746 

5.483 

3.617 

2.111 

.933 

.080 

-.441 

-.380 

.182 

1.029 

2.099 

3.330 

4.682 

6.149 

9.471 

11.313 

13.225 

15.097 

16.848 

18.401 

19.703 

21.578 

22.181 

22.552 

22.678 

22.512 

22.026 

21.207 

18.676 

17.049 

15.233 

13.224 

11.008 

8.574 

5.926 

.195 

-2.732 

-5.615 

-8.413 

-11.118 

-13.764 

-16.376 

-21.550 

-24.075 

-26.518 

-28.868 

-31.116 

-33.344 

-35.585 

-40.220 

-42.612 

-45.064 

-47.577 

-50. 196 

-52.990 

-56.032 

-63.123 

-67.227 

-71.679 

-76.465 

-81.568 

-87.004 

-92.800 

-105.588 

•112.601 

-119.990 

-127.708 

-135.705 

-143.951 

-152.401 

-169.847 

-178.719 

-187.590 

-196.342 

-204.888 

-213.137 

-221.043 

-235.713 

-242.419 

-248.630 

-254.281 

-259.285 

-263.571 

-267.092 

-271.768 

-272.912 

-273.258 

-272.794 

-271.510 

-269.405 

-266.492 

-258.415 

-253.352 

-247.661 

-241.359 

-234.460 

-226.978 

-218.929 

-201.346 

-191.959 

-182.279 

-172.371 

-162.290 

-152.063 

-141.724 

C-31 


102271.97 

85254.922 

61712.242 

38127.859 

19438.902 

7606.731 

1872.094 

82.094 

79.110 

467.492 

682.338 

677.714 

580.431 

496.387 

452.865 

426.089 

386.435 

321.962 

243.751 

171.046 

115.775 

77.590 

48.926 

25.573 

8.045 

3.290 

15.845 

43.240 

80.056 

115.198 

137.704 

141.215 

125.611 

102.317 

85.953 

84.892 

92.603 

93.179 

75.131 

41.146 

10.423 

-.594 

7.738 

20.772 

20.074 

3.110 

-18.986 

-37.877 

-59.394 

-98.984 

-161.050 

-228.572 

-269.827 

-262.812 

-210.362 

-131.315 


RUN  3 


H 


LONG  TERN  MTF  OE  THE  ATMOSPHERE 


1787646.2 

1681069.7 

1393045.5 
1005998.1 
619852.25 
313950.31 
119888.44 
27470. 137 
402.456 
2417.552 
10098.687 
14132.754 
14046.871  - 
12305.914 
10905.391 
10315.865 
9969.502 
9217.563 
7800.448 
5997.551 
4295.559 
2977.896 
2042.164 
1302.396 
665.160 
198.042 
115.065 
574.897 
1509.727 
2753.863 
3949.301 
4716.795 


1785938.5 
1653680.0 

1348034.5 

955719.50 
576025.12 
282955.44 
103417.89 
21334.129 

-303.787 

3378.699 

10869.766 

14303.482 

13861.637 

12085.207 

10787.613 

10272.070 

9904.924 

9072.695 

7588.061 

5769.643 

4111.550 

2846.650 

1945.943 

1223.705 

595.144 

165.853 

142.460 

666.563 

1653.105 

2908.368 

4072.652 

4761.070 


1780822.0 

1623515.7 

1301650.5 
905671.88 
533615.63 
254291.84 
88496.031 
16138.961 
-632.636 
4383.441 
11568.992 
14414.184 
13659.049 
11876.045 
10683.832 
10229.637 
9831.883 
8914.779 
7367.139 
5542.176 
3930.892 
2719.449 
1850.387 
1146.004 
529.190 
.  141.626 
181.760 
768.380 
1804.443 
3064.426 
4189.171 
4802.537 


1772323.0 

1590721.7 

1254097.7 
856036.62 
492720. 06 
227455.41 
75055.609 
11803.486 
-641.394 
5405.696 

12192.943 

14469.088 

13444.641 

11680.775 

10595.934 

10189.648 

9751.088 

8746.398 

7140.500 

5316.739 
3753.640 
2595.197 
1754.595 

1066.740 
465.733 
123.358 
231.355 
878.899 

1961.870 

3221.823 

4300.440 

4830.652 


1760482.2 

1555452.2 

12 05582.7 

806989.50 
453422.19 
202 432.37 
63023.898 

8250.271 

-383.887 

6422.125 

12738.100 

14472.088 

13222.510 

11500.193 

10522.670 

10153.771 

9663.443 

8570.201 

6911.128 

5096.350 

3581.413 

2474.122 

1658.403 

985.315 

403.782 

109.738 

288.438 

996.327 

2123.139 

3378.795 

4405.067 

4856.926 


1745356.2 

1517872.0 

1156312.0 

758696.62 

415794.19 

179200.09 

52326.117 

5403.681 

91.405 

7413.246 

13202.928 

14426.656 

12995.064 

11333.475 

10461.641 

10116.049 

9568.273 

8387.496 

6681.239 

4883.546 

3416.695 

2357.396 

1563.214 

902.280 

344.035 

100.005 

351.058 

1118.234 

2284.365 

3532.684 

4500.584 

4869.002 


1727018.2 

1478155.7 

1106490.0 

711314.37 

379895.25 

157726.50 

42886.648 

3194.004 

741.663 

8363.357 

13587.967 

14336.590 

12764.639 

11179.262 

10408.777 

10074.812 

9463.922 

8199.098 

6452.605 

4679.479 

3261.203 

2246.267 

1471.398 

820.533 

288.305 

95.835 

418.864 

1243.859 

2443.511 

3679.607 

4584.672 

4872.316 


1705554.7 

1436483-7 

1056319.0 

664991.00 

345769.50 

137972.31 

34627.336 

1550.062 

1528.023 

9261.135 

13895.924 

14208.873 
12533-646 
11036.441 
10360.971 
10025.732 
9347. 660 
8003.727 
6224.679 
4484.902 
3115.270 
2141.520 
1384.535 
741.631 
238.959 
100.022 
492.748 
1374.530 
2599.673 
3818.815 
4655.909 
0.000 


RUN  3 


SHORT  TERM  MTF  OF  THE  ATMOSPHERE 


1787646.2  1782230.5  1766084.5  1739515.7  1703024.7  1657288.7  1603146.0  1541566.2 

1473625.7  1400474.7  1323307.2  1243328.2  1161721.0  1079618.5  998078.87  918059-38 

840402.12  765817.87  694879.87  628022.75  565543.00  507605.81  454256.31  405433.19 

360985.50  320688.31  284264.31  251400.25  221765. 8T  195025.34  170870.94  148997.12 

129137.59  111070.73  94606.984  79596-.812  65927.328  53517.461  42311.930  32274.168 

23380.426  15610.396  8942.182  3346.300  .1219.216  -4811.179  >7502.373  >9380.240 

>10546.38  >11113.26  -11201.05  -10931.67  -10426.68  -9801.781  -9160.836  -8593.234 

-8169. 332  -7938.653  -7929.545  -8149.484  -8586.289  -9212.072  -9983*271  -10851.11 

-11757.62  -12632.94  -13415.64  -14044.73  -14467.55  -14641.76  -14538.06  -14142.59 

-13453.94  -12484.39  -11256.45  -9801.342  -8157.185  -6367.749  -4480.472  -2544.378 

-608. 494  1281.388  3084.422  4768.081  6307.752  7687.098  8897.814  9937.693 

10810.080  11521.584  12082.672  12505.729  12805.203  12996.031  13092.340  13108.154 

13050.945  12930.121  12748.908  12508*.  396  12208.049  11845.840  11419.312  10926.178 
10364.719  9733.762  9034.266  8268.719  7442.167  6561.742  5637.688  4682.184 

3709.214  2734.493  1773.688  843.280  -41.107  -864.968  -1615.026  -2279.377 

-2848.235  -3314.277  -3672.706  -3920.965  -4057.905  -4092.003  -4029.241  -3880.464 

-3652.440  -3358.165  -3009.818  -2619.146  -2198.236  -1759.041  -1314.564  -877.723 

-460.740  -73.849  275.512  582.965  847.738  1070.736  1254.070  1400.482 

1512-621  1593.326  1646.618  1679.315  1699.565  1716.996  1741.529  1781.269- 

1842.491  1927.659  2037.674  2171.445  2327.517  2504.442  2700.047  2912.192 

3134.235  3360.960  3583.376  3791-586  3975.979  4128.095  4241.121  4309-537 

4329.185  4296.613  4209.160  4064.646  3862.159  3602.481  3287.504  2921.451 

2509.146  2055.966  1567.871  1050.672  510.791  -45-375  -610.600  -1177.560 

-1738.515  -2286.165  -2813-756  -3314.929  -3783.891  -4215.488  -4604.079  -4946.930 

-5240.757  -5479-330  -5661.122  -5784.210  -5848.120  -5853-926  -5804.354  -5703.984 

-5557.754  -5372.167  -5153-277  -4906. 913  -4638.929  -4355.080  -4062.106  -3767.456 

-3479-262  -3206.071  -2955.256  -2731.747  -2538.409  -2375.226  -2241.366  -2135.940 

-2058.653  -2010.976  -1994.633  -2011.837  -2062.924  -2146.335  -2260.073  -2400.602 

-2568.284  -2762.475  -2987.378  -3247.971  -3549.952  -3897.104  -4292.535  -4736.799 

-5229.526  -5770.157  -6358.749  -6995.352  -7679.693  -8409.422  -9180.654  -9986.021 

-10816.28  -11660.98  -12509.44  -13349.93  -14171.96  -14963.88  -15714.55  -16411.86 

-17041.85  -17602.25  -18072.99  -18458.83  -18744.95  -18938.84  -19031.71  0.000 


C-33 


WANDER 


5.700 

-22.800 

-*2.750 

-25. 

650 

2.850 

11. *00 

19.950 

8. 

,550 

-5.700 

5.700 

-5.700 

19. 

950 

31.350 

*2.750 

*2.750 

28. 

,500 

57.000 

71.250 

51.300 

5*. 

150 

*2.750 

51.300 

31.350 

25. 

,650 

*5.800 

31.350 

*2.750 

37. 

050 

17.100 

31-350 

51.300 

51. 

300 

51.300 

-34.200 

-25.650 

-42. 75C 

14.250 

11.400 

2.850 

-11.400 

0.000 

31.350 

42.750 

17.100 

39.900 

54.150 

54.150 

51.300 

*8. *50 

31.350 

31.350 

39.900 

45.600 

51.300 

57.000 

48.450 

51.300 

5*. 150 

31.350 

31.350 

LONG  TERM  MTF  OF  THE  ATMOSPHERE  CALCULATED  FROM  FIT. 

1891009.7  1890673*7  1889943.5  1888914.5  1887626.2  1886104.0  1884365.2  1882423.2 
1880289.2  1877971.7  1875479.5  1872818.5  1869995.5  1867015.5  1863884.0  1860605.0 


1857183.2  1853622.5  1849926.7  1846099.0 

1825096.2  1820542.7  1815879.2  1811105.7 

1785691.2  1780315.2  1774844.2  1769284.5 

1740182.7  1734115.5  1727973.0  1721751.5 

1689569.7  1682930.7  1676225.2  1669455.7 

1634729.2  1627618.0  1620454.7  1613240.2 

1576454.5  1568963.0  1561432.7  1553861.2 
1515455.0  1507674.5  1499858.0  1492016.2 

1452398.5  1444406.0  1436388.7  1428343.5 

1387886.2  1379750.7  1371600.5  1363437.5 
1322484.0  1314268.0  1306051.5  1297829.5 

1256684.7  1248457.0  1240227.5  1232008.5 

1190959.7  1182761.5  1174583.5  1166407.5 

1125710.7  1117606.5  1109516.2  1101438.2 
1061314.0  1053342.2  1045387.7  1037455.0 

998095.37  990293.1 2  982512.37  974754.37 

936343.12  928735.75  921162.87  913611.87 

876293.62  868924.25  861576.37  854261.37 
818158.00  811037.87  803942.25  796886.5 0 
762096.88  755252.37  748426.87  741643.75 

708261.87  701696.12  695161.75  688663.25 

656743.25  650470.00  644239-75  638041.75 

607620.62  601654.62  595723.50  589828.1.3 

560945.87  555285.12  549661.87  544080.12 

516736.62  511388.31  506065-31  500792.94 

475000.37  469946.62  464947.37  459979.06 

435710.87  430970.81  426266.37  421603.69 
398847.44  394397.44  389996.12  385633.87 

364342.37  360201.00  356079.62  352009.87 

332150.50  328291.94  324459-94  320666.50 
302201.31  298605.44  295053-50  291525.06 

274400.50  271072.37  267778.56  264510.81 
248667.28  245595.34  242543.69  239533-09 
224909.78  222070.94  219269-50  216490.03 
203029.56  200418.72  197843.59  195285.31 
182925.34  180531.09  178171.41  175825.41 
164501.84  162314.87  160145.66  158006.78 
147658.41  145656.25  143677.06  141724.53 

132288.50  130467.34  128664.27  126884.50 
118301.39  116644.02  115008.75  113393.50 
105599-34  104097.55  102612.91  101145.81 
94089.594  92728.625  91388.469  90064.375 

83684.406  82456.281  81243.469  80051.141 

74296.406  73189-219  72098.656  71020.828 
65843.750  64849.727  63867.617  62899.516 
58250.695  57358.977  56476.914  55609.109 
51442.547  50644.602  49856.047  49078.242 
45352.844  44638.516  43934.328  43239.312 
39913.656  39277.430  38648.859  38029.437 
35065.687  34499-312  33939.727  33389.891 
30756.395  30251.926  29755.898  29267-418 
26929.781  26482.945  26042.484  25610.012 
23538.930  23144.129  22754.520  22372.191 
20542.344  20193.496  19849-301  19509.816 


1842142.5  1838061.5  1833858.2  1829535.7 

1806226.2  1801243.5  1796158.0  1790973.5 

1763634.2  1757897.2  1752074.7  1746170.5 

1715456.5  1709090.5  1702654.0  1696144.2 
166262 9.0  1655739.7  1648793.5  1641790.0 

1605973.5  1598664.2  1591303.2  1583901.2 
1546252.0  1538608.7  1530923.5  1523206.2 

1 484147.7  1476247.2  1468323-7  1460372.7 

1420288.7  1412216.7  1404120.0  1396012.2 
1355268.0  1347079-7  1338888.2  1330692.2 

1289602.2  1281374.2  1273145.2  1264918.2 

1223780.5  1215568.0  1207356.0  1199152.5 

1158243.7  1150094.2  1141951.7  1133824.5 

1093384.2  1085340.0  1077312.7  1069301.2 

1029538.7  1021647.2  1013772.0  1005920.9 
967027.00  959314.87  951628.37  943972.75 

906096.75  898598.12  891138.62  883697.50 
846973.00  839726.75  832501.62  825316.50 

789860.25  782864.75  775911.25  768985.87 

734895.50  728183-62  721508.00  714866.12 

682208.37  675786.62  669402.12  663055.25 

631883.62  625762.62  619674.87  613628.87 

583972.12  578163.25  572383-50  566648.12 

538533.25  533024.75  527555.12  522124.44 
495557.06  490361.81  485199-56  480084.50 

455048.87  450162.94  445298.69  440489.19 
416979.56  412390.81  407834.25  403319.81 
381307.31  377009.87  372754.81  368532.19 

347966.37  343958.81  339993-25  336051.06 
316899.69  313174.25  309481.25  305826.50 
288037.19  284575.87  281155.00  277757.69 
261280.09  258080.72  254915-69  251773-03 
236545.34  233597.62  230670.50  227776.75 

213740.12  211022.44  208329.00  205664.81 

192757.75  190261.00  187790.09  185344.97 
173509.00  171217.62  168955.16  166715.  1 2 
155885.47  153796.12  151722.16  149675.72 
139792.72  137881.31  135994.25  134129.44 
125127.94  123391.08  121674.75  119975.94 
111795.06  110216.47  108658.14  107122.03 
99703. 812  98272.578  96860.000  95467.234 
88754.984  87463.094  86186.969  84926.531 
78868.438  77704.484  76551.312  75416.656 
69955.797  68908.219  67873.937  66853.547 
61944.484  61003.086  60073.336  59154.664 
54753.156  53908.703  53074.914  52251.578 
48312.320  47558.062  46812.523  46077.492 
42555.102  41880.719  41216.000  40560.086 
37419.922  36818.562  36225.289  35641.586 
32847.734  32312.836  31786.176  31267.324 
28785.000  28311.184  27843.496  27383.687 
25183.141  24762.289  24348.781  23941.172 
21993.859  21622.652  21256.348  20897.969 
19176.762  18849.121  18527.766  18208.086 


RUN  3 


SHORT  TERM  MTF  OF  THE  ATMOSPHERE  CALCULATED  FROM  FIT. 


2020644.0  2020480.2  2020131.2  2019645.0 

2015672.5  2014621.0  2013495.0  2012297.5 

2005329.5  2003755.0  2002124.2  2000438.5 
1991238.0  1989251.2  1987219.0  1985141.0 

1974104.7  1971774.5  1969403.7  1966995.7 

1954395.5  1951768.0  1949108.0  1946412.7 

1932455.5  1929571.5  1926657.0  1923712.0 

1908566.5  1905456.2  1902319.2  1899156.5 

1882967.2  1879656.7  1876324.0  1872967.7 

1855858.7  1852374.5  1848867.5  1845342.5 

1827426.2  1823789.2  1820132.7  1816455.0 

1797831.2  1794059.0  1790270.2  1786466.0 

1767227.5  1763335.7  1759433*0  17555T6.5 

1735743.5  1731753.2  1727750.0  1723739.0 

1703511.5  1699431.0  1695347.0  1691250.0 

1670640.2  1666491.5  1662334.5  1658168.7 

1637240.2  1633033.7  1628819.7  1624600.2 

1603410.5  1599158.0  1594899.0  1590634.7 

1569246.5  1564954.5  1560663.0  1556364.7 
1534830.0  1530518.5  1526198.5  1521877.7 

1500244.2  1495914.7  1491579.0  1487246.5 
1465557.0  1461224.5  1456880.7  <1452542.0 

1430848.7  1426512.0  1422173.0  1417834.7 

1396171.5  1391842.0  1387518.0  1383191.7 
1361590.0  1357279.5  1352968.8  1348658.5 

1327182.5  1322872.0  1318584.2  1314302.2 

1292932.7  1288674.0  1284407.5  1280155.7 

1258952.7  1254715.7  1250498.7  1246279.5 

1225257.7  1221068.7  1216883.0  1212705.7 

1191900.2  1187746.0  1183609.0  1179480.5 

1158897.7  1154808.2  1150708.0  1146631.5 

1126298.5  1122260.7  1118220.2  1114191.2 

1094135.7  1090140.2  1086164.7  1082184.2 

1062421.2  1058490.0  1054569.0  1050648.0 

1031187.2  1027322.1  1023451.4  1019603.7 
1000455.0  996645.12  992855.50  989063.62 

970243.87  966499.50  962777.12  959046.37 

940564.12  93689.0.63  933240.00  929577.37 

911439.62  907844.00  904245.00  900665.75 

882883.37  879348.75  875823.12  872314.25 

854891.62  851435.62  847981.25  844540.87 
827486.00  824098.37  820725.50  817362.50 

800669.87  797359*62  794056.62  790760.12 

774448.12  771208.25  767987.37  764774.12 

748828.37  745666.25  742512.00  739381.62 
723809-00  720720.62  717648.12  714580.12 

699389.12  696381.25  693378.13  690387.87 

675574.50  672642.25  669711.12  666797.12 

652358.25  649504.50  646651.75  643808.62 

629748.37  626963.00  624187.37  621417.50 

607726.37  605019.50  602314.75  599619.87 

586291.87  583657.37  581024.87  578410.62 

565457.12  562888.62  560334.38  557790.12 

545194.12  542700.62  540213.87  537744.87 
525501.00  523083.37  520669.44  518272.56 
506387.94  504040.25  501694.25  499352.75 


2019042.7  2018336.5  2017535.2  2016645.5 

2011031.5  2009699.5  2008303.7  2006846.5 

1998699.7  1996909.0  1995067.7  1993177.0 
1983019.0  1980853-7  1978645.5  1976395.7 

1964548.7  1962064.7  1959544.0  1956987.7 

1943685.2  1940925.7  1938134.2  1935309.5 
1920740.0  1917738.2  1914709.0  1911652.0 

1895966.7  1892754.5  1889515.2  1886253-5 

1869589.7  1866190.5  1862767.2  1859323-5 
1841799.0  1838233.5  1834650.2  1831047.5 

1812764.7  1809058.0  1805331.0  1801589.7 

1782648.5  1778812.2  1774964.2  1771104.2 

1751586.2  1747644.0  1743689.5  1739724.2 

1719710.7  1715677.2  1711630.7  1707575.0 
1687144.5_1683032.2  1678909.0  1674778.5 
1653999-7  1649819.7  1645633-0  1641438.2 
1620373.0  1616142.0  1611902.2  1607658.2 

1586369.2  1582093-2  1577813-5  1573532.2 

1552067.5  1547760.5  1543455.7  1539142.0 
1517552.0  1513231.0  1508901.7  1504575.5 

1482909.7  1478570.5  1474235.5  1469896.0 

1448202.7  1443863*7  1439525.7  1435186.5 
1413502.0  1409167.7  1404834.7  1400503.5 

1378869.2  1374547.5  1370224.7  1365906.5 

1344351.2  1340054.0  1335752.2  1331458.0 

1310020.2  1305741.7  1301467.0  1297196.2 

1275906. 3  1271662.5  1267418.5  1263187.2 

1242064.5  1237861.0  1233646.0  1229452.7 
1208535-7  1204368.7  1200202.2  1196046.0 

1175356.5  1171231.0  1167117.7  1163006.5 
1142551.0  1138477.2  1134417.2  1130350.5 

1110159.5  1106143.2  1102132.0  1098133-0 

1078220.5  1074255.7  1070308.0  1066356.0 
1046741.0  1042842.0  1038953.9  1035062.7 
1015752.9  1011924.7  1008092.4  1004271.1 

985281.12  981513.00  977746.12  973989.87 

955329.62  951627.87  947932.50  944245.00 

925930.25  922291.25  918667.12  915047.50 

897083.37  893525.12  889960.12  886411.37 

868811.12  865312.87  861828.13  858352.87 

841114.62  837694.25  834283-87  830875-75 

814002.50  810653.25  807314.50  803994.62 

787491.12  784213.37  780947.25  777696.00 

761564.62  758367.00  755177.00  751994.87 

736243.50  733125.25  730003.62  726902.37 

711516.50  708473-62  705439.50  702414.12 

687406.25  684438.00  681474.37  678516.00 

663892.25  660996.88  658106.50  655221.88 

640979.12  638162.87  6353«8.50  632543-62 

618661.12  615914.62  613177.87  610447.25 

596938.87  594264.00  591595.37  588940.63 

575802.62  573200.50  570609.00  568028.37 
555248.75  552725.37  550204.6 2  547698.00 
535279.00  532819.25  530374.62  527936.12 
515871.94  513488.06  511108.56  508754.81 

497029.50  494715.37  492420.44  490106.44 
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I  10  —  output  specifier 

I  W1  JCV2  TNI  YW2  —  Quantex  window  parameters 
I  NL  NU  NA  --  Sean  range; No.  of  flalda  averaged 
I  SCALE  —  aoaling  parameter  (mieroradiana/plxel) 
I  OOPT  —  diameter  of  option  (metera) 

I  OBS  —  obaeuration  ratlo- 
I  WAV  —  wavelength  of  laaer  light  (metera) 
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APPENDIX  D.  SOURCE  CODE  LISTINGS  FOR  DATA  ANALYSIS  PROGRAMS 


This  appendix  contains  the  source  code  listings  of  the  main  modules  of 
MTFLTV,  MTFLD,  and  PREDICT,  of  each  secondary  module,  and  of  the  command  file 
GOMP&LINK .  CMD . 


ooo  non  nnn  on  non 


FTN7X.L.Q 
$FIL£S(0 ,2) 

PROGRAM  MTFLTV 
C 

C  This  program  calculates  the  path  weighted  Cn  parameters  from  short  term, 

C  long  term  and  wander  data  observed  by  a  vidicon  looking  through 

C  a  telescope  directly  at  a  HeNe  laser  source  some  distance  away. 

C 

0  DLTA(512)*  The  data  used  for  long  term  calculations 
C  DSTAC512)*  The  data  used  for  short  term  calculations 

C  DLTK512)*  The  calculated  diffraction  function  of  the  designator  optics 
C  DSTK512)*  The  source  distribution  funciton  of  the  designator 
C  DWD(512)*  The  wander  data 
C  DWDI ( 5 1 2 ) s  *••**  not  presently  used  •••*» 

C  WAVs  Wavelength  of  the  designator  radiation 
C  OOPTs  Diameter  of  the  designator  optics 
C  0BS»  Obscuration  ratio  of  the  designator  optics 
C  SCALEs  Scale  of  viewing  telescope  in  microradians  per  TV  line 
C  SSIGMAs  •*»*•  not  used  in  this  program  •*•** 

C  RANGEs  The  distance  from  the  designator  to  target 
C  KSs  A  parameter  that  indicates  the  highest  term  in  the  calculated 
C  data  that  has  significance. 

C  1DATE ( 20 ) s  Information  on  the  time  and  date  of  run 
C  ICMTS(20)s  Comments  about  the  run 

C  IWEATH(20)s  Comments  about  the  weather  at  time  of  run 
C  LUs  Logical  unit  number  for  output  terminal 
C  Lis  Logical  unit  number  for  input  terminal 
C  NAs  Total  number  of  video  frames  averaged 

C  NLs  The  smallest  pixel  number  in  a  scan  line  that  is  summed 

C  NCs  The  largest  pixel  number  in  a  scan  line  that  is  summed 

C  10s  Outputs  the  data  and  all  calculations  to  files  30  to  50 

C 

COMMON  /PARAMS/WAV , DOPT , OBS , SCALE , RANGE , EXT , A , SSIGMA 
COMMON  /WINDOW/XW 1 , XW2 , TW 1 , YW2 , NFLDS 

Input  general  data  about  the  run 

CALL  GINPUT( 1 ,10) 

Input  data  from  the  video  digitizer  and  set  up  the  long  term, 
short  term  and  wander  data  blocks. 

Currently,  the  number  of  data  points  is  set  to  256. 

NPTSs256 

CALL  DINPUT(NPTS,DLTA,DSTA,DWD) 

IF(IO.EQ.I)  CALL  OUTPUT ( 1 ,NPTS,DLTA) 

IFdO.EQ.  1 )  CALL  OUTPUT ( 2 ,  NPTS ,  DSTA ) 

IF(IO.EQ.I)  CALL  OUTPUT < 18, NFLDS, DWD) 

Center  the  long  and  short  term  line  spread  functions 

NsNPTS 

CALL  LSFS(N.DLTA) 

IF(IO.EQ.I)  CALL  OUTPUT ( 3, N.DLTA) 

NsNPTS 

CALL  LSFS(N.DSTA) 

IFdO.EQ. 1)  CALL  0UTPUT(4,N,DSTA) 

Take  the  fast  Fourier  transforms  of  the  long  and  short  term 
line  spread  functions  to  obtain  the  optical  system  MTFs 


oor>  non  nn  nn  non nn  on  no  o o o  On  oo 


NsNPTS  • 

CALL  FFT ( 0 , M , OLTA  > 

IF(IO.EQ.I)  CALL  OUTPUT (5, N.DLTA) 

NsNPTS 

CALL  FFT(0,N,DSTA) 

IF(IO.EQ.I)  CALL  0UTPUT(6,N,DSTA) 

Calculate  the  diffraction  pattern,  the  line  spread  function, 
and  the  MTF  of  the  telescope  optics 

NaNPTS 

CALL  OPTICS (N, DLTI) 

IF(IO.EQ.I)  CALL  OUTPUT ( 10, N, DLTI) 

CALL  LSF ( N , OLTI ) 

IF(IO.EQ.I)  CALL  OUTPUT ( 1 1 , M , DLTI ) 

NaMPTS 

CALL  FFT(0,N,DLTI) 

IFCI0.EQ.1)  CALL  OUTPUTC 12, N, DLTI) 

Obtain  the  MTF  of  the  TV  system 

NsNPTS 

CALL  TVMTF ( N , OSTI ) 

IF(IO.EQ.I)  CALL  OUTPUT ( 14, N.DSTI) 

Multiply  the  telesoope  optics  MTF  by  the  TV  system  MTF  to 
obtain  the  combined  MTF 

CALL  MUL ( N , OSTI , DLTI ) 

IF(IO.EQ.I)  CALL  OUTPUT ( 15, N.DSTI) 

Divide  the  overall  system  MTFs  for  long  and  short  term  data  by 
the  oomblned  telesoope  optics  and  TV  system  MTF  to  obtain  the 
MTF  of  the  atmosphere 

CALL  DI VIDE ( N , DSTA , DSTI ) 

IF(IO.EQ.I)  CALL  0UTPUT(16,N,DSTA) 

CALL  DIVIDE (N.DLTA, DSTI) 

IF(IO.EQ.I)  CALL  OUTPUT ( 17, N.DLTA) 

Calculate  the  Cn*#2  values  from  the  long  and  short  term  data 
by  fitting  the  atmospheric  MTFs  to  Fried's  model 

CALL  FIT ( 0 , CNSQLT, COEFFLT , N , DLTA ) 

CALL  FIT ( 1 , CNSQST , COEFFST , N , DSTA) 

Caloulate  Cn**2  from  the  laser  spot  wander  data 

CALL  HANDER < CNSQWA , NFLDS , DWD ) 

Output  the  Cn**2  values  to  a  file  and  to  the  operator's  terminal 

0PEN(905,FILE> '0UTPUTCM.DAT' ) 

WRITE(905 ,100)  CNSQST, CNSQLT, CNSQWA 
CLOSE (905) 

100  F0RMAT('  Cn**2  measured  from  short  term  *  * , E 1 1 . 4 , / , 

#  '  Cn**2  measured  from  long  term  *  '.E11.4,/, 

#  '  Cn**2  measured  from  wander  *  ' , E 1 1 . 4 ) 

SCALE* 10.0*SCALE 


no  n 


Recalculate  the  long  and  short  tern  atmospheric  MTFs 
using  the  values  of  Cn**2  just  found  • 

CALL  AMTF  <  0 , CNSQLT , COEFFLT , NPTS , DLTA  > 

CALL  AMTF( 1 , CNSQST , COEFFST , NPTS , DSTA) 

XF(IO.EQ. 1 )  CALL  OUTPUT ( 20 , NPTS, DLTA) 

IF(IO.EQ.I)  CALL  OUTPUT(21 , NPTS, DSTA) 

STOP 
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$FILES(0,2> 

PROGRAM  MTFLD 
C 

C  This  program  calculates  the  path  weighted  Cn  parameters  from  short  term, 

C  long  term  and  wander  data  observed  by  a  vidicon  looking  at  the 

C  spot  pattern  produced  by  a  laser  designator 
C 

C  DLTA(512)«  The  data  used  for  long  term  calculations 
C  DSTA(512)»  The  data  used  for  short  term  calculations 

C  DLTK512)*  The  calculated  diffraction  function  of  the  designator  optics 
C  DSTI(512)>  The  source  distribution  funclton  of  the  designator 
C  DWD(512)s  The  wander  data 
C  DWDH512)*  ***••  not  presently  used  *•**• 

C  WAVs  Wavelength  of  the  designator  radiation 
C  DOPT*  Diameter  of  the  designator  optics 
C  OBS*  Obscuration  ratio  of  the  designator  optics 
C  SCALEs  Scale  of  viewing  telescope  in  microradians  per  TV  line 
C  SSIGMAs  •*»*»  not  used  in  this  program  »»*»» 

C  RANGEs  The  distance  from  the  designator  to  target 
C  KSs  a  parameter  that  indicates  the  highest  term  in  the  calculated 
C  data  that  has  significance 

C  IDATE(20)s  Information  on  the  time  and  date  of  run 
C  ICMTS(20)s  Comments  about  the  run 

C  IWEATH(20)s  Comments  about  the  weather  at  time  of  run 
C  LUs  Logical  unit  number  for  output  terminal 
C  Lis  Logical  unit  number  for  input  terminal 
C  NAs  Total  number  of  video  frames  averaged 

C  NLs  The  smallest  pixel  number  in  a  scan  line  that  is  summed 

C  NCs  The  largest  pixel  number  in  a  scan  line  that  is  summed 

C  IQs  Outputs  the  data  and  all  calculations  to  files  30  to  50 

C 

COMMON  /PARAMS/ WAV , DOPT , OBS , SCALE , RANGE , EXT , A , SSIGMA 
COMMON  /WINDOW/XW 1 , XW2 ,  YW 1 , YW2 , NFLDS 

Input  general  data  about  the  run 

CALL  GINPUT(2,I0) 

Input  data  from  the  video  digitizer  and  set  up  the  long  term, 
short  term  and  wander  data  blocks. 

Currently,  the  number  of  data  points  is  set  to  256. 

NPTSs256 

CALL  DINPUT ( NPTS , DLTA , DSTA , DWD) 

IF(IO.EQ.t)  CALL  0UTPUT( 1 , NPTS, DLTA) 

IF(IO.EQ.I)  CALL  OUTPUT (2, NPTS, DSTA) 

IF(IO.EQ.I)  CALL  0UTPUT( 18, NFLDS, DWD) 

Center  the  long  and  short  term  line  spread  functions 

NaNPTS 

CALL  LSFS(N.DLTA) 

IF(IO.EQ.I)  CALL  OUTPUT ( 3 , N , DLTA ) 

NsNPTS 

CALL  LSFS(N ,DSTA) 

IFCIO.EQ.I)  CALL  OUTPUT ( 4, N, DSTA) 

Take  the  fast  Fourier  transforms  of  the  long  and  short  term 
line  spread  functions  to  obtain  the  optical  system  MTFs 


NsNPTS 

CALL  FFTCO.N.DLTA) 

IF(IO.EQ.I)  CALL  OUTPUTC5.N.DLTA) 

NsNPTS 

CALL  FFT(0,N,DSTA) 

IFdO.EQ.  1)  CALL  OUTPUT C  6 ,  N , DSTA) 

Calculate  the  laser  source  function  assuming  it  has  a  Gaussian 
energy  distribution 

NsNPTS 

CALL  SOURCE (N, DSTI) 

IFdO.EQ. 1)  CALL  OUTPUTC  7,  N,  DSTI) 

Calculate  the  line  spread  function  and  the  MTF  of  the  source 
CALL  LSFCN.DSTI) 

IFdO.EQ.  1)  CALL  OUTPUT C  8 ,  N , DSTI ) 

NsNPTS 

CALL  FFTCO.N.DSTI) 

IFdO.EQ.  1)  CALL  OUTPUT(9,N,DSTI) 

I 

Calculate  the  diffraction  pattern,  the  line  spread  function, 
and  the  MTF  of  the  telescope  optics 

NsNPTS 

CALL  OPTICS (N.DLTI) 

IFdO.EQ.  1)  CALL  OUTPUTC 10, N.DLTI) 

CALL  LSF( N.DLTI) 

IFdO.EQ.  1)  CALL  OUTPUT(  1 1 , N.DLTI) 

NsNPTS 

CALL  FFT(0, N.DLTI) 

IFdO.EQ.  1)  CALL  OUTPUT ( 12, N.DLTI) 

Multiply  the  source  and  optics  MTFs  to  obtain  their  combined  MTF 

CALL  NUL(N.DSTI.DLTI) 

IFdO.EQ.  1)  CALL  OUTPUT ( 13, N, DSTI) 

Obtain  the  MTF  of  the  TV  system 

CALL  TVMTF( N.DLTI) 

IFdO.EQ. 1)  CALL  OUTPUT (  1 4 , N , DLTI) 

Multiply  the  source-optics  MTF  by  the  TV  system  MTF  to 
obtain  their  combined  MTF 

CALL  MULCN, DSTI, DLTI) 

IFdO.EQ.  1)  CALL  OUTPUTC  15, N, DSTI) 

Divide  the  overall  system  MTFs  for  long  and  short  term  data  by 
the  combined  source-optlcs-TV  system  MTF  to  obtain  the 
MTF  of  the  atmosphere 

CALL  DIVIDECN, DSTA, DSTI) 

IFdO.EQ. 1)  CALL  OUTPUTC  1 6 ,  N , DSTA) 

CALL  DI VIDE CN.DLTA, DSTI) 

IFCIO.EQ. 1 )  CALL  OUTPUTC 17, N.DLTA) 

Calculate  the  Cn**2  values  from  the  long  and  short  term  data 
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by  fitting  the  ataospheric  MTFs  to  Fried's  model 

CALL  FIT ( 0 , CNSQLT , COEFFLT , N , DLT A ) 

CALL  FIT( 1 , CNSQST , COEFFST , N , DSTA) 

Calculate  Cn**2  from  the  laser  spot  wander  data 

CALL  WANDER ( CNSQWA , NFLDS , DWD) 

Output  the  Cn**2  values  to  a  file  and  to  the  operator's  terminal 

OPEN ( 905 , FILE* • OUTPUTCN . DAT • ) 

WRITE<905, 100)  CNSQST, CNSQLT, CNSQWA 
CLOSE (905) 

100  FORMAT ( *  Cn**2  measured  from  short  term  *  ',E11.4,/, 

#  '  Cn**2  measured  from  long  term  *  ',E11.4,/, 

#  '  Cn**2  measured  from  wander  *  ',E11.4) 

SCALE* 10.0*SCALE 

Recalculate  the  long  and  short  term  atmospheric  MTFs 
_  using  the  values  of  Cn**2  just  found 

CALL  AMTF ( 0 , CNSQLT , COEFFLT , NPTS , DLTA ) 

CALL  AMTF( 1, CNSQST, COEFFST, NPTS, DSTA) 

IF(IO.EQ.I)  CALL  0(JTPUT<  20,  NPTS,  DLTA) 

IF(IO.EQ.I)  CALL  0UTPUT(21 , NPTS, DSTA) 

STOP 

END 
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|FILES(0,2) 

PROGRAM  PREDICT 
C 

C  Predicts  the  performance  of  an  optical  projection  system  in 
C  a  turbulent  atmosphere. 

C 

COMMON  /PARAMS/WAV , DOPT , OBS , SC ALE , RANGE ,EXT,A, SSIGMA 
.COMMON  /PPARMS/LSOURCE.LOPTICS.CNSQ 
REAL  FN1(512),FN2(512) 

DATA  IU,IOU/5,8/ 

C 

C  Input  general  program  parameters 

C 

LNORMsO 

NPTS*256 

CALL  GINPUT(3,IO) 

IF(LSO(JRCE  .ME.  3)  THEN 

IF(LSOURCE  .EQ.  1)  THEN 
NsNPTS 

CALL  SOURCE (N,FN1) 

IF (10  .EQ.  1)  CALL  POUTPUT ( 1 , N , FN 1 ) 

ELSE  IF(LSOURCE  .EQ.  2)  THEN 
0PEN(901 , FILE= 'SOURCE. DAT' ) 

READ (901  ,*)  (FN1(I),  Isl.NPTS) 
a0SE(901) 

ENDIF 

C 

C  Convert  PSF  of  source  to  LSF 
C 

CALL  LSF(NPTS.FNI) 

IF (10  .EQ.  1)  CALL  P0UTPUT(2,NPTS,FN1 ) 

C 

C  Calculate  the  source  MTF 
C 

NsNPTS 

CALL  FFT ( 0 , N , FN 1 ) 

IF (10  .EQ.  1)  CALL  POUTPUT ( 3, N.FN1 ) 

IF(LOPTICS  .EQ.  1)  THEN 
C 

C  Calculate  the  optics  diffraction  pattern  if  not  input  from  file 
C 

NsNPTS 

CALL  OPTICS(N,FN2) 

IF(I0  .EQ.  1)  CALL  POUTPUT ( 4 f N , FN2 ) 

ELSE  IF(LOPTICS  .EQ.  2)  THEN 
C 

C  Input  optics  diffraction  pattern  from  data  file  if  not  calculated 
C 

OPEN (901, FILE* ’OPTICS.DAT’) 

READ (901 ,  •)  (FN2(I),  I*1,NPTS) 

CLOSE (901) 

ENDIF 

C 

C  Convert  diffraction  pattern  (PSF)  of  optics  to  an  LSF 
C 

CALL  LSF(NPTS,FN2) 

IF( 10  .EQ.  1)  CALL  P0UTPUT(5 ,NPTS, FN2) 

C 

C  Calculate  the  optics  MTF 
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NsNPTS 

CALL  FFT(0,N,FN2) 

Normalize  the  optics  MTF 

DIV*FN2( 1 ) 

00  1*1, N 

FN2(I)sFN2(I)/DIV 
END  00 

IFUO  .EQ.  1)  CALL  POUTPUT ( 6, N,FN2) 

Multiply  the  source  and  optics  MTFs  to  obtain  their  combined  MTF 

CALL  MUL(N,FN1,FN2) 

IF (10  .EQ.  1)  CALL  P0UTPUT ( 7 . N , FN 1) 

ELSE 

Input  source-optics  MTF 

.  0PEN(901 .FILE* ’SRCA0PT.DAT' ) 

READC901 ,  •)  (FNKI),  1*1, NPTS) 

CLOSE (901) 

EN0IF 

Calculate  the  atmospheric  MTF  —  long  term 

CALL  AMTF(0,SQCN,-1 . , NPTS,FN2) 

IF (10  .EQ.  1)  CALL  POUTPUT (8, NPTS, FN2) 

Multiply  the  source-optics  MTF  by  the  atmospheric  MTF  to  obtain 
the  overall  MTF  —  long  term 

CALL  MUL ( NPTS , FN2 , FN 1 ) 

IF(I0  .EQ.  1)  CALL  POUTPUT (9, NPTS, FN2) 

Calculate  the  Image  LSF  by  taking  the  inverse  Fourier  transform 
of  the  overall  MTF  —  long  term 

NaNPTS 

CALL  FFT( 1 , N,FN2) 

IF (10  .EQ.  1)  CALL  P0UTPUT( 10,N,FN2) 

Convert  the  image  LSF  to  a  PSF  —  long  term 

CALL  ABEL(N,FN2) 

CALL  POUTPUT ( 1 1 ,  N , FN2) 

Calculate  the  radially  integrated  power  —  long  term 

CALL  INT(LN0RM, N,FN2) 

CALL  POUTPUT (1 2, N,FN2) 

Calculate  the  atmospheric  MTF  —  short  term 

CALL  AMTF(1,SQCN,-1. ,NPTS,FN2) 

IF(I0  .EQ.  1)  CALL  P0UTPUT( 1 3, NPTS, FN2) 

Multiply  the  source-optics  MTF  by  the  atmospheric  MTF  to  obtain 
the  overall  MTF  —  short  term 
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CALL  MUL (  NPTS , FN2 ,  FN 1 ) 

IF(IO  .EQ.  1)  CALL  POUTPUT( 1 4,NPTS,FN2) 

Calculate  the  Image  LSF  by  taking  the  Inverse  Fourier  transform 
of  the  overall  NTF  —  short  term 

NaNPTS 

CALL  FFT(1,N,FN2) 

IF (10  .EQ.  1)  CALL  P0UTPUT( 15.N.FN2) 

Convert  the  image  LSF  to  a  PSF  —  short  term 

CALL  ABEL(N,FN2) 

CALL  POUTPUT ( 1 6 ,  N , FN2) 

Calculate  the  radially  integrated  power  --  short  term 

CALL  INT ( LNORM , N , FN2 ) 

CALL  P0UTPUT( 1 7 *  M , FM2 ) 

Output  parameters  to  data  file 

OPEN (901 , FILE* ' P ARAMS . DAT • ) 

WRITE (901,1000)  WAV , D0PT , 0BS , SCALE , RANGE , CNSQ , EXT , A , SSIGMA 
1000  FORMAT (// , '  Wavelength  a  ' ,E12.5,/,'  Diameter  of  optics  a 

#  E12.5,/, '  Obscuration  ratio  a  ',E12.5,/,'  Scale  a  ',E12.5, 

#  '  Range  a  • tE12.5,/,'  Cn»»2  a  *,£12.5/,’  Extinction  coe', 

#  'fflclent  a  '.£12.5,/,'  Amplitude  a  ',£12.5,/,'  Sigma  a 

#  E12.5) 

CLOSE (901) 

STOP 

END 
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SUBROUTINE  ABEL(NPTS.FN) 

Converts  a  line  spread  function  to  a  radially  symmetric 
point  spread  function 

REA L  FN(NPTS) 

A*FN(1) 

Ba  1 . E-6*  A 

FN( 1)a.2*(FN(1 )-FN(2) ) 

DO  lal , NPTS-2 
CaFN(I*1) 

FN(I*1 )* . 4*A*.2*C-. 6*FN( 1*2) 

AaC 

IF  (ABS(A)  .GT.  B)  N»I 
END  DO 
DO  I*1,N 

FN(I)aFN(I)/(2 .*SQRT( (FLOAT( I)-. 9)**2-CFL0AT(I-1 ) ) *  *2 ) ) 
11*1*1 
DO  JaII,N 

FN(I)*FN(I)*FN( J)/SQRT( (FLOAT( J)-. 9)**2-(FLOAT(I-1 ) )**2) 
END  DO 

FN(I)aFN(I)/3. 14159 
END  DO 
RETURN 
END 
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REAL  FUNCTION  AIRY(X) 

C 

C  Calculate  the  value  of  the  Airy  function  of  X. 

£  AIRY(X)  a  2*J1(X)/X  where  JKX)  is  the  Bessel  function  of  order  one. 

IF(X  .EQ.  0.)  THEN 
AIRYal. 

RETURN 

ENDIF 

Z»0 

IF(X  .LE.  15.)  THEN 
J*20.+10.»X«X»X/3. 

ELSE 

J*90.+X/2. 

ENDIF 

IF(X  .GE.  5.)  THEN 
Hal .4*X+60./X 

ELSE 

WaX+6 . 

ENDIF 

IaMAXI (H, (3.+X/4. ) ) 

DO  MsI,J,3 
V.1.E-28 
SaO. 

TaS 

Fa 1 . -4 . •<  FL0AT<  M> /2 . -FL0AT( M/2 ) ) 

MMaM-2 
DO  Kal.MM 

U*2. »FL0AT(M-K) »V/X-T 

TaV 

VaU 

IF((M-K-2)  .EQ.  0)  BaU 
Fa-F 
FFal.-F 
SaS+V*FF 
END  DO 
Ua2.»V/X-T 
SaS+U 
BaB/S 

IF((ABS(B-Z)-ABS<B»1.£-6)>  .LT.  0.)  GOTO  10 
ZaB 
END  DO 

10  AIRY»2.»B/X 
RETURN 
END 
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SUBROUTINE  AMTF(R,EX, J) 

Calculate  the  long  and  short  term  atmospheric  MTF. 

The  result  is  then  multiplied  by  the  extinction  factor  EXP(-EXT»RANGE). 
DIMENSION  R(512) 

COMMON/PARAMS/WAV , DOPT , OBS , SCALE , RANGE , SQCN , EXT , A , B 
AAa21.6*SQCN»RANGE»WAV**(-. 3333333) 

IF  (EX  .LE.'O.)  EX* EXP ( -EXT*  RANGE ) 

D*0. 

R(1)*£X 
L»0 

DO  1*2,512 

IFU)  10,5,15 

5  FF*976. 56* (1-1) /SCALE 

ARGaAA*(FF**1. 666667) 

1F( J.EQ. 1 )  ARG*ARG*(1-(FF»WAV/0PT)*«. 3333333) 

IF( ARG  .GT.  13.)  Lai 
IF( (ARG-D)  .LT.  0.)  La-1 
R ( I ) a  EX*EXP ( -ARG ) 

DaARG 
Mai 

GO  TO  20 

0  R(I)aR(M)»(COS((I-M)*3.11H59/M))»*2 

IF(I.LT.3*M/2)  GOTO  20 
Lai 

15  R(I)aO. 

20  END  DO 
RETURN 
END 
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SUBROUTINE  BITINV(K,L,M) 

Finds  the  bit  inversion  of  an  (M)  bit  number  (K)  to  (L) 
This  is  used  in  the  Fourier  transform  subroutine. 

L*0 

DO  1*0, M-1 

IF(8TEST(K,I) )  THEN  L*IBSET(L,M-I) 

END  DO 
RETURN 
END 
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SUBr  mjtIHE  DIGITIZE 

Digitize  a  frame  of  video  with  the  Quantex  DS-30  (Store  Input  function) 

E  INTEGER  CODE (6) 

DATA  CODE/OOOOO IB, 03 1040B, 00 1056B , 033500B, 007452B , Q30500B/ 

DATA  LUDS30/24/ 

CALL  EXEC (2,21 00B+LUDS30 , CODE , - 1 1 ) 

RETURN 

END 
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CHARACTER* 16  FUNCTION  DIGITS (N , BASE) 

C 

C  Convert  integer  N  to  string  of  digits  in  base  BASE  (<  17) 

C 

INTEGER  N,BASE 

CHARACTER  DIGIT* 16, A*  16 

DATA  DIGIT/'0123456789ABCDEF»/ 

IF(N  .EQ.  -32768)  THEN 
HaO 
ELSE 

HsIABS(N) 

ENDIF 

C 

C  Determine  last  digit 

C 

La HOD (H, BASE) ^1 
A  < 1 : 1 )*DIGIT(L:L) 

Ha  M/ BASE 

IF(N  .LT.  0)  M»M>2*( 16384/BASE) 
la  1 

C  . 

C  Determine  rest  of  digits 

C 

DO  WHILE  (H  .NE.  0) 

IaI+1 

L*M0D(M, BASE)+1 
A(I:I)aDIGIT(L:L) 

HaH/BASE 
END  DO 
C 

C  Invert  string  A 

C 

KaI+1 

DIGITSa’  • 

DO  Ja1,I 
KaK-1 

DIGITSC J: J)aA(K:K) 

END  DO 
RETURN 
END 
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SUBROUTINE  OINPUT( NPTS , DLTA , DSTA , OVD) 

C 

C  This  subroutine  inputs  data  from  the  video  digitizer  and  processes 

C  it  to  provide  the  long  term,  short  term  and  wander  data  for 

C  further  processing. 

C 

COMMON  /WIND0W/XW1 ,XW2, YW1 , YW2 , NFLDS 

COMMON  /PARAMS/WAV , DOPT , OBS , SCALE , RANGE , EXT , A , SSIGMR 

REAL  DLTA (  NPTS ) ,DSTA( NPTS ) , DWD ( NFLDS ) , DLTI ( 5 1 2 ) , DSTI (512) 

INTEGER  XW1 , XW2.YW1 ,YW2,BEEP 
DATA  BEEP/003407B/ 

C 

C  Initialize  arrays 

C 

N2.NPTS/2 
DO  I*  \  NPTS 
DLTA(I)aO. 

DSTA ( I ) *0. 

DLTI(I)aO. 

OSTI(I)aO. 

END  DO 

DO  la 1, NFLDS 
DWD(I)aO. 

END  DO 
C 

C  Perform  summations  of  video  data  over  NFLDS  fields 

C 

DO  la  1, NFLDS 
VINTaO. 

BASEaO. 

WRITE ( 1 , * ( A2) ' )  BEEP 
WRITE ( 1 ,•) 

WRITE  ( 1  ,*)  •  Advance  to  field  #' ,1, '  out  of, NFLDS 
WRITE ( 1  ,•)  *  Enter  1  when  ready  or  0  to  stop:  ' 

READO,*)  INPUT 
IF (INPUT  .EQ.  0)  STOP 
C 

C  Digitize  one  video  field 

C 

CALL  DIGITIZE 
C 

C  Obtain  vertical  video  intensity  profile  by  integrating 

C  horizontal  scan  lines 

C 

DO  Ja 1 , NPTS 

CALL  VIDEODATA (J,Z) 

DLTI(J)aZ 

VINTaVINT+Z 

IF( J  .LE.  32  .OR.  J  .GE.  NPTS-32)  BASEaBASE+Z 
END  DO 

BASEaBASE/64. 

VINTaVINT-BASE»NPTS 

SUMaO. 

DO  Ja 1 , NPTS  • 

SUMaDLTI ( J ) +SUM-BASE 
IFCSUM  .GT.  VINT/2.)  GOTO  10 
END  DO 
JaN2 

10  MIDPTaJ 
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’ SUBROUTINE  DI VIDE ( N, ARRAY 1 , ARRAY2) 

Divide  the  first  N  points  of  ARRAY 1  by  ARRAY2  and  store  the 
result  in  ARRAY  1.  If  the  value  of  the  divisor  is  zero, 
then  the  result  of  the  division  is  set  to  zero. 

DIMENSION  ARRAY 1 <N) ,ARRAY2(N) 

DO  1*1 ,N 

IF(ARRAY2(I>  .NE.  0.)  THEN 

ARRAY 1 ( 1 ) aARRAY 1 ( I ) / ARRAY2 ( D 

ELSE 

ARRAYKDsO. 

ENDIF 
END  DO 
RETURN 
END 
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SUBROUTINE  FFT(LINV,NPTS,FN) 


Fourier  transform  subroutine 

IMPLICIT  COMPLEX  (C) 

COMPLEX  CE(512) , CG( 1024) 

REAL  FN(NPTS) 

DATA  1ST, PI/0, 3. 1415927/ 

NN>2*NPTS+2 

NEXP»NINT(ALOG(FLOAT(NPTS) )/AL0G(2. ) )+1 
TH*PI/NPTS 
CG( 1  )*FN( 1 ) 

DO  I*2,NPTS 

CG(I)»FN(I) 

«  CG(NN-I)aCGd) 

END  DO 

CG(NPTS+1)sCG(NPTS) 

IF (1ST  .EQ.  0)  THEN 
ISTs  1 
A*0. 

nn  T*  i  i|pt<? 

CE(I)*CEXP(CMPLX(0. ,A>) 

AaA+TH 

CE?NPTS/2+1}*(0.,1.) 

END  IF 

DO  Msl.NEXP 

Ka2**(NEXP-M) 

La2**(M-1) 

DO  LLal.L 
KCaLL-1 

CALL  BITINV  <KC,N,N£XP-1) 
CDsCE(N+1 ) 

IF(LINV.EQ. 1)  CDaCONJG(CD) 

KLa2*KC»K+1 

KUsKL+K-1 

DO  IsKL.KU 

CAaCG(I)+CG(I*K)»CD 
CG(I+K)aCG(I)-CG(I+K)*CD 
CG(I)sCA 
END  DO 
END  DO 
END  DO 
NNa2»NPTS 
DO  Ial.NN 
JaI-1 

CALL  BITINV  (J,N,NEXP) 

IF( J  .LT.  N)  THEN 
NaNf  1 

CAsCG(I)/(32.,0.) 

CG(I)aCGCN) 

CG(N)aCA 

ELSE 

IF( J  .EQ.  N)  CG(I)sCG(I)/(32. ,0. ) 
ENDIF 
END  DO 
DO  I*1,NPTS 

FN(I)sCG(I) 

IF(FN(I)  .GE.  1 .E-6*FN( 1 ) )  LASTPTsI 
END  DO 
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RETURN 

END 
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’SUBROUTINE  FIT ( CNSQ , COEFF , NPTS , MTF ) 

C 

C  Find  function  parameters  (Cn**2  and  MTF  coefficient)  by  performing 
C  a  least  squares  fit  of  the  MTF  to  Fried's  model: 

C  MTF(X)  *  COEFF  •  EXP(CNSQ«X) 

C 

COMMON/PARAMS/ WAV , DOPT , OBS , SCALE , R ANGE , EXT , A , B 
REAL  MTF (NPTS) 

SUMXeO. 

SUMTaO. 

SUMX2«0. 

SUMXTeO. 

POINTS. 0. 

DO  I* 1, NPTS 

IF(MTFCI)  .LE.  0.)  GOTO  10 
FF*976 . 56»FLOAT< I- 1 ) /SCALE 
X*-21 .6*FF**( 1 .666667>*RANGE#WAV##(-. 3333333) 

T*ALOG(MTF(I) ) 

POINTS.  POINTS-*’ 1 . 

SUMX. SUMX+X 
SUMX2*SUMX2+X*X 

SUMY.  SUMY-*! 

SUMXYeSUMXYfXn 
IF(I  .GE.  4)  THEN 

IF(MTF(I)  .GE.  MTF (1-1 )  .OR.- MTF(I)  .LT.  .1»NTF(1))  GOTO  10 
ENDIF 
END  00 

10  SQCNs(SUMXY*P0INTS-SUMT*SUMX)/(SUMX2*P0INTS-SUMX*SUMX) 

COEFF . EXP ( ( SUMX2*SUMY -SUMXY*SUMX ) / ( SUMX2*  POINTS-SUMX*SUMX )  ) 

RETURN 

END 
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SUBROUTINE  GINPUT(LPROG, 10) 

This  subroutine  inputs  the  general  data  needed  to 
process  the  run. 

WAV  -  wavelength  of  source  in  micrometers 
SCALE  -  scale  of  the  systea  in  micrometers  per  frame 
line  or  field  line 

DOPT  -  diameter  of  the  objective  of  the  optics  in  meters 
OBS  -  ratio  of  the  diameter  of  the  obscuration  to  the 
diameter  of  the  objective 

SSIGMA  -  the  divergence  of  the  source  in  micrometers 
RANGE  -  the  distance  from  the  optics  to  the  target 
XW1.XW2  -  left  and  right  of  data  window 
TW1,YW2  -  bottom  and  top  of  data  window 

COMMON  /HEADER/HEADER 1 , HEA0ER2 , D ATE 
COMMON  /WINDOW/XW 1 , XW2 , YW 1 , YW2 , NFLDS 
COMMON  /PARAMS/WAV, DOPT , OBS , SCALE , RANGE ,EXT,A, SSIGMA 
COMMON  /PPARMS/LSOURC , LOPTICS , CNSQ 
INTEGER  XW1 ,XW2, YW1 , YW2 
CHARACTERS  HEADER  1 ,  HEADER2 ,  DATE 
10*0 

WRITE ( 1 ,  •)  '  Enter  1  to  input  parameters  interactively:' 
READ<1,»)  10 
IF(I0  .EQ.  0)  THEN 
IU*5 

OPEN <IU, FILE* 'INPUT.DAT') 

READdU,*)  10 
READdU,*)  WAV 
WAV* 1 . E-6*WAV 
READdU,  *)  DOPT 
READdU,*)  OBS 
READdU,*)  SCALE 
READdU ,*)  RANGE 
IF(LPR0G  .LE.  2)  THEN 
READdU,  *(A80)')  HEADER1 
READdU,  *(A80)')  HEADER2 
READdU,  '(A80)')  OATS 
READ(IU,»)  XW1 ,XW2, YW2, YW1 
READdU,*)  NFLDS 
ENDIF 

IF(LPR0G  .EQ.  3)  THEN 
READdU,'  (A80)')  HEADER) 

READdU,  '(A80) ')  HEADER2 
READdU,*)  CNSQ 
READdU,*)  EXT 
READdU,*)  A 
READdU,*)  LSOURCE 
READdU,*)  LOPTICS 
ENDIF 

IF(LPR0G  .EQ.  2  .OR.  LPROG  .EQ.  3)  THEN 
READdU,*)  SSIGMA 
ENDIF 
CLOSEdU) 

ELSE 

Interactive  input  follows. 

IU*  1 

WRITEC1,*)'  Enter  1  for  complete  data  output:' 


READ  (III,*)  10 

WRITEd,*)  •  Enter  the  wavelength  in  micrometers :  ' 

READdU,*)  WAV 
WAVsl .E-6*WAV 

WRITEd,*)  •  Enter  diameter  of  the  optica  in  meters:  ' 
READdU,*)  DOPT 

WRITEd,*)'  Enter  obscuration  ratio  of  optics:  ' 

READ (IU,*)  OBS 

WRITEd,*)'  Enter  pixel  resolution  scale  factor  in  ', 

#  'microradians  per  pixel:' 

READdU,*)  SCALE 

.  WRITEd,*)'  Enter  the  range  in  meters:  ' 

READdU,*)  RANGE 
IF(LPR0G  .LE.  2)  THEN 

WRITEd,*)'  Enter  comments  about  run  (for  header):  ' 

READdU,  '<A80>’>  HEADER  1 

WRITEd,*)'  Enter  weather  Information:  ' 

READdU,  '(A80)')  HEADER2 

WRITEd,*)'  Enter  date  and  time  of  run:  ' 

READdU,  '(A80>»)  DATE 

WRITEd,*)'  Enter  left  pixel  limit  of  image  window:  ’ 
READdU,*)  XW1 

WRITEd,*)'  Enter  right  pixel  limit  of  image  window:  • 
READdU,*)  XW2 

WRITEd,*)'  Enter  top  pixel  limit  of  image  window:  ' 
READdU,*)  YW2 

WRITEd,*)'  Enter  bottom  pixel  limit  of  image  window:  • 
READdU,*)  YW1 

WRITEd,*)'  Enter  number  of  fields  integrated:' 

READdU,*)  NFLDS 
END  IF 

IF(LPR0G  .EQ.  3)  THEN 

WRITEd,*)'  Enter  information  about  source:' 

READ( 1 , ' (A80) • )  HEADER 1 

WRITEd,*)'  Enter  Information  about  optics:' 

READ( 1 , ' (A80) ' )  HEADER2 

WRITEd,*)  '  Enter  Cn**2  in  meters**4-2/3) : ' 

READd,*)  CNSQ 

WRITEd,*)  '  Enter  extinction  coefficient  in  meters**(-1 ) : ' 
READd,*)  EXT 

WRITEd,*)  •  Enter  maximum  of  source  in  Watts/steradian: ' 
READd,*)  A 

WRITEd,*)  '  Enter  source  option:  1  —  Gaussian;  2  —  Input; 

#  '  3  --  Source-optics  fn.  input.' 

READd,*)  LSOURCE 

WRITEd,*)  '  Enter  optics  diffraction  option:  1  —  Computed; 

#  '  2  —  Input.' 

READd,*)  LOPTICS 

ENDIF 

IFCLPROG  .EQ.  1  .OR.  LPROG  .EQ.  3)  THEN 
WRITE (IU, 100) 

FORMAT ( '  Enter  sigma  of  the  designator  laser.',/,'  This  ', 

#  'is  the  angle  in  radians  where  the  intensity  falls  off  to',/ 
f  '  1/e  of  the  peak  intensity.  ') 

READdU,*)  SSIGMA 
ENDIF 
ENDIF 
RETURN 
END 
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SUBROUTINE  INT(LNORM,NPTS,FN> 

Given  the  point  spreed  function,  calculate  the  power  inside 
a  circle  as  a  function  of  radius. 

Will  normalize  the  power  function  to  the  total  power  if  LNORMsI. 

COMMON  /PARAMS/ WA V , OPT , OBRA , SCALE , RANGE , EXT , A , B 
REAL  FN(NPTS) 

Calculate  the  power  function 

R2s (3CALE*RANGE ) **2 
FN ( 1 ) * . 7853982E- 1 2«FN ( 1 ) »R2 
DO  I.2.NPTS 

FN(I)*FN(X-1)«.6.28315E-12«FN(I)»R2»FLOAT(I-1) 

END  DO 

Normalize  the  power  function 

IFCLNORM  .EQ.  1)  THEN 
DO  Isl.NPTS 

FN(I)»FN(I)/FN(NPTS) 

END  DO 
ENDIF 
RETURN 
END  ’ 
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SUBROUTINE  LSF(NPTS,FN) 


Convert  a  point  sprtad  function  to  a  line  spread  function. 

REAL  FN(NPTS) 

DO  lal.NPTS 
ZsFN(I) 

J«1 

10  HaSQRT(FL0AT<(I-1)”2*J»J)) 

IlsZFIX(l) 

RlaFLOAT(IR) 

RFsR.RI 

Z*Z+2 . • ( < 1 . -RF ) *FN ( I R+ 1 ) ♦ RF»FN ( IR*2 ) ) 

JsJ+1 

IF(IR  .LT.  NPTS)  GOTO  10 
FN(I)sZ 
END  DO 
RETURN 
END 


1 
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SUBROUTINE  LSFS(NPTS.FN) 

Shift  a  lina  sprtad  function  so  that  it  is  centered  at  the  origin. 

REAL  FN(NPTS) ,T(512) 

N2*NPTS/2 

A«0. 

00  I»1 , NPTS 
TU).FN(I) 

A*FN(I)+A 
END  DO 
3*0. 

DO  1*1 , NPTS 
S*SVT(I) 

IF(S  .GT.  A/2.)  GOTO  10 
END  DO 
I*N2 

10  DO  J«1,N2 

K1.I+J-1 

K2*I-J*1 

IF(K1  .GT.  NPTS)  K1*K1-NPTS 
IF(K2  .LT.  1)  K2*K2fNPTS 
FN( J)*T(K1 )>T(K2) 

FN( J+N2)*0. 

END  DO 
RETURN 
END 


I 
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SUBROUTINE  MUL ( N , ARRAY 1 , ARRAY2 ) 

Calculate  the  product  of  the  first  N  points  of  ARRAY  1  and  ARRAY2 
and  stora  the  result  in  ARRAY1. 

DIMENSION  ARRAY 1 (N) ,ARRAY2(N) 

00  I«1,N 

ARRAY 1 ( I ) * ARRAY 1 ( I ) *ARRAY2(I ) 

END  DO 
RETURN 
END 
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SUBROUTINE  OPTICS (NPTS.FN) 

Calculate  the  optics  diffraction  pattern  (point  spread  function) 

COMMON  //WAV , DOPT , OBS , SCALE , RANGE , SQCN,EXT,A, SSIGMA , IDATE (30), 

#  IHEAD(80) 

REAL  FN(NPTS) 

INTEGERS  IDATE, IHEAD 
OB*OBS*OBS 
OBAs ( 1 . -OB) **2 
0*3. 14159E-6*SCALE*DOPT/WAV 
NsNPTS 
DO  1*1, N 

RsQ*FLOAT ( I- 1 ) 

IF(R  .LE.  60.)  THEN 

FN(I)*((AIRY(R)-0B»AIRT(R«0BS))»»2)/0BA 
NPTSal 

ELSE 

FN(I)*0. 

ENDIF 
END  DO 
RETURN 
END 


n oo  n  o  o 


FTM7X.L.Q 

SUBROUTINE  OUTPUT(LHEAD,NPTS, ARRAY) 

Output  intermediate  and  final  results  to  files 

COMMON  /PARAMS/WAV , DOPT , OBS , SCALE , RANGE , EXT, A , SSIGMA 
COMMON  /HEADER/HEADER  1 , HEADER2 , DATE 
REAL  ARRAT(NPTS) 

CHARACTER  UNIT*5 , FILE*1 1 , DIGITS* 1 6 
CHARACTERS  HEADER  1  ,HEADER2 ,  DATE ,  TITLE 

Specify  TITLE  string,  i.e.  the  title  of  this  output  data 

IFCLHEAD  .EQ.  1)  THEN 

TITLE*'  Optical  System  LSF  —  Long  Tara' 

ELSE  IF(LHEAD  .EQ.  2)  THEN 

TITLE*'  Optical  Systea  LSF  —  Short  Tern' 

ELSE  IFUHEAD  .EQ.  3)  THEN 

TITLEs'  Centered  LSF  —  Long  Tara' 

ELSE  IFUHEAD  .EQ.  4)  THEN 

TITLEs'  Centered  LSF  —  Short  Tera' 

ELSE  IFUHEAD  .EQ.  5)  THEN 

TITLEs'  Optical  Systea  MTF  —  Long  Tera' 

ELSE  IFUHEAD  .EQ.  6)  THEN 

TITLEs*  Optical  Systea  MTF  —  Short  Tera' 

ELSE  IFUHEAD  .EQ.  7)  THEN 

TITLEs'  Source  Distribution  Function' 

ELSE  IFUHEAD  .EQ*  8)  THEN 

TITLEs'  LSF  of  Source  Distribution  Function’ 

ELSE  IFUHEAD  .EQ.  9)  THEN 

TITLEs'  Fourier  Transfora  of  Source  Distribution  Function' 
ELSE  IFUHEAD  .EQ.  10)  THEN 

TITLEs'  Diffraction  Pattern  of  Telescope  Optics' 

ELSE  IF(LHEAD  .EQ.  11)  THEN 

TITLEs'  LSF  of  Telescope  Optics' 

ELSE  IFUHEAD  .EQ.  12)  THEN 

TITLEs'  MTF  of  Telescope  Optics' 

ELSE  IF(LHEAD  .EQ.  13)  THEN 

TITLEs’  Source  and  Optics  Calibration  Function' 

ELSE  IFUHEAD  .EQ.  14)  THEN 
TITLEs'  MTF  of  TV  Systea' 

ELSE  IFUHEAD  .EQ.  15)  THEN 

TITLEs'  MTF  of  the  Overall  System' 

ELSE  IFUHEAD  .EQ.  16)  THEN 

TITLEs'  MTF  of  the  Atmosphere  —  Long  Term' 

ELSE  IFUHEAD  .EQ.  17)  THEN 

TITLEs’  MTF  of  the  Atmosphere  —  Short  Term' 

ELSE  IF(LHEAD  .EQ.  18)  THEN 

TITLE*'  Wander  of  the  Laser  Spot  Image' 

ELSE  IFUHEAD  .EQ.  19)  THEN 

TITLEs’  Fourier  Transform  of  the  Laser  Spot  Image  Wander' 
ELSE  IFUHEAD  .EQ.  20)  THEN 

TITLEs’  MTF  of  the  Atmosphere  Calculated  from  Fitted  Data'// 

#  •  —  Long  Term' 

ELSE  IFUHEAD  .  EQ.  21)  THEN 

TITLEs’  MTF  of  the  Atmosphere  Calculated  from  Fitted  Data’// 

#  '  —  Short  Term' 

ENDIF 

C  Specify  I/O  unit  number  and  file  name 


IOU«900+LHEAD 
UNIT»DIGITS(LHEAD, 10) 

FILE* 'OUTPUT* //UNIT 

Output  data  to  flla 

OPEN(IOU,FILE*FILE) 

WRITE ( IOU, •)  TITLE 
WRITE (IOU,»)  HEADER  1 
WRITE ( IOU, •)  HEADER2 
WRITE(IOU,*)  DATE 

WRITE (IOU, '(8(1X*F9.3))')  (ARRAY(I),  I*1,NPTS) 
CLOSE (IOU) 

RETURN 
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SUBROUTINE  POUTPUTC LHEAD , NPTS , ARRAY) 

Output  Intermediate  and  final  results  to  files 

COMMON  /PARAMS/WAV , DOPT , OBS , SCALE , RANGE , EXT , A , SSIGMA 
COMMON  /HEADER/HEADER  1 , HEADER2 , DATE 
REAL  ARRAY (NPTS) 

CHARACTER  UNIT*5 ,FILE*1 1 , DIGITS* 1 6 
CHARACTER*80  HEADER1 .HEADER2, DATE, TITLE 

Specify  TITLE  string,  i.e.  the  title  of  this  output  data 

IF (LHEAD  .EQ.  1)  THEN 

TITLES'  Gaussian  source  function' 

ELSE  IFCLHEAD  .EQ.  2)  THEN 

TITLEs'  Source  line  spread  function' 

ELSE  IF( LHEAD  .EQ.  3)  THEN 

TITLEs'  MTF  of  the  source' 

ELSE  IF (LHEAD  .EQ.  4)  THEN 

TITLEs'  Optics  diffraction  pattern' 

ELSE  IF( LHEAD  .EQ.  5)  THEN 

TITLEs'  Optics  line  spread  function' 

ELSE  IF (LHEAD  .EQ.  6)  THEN 

TITLEs'  MTF  of  the  optics* 

ELSE  IF (LHEAD  .EQ.  7)  THEN 

TITLEs'  Combined  source  and  optics  MTF' 

ELSE  IF (LHEAD  .EQ.  8)  THEN 

TITLEs'  Atmospheric  MTF  —  long  term' 

ELSE  IF (LHEAD  .EQ.  9)  THEN 

TITLEs*  Combined  source,  optics,  and  atmospheric  MTF  —  •// 

#  'long  term' 

ELSE  IF (LHEAD  .EQ.  10)  THEN' 

TITLEs '  Image  line  spread  function  — -  long  term' 

ELSE  IF (LHEAD  .EQ.  11)  THEN 

TITLEs'  Image  point  spread  function  ~  long  term' 

ELSE  IF (LHEAD  .EQ.  12)  THEN 

TITLEs'  Radially  integrated  power  of  source  —  long  term’ 
ELSE  IF (LHEAD  .EQ.  13)  THEN 

TITLEs'  Atmospheric  MTF  —  short  term' 

ELSE  IF( LHEAD  .EQ.  14)  THEN 

TITLEs'  Combined  source,  optics,  and  atmospheric  MTF  —  '// 

#  'short  term' 

ELSE  IF (LHEAD  .EQ.  15)  THEN 

TITLEs'  Image  line  spread  function  —  short  term’ 

ELSE  IF (LHEAD  .EQ.  16)  THEN 

TITLEs'  Image  point  spread  function  —  short  term' 

ELSE  IF (LHEAD  .EQ.  17)  THEN 

TITLEs'  Radially  integrated  power  of  source  — -  short  term* 
ENDIF 

Specify  I/O  unit  number  and  file  name 

I0Us900+LHEAD 
UNITsDIGITS(LHEAD, 10) 

FILEs 'OUTPUTV/UNIT 

Output  data  to  file 

0PEN( IOU , FILEsFILE) 
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WRITE (IOU,»)  TITLE 
WRITE(IOU,»)  HEADER1 
WRITE (IOU,»)  HEADERS 
WRITE (IOU,»)  DATE 

WRITEdOU.dadX.Fg.B))1)  (ARRAY(I),  I*1,NPTS> 
CLOSE (IOO) 

RETURN 

END 


,L,Q, J, S 

SUBROUTINE  SOURCE(NPTS.FN) 

Calculates  the  values  of  a  Gaussian  source  function  having 
a  aean  of  zero,  a  spread  of  SSIGMA,  and  an  aaplitude 
equal  to  A. 

COMMON  /PARAMS/ WA V, DOPT , OBS , SCALE , RANGE , EXT , A , SSIGMA 
REAL  FN(NPTS) 

DO  Is 1 , NPTS 

X»SCAL£»FL0AT(I-1) 

ARGsX*X/ ( 2  *  *SSIGMA*SSIGMA) 

IFCARG  .LE.  70.)  THEN 
FN ( I ) sA*EXP ( «ARG) 

NPTS- 1 

ELSE 

FN(I)*0. 

ENDIF 
END  DO 
RETURN 
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SUBROUTINE  TVMTF C  NPTS , FN ) 

C 

C  Obtain  ttaa  TV  systen  MTF . 

C 

C  Not*:  For  now  this  routine  creates  a  step  function  with 
C  the  first  256  values  of  FN  set  to  1  and  the 

C  rest  set  to  0. 

C 

REAL  FN(NPTS) 

11*256 

DO  1*1 , NPTS 

XF(I  .LT.  N)  THEN 
FN(I)*1 . 

ELSE 

FN(I)*0. 

ENDIF 
END  DO 
RETURN 
END 
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SUBROUTINE  VIDEODATA ( ROW, ROWSUM) 

Read  in  a  row  of  digitized  video  data  and  sum  the  data 

COMMON  /WINDOW/ XW1,XW2,YW1,YW2,NFLDS 
INTEGER  ZDATC512) , ROW , QROW , QCOL , XW 1 ,XW2,YW1 ,YV2 
INTEGER  C0DE(9) ,L0  ADD, HI  ADD.WCOUNT 
DATA  LUDS 30/24/ ,LUUPIB/207 

DATA  C0DE/00000 IB , 036442B , 0 , 0 , 0 , 0 , 0 , 050 1 40B , 007400B/ 
NPIX.XW2 -XV1+1 
QR0W.R0W+TW2-1 
QCOLsXWI 

LO  ADD«IAND(QCOL-1 ,255) 

Hi"  ADD*2* ( QROW. 1 ) +IAND ( QCOL- 1 , 25  6  > /2S6 
WC0UNT*NPIX+1 

C0DE(3)*001040B+IAND(L0  ADD, 15) 

CODE ( 4 ) *030 1 OOB+I AND ( LO” ADD,240)»16*IAND(HI  ADD,  15) 
CODE(5)*050140B+IAND(HI"ADD,240)»16«.IAND(HI”XDD,3840)/256 
CODE  ( 6 )  *001 440B*IAND<  WCBUNT ,  1 5 ) 

CODE(7)*030100B+IAND(VC0UNT,240)»16+IAND(VC0UNT,3840)/256 
CALL  ABRT(LUHPIB,2) 

CALL  EXEC (2,21 00B+LUDS30 , CODE ,-17) 

CALL  EXEC (1,21 00B+LUDS30 , ZDAT , -NPIX) 

ROWSUMsO. 

DO  1*1 ,NPIX,2 

ROWSUM»ROWSUM*IAND(ISHFT(ZDAT(I) ,0) ,377B)+ 

#  IAND(ISHFT(ZDAT(I),-8),377B) 

END  DO 
RETURN 
END 
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FTN7X,L,Q,J,S 

SUBROUTINE  WANDER (NPTS.FN) 

Calculate  the  Cn*#2  from  the  wander  of  the  laser  spot  image 

COMMON  //WAV , DOPT , OBS , SCALE , RANGE , SQCN , EXT , A , SSIGMA , IDATE ( 30 ) t 
#  IHEAD(SO) 

REAL  FN(NPTS) 

XF(NPTS  .LE.  1)  RETURN 
SUM«0. 

SUMSQsO. 

Calculate  the  sum  of  and  the  sum  of  the  squares  of  the  wander  data 

DO  1*1 , NPTS 

SUM*SUM*FN(I) 

SUMSQ*SUMSQ+FN  ( I )  »FN  ( I) 

END  DO 

Calculate  the  variance  of  the  wander  data 

VAR*(SUMSQ-SUM»SUM/FLOAT(NPTS>  >/FL0AT(NPTS-1 ) 

Calculate  the  Cn**2  using  the  variance  of  the  wander  data 

SQCNsO . 956E-1 2»SCALE»SCALE»VAR*D0PT»» ( 0.3333333 ) /RANGE 
RETURN 
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**  COMP  &L INK. CMD 


aa 

aa 

“  Command  file  for  eompillng  and  linking  programs 

••  MTFLTV.MTFLD  and  PREDICT 

aa 

M 


•  • 

••  Compilt  all  of  tho  program  modulas 

*• 

FTM7X  /MTF/ ABEL. FTN  0  - 
FTH7X  /MTF/AIRT.FTM  0  - 
FTM7X  /MTF/AMTF.FTM  0  - 
FTM7X  /MTF/BITIMV.FTM  0  - 
FTM7X  /MTF/DIGITIZE . FTN  0  - 
FTM7X  /MTF/DIGITS.FTN  0  - 
FTM7X  /MTF/DINPUT.FTN  0  -  . 

FTM7X  /MTF/DI VIDE. FTM  0  - 
FTM7X  /MTF/FFT.FTM  0  - 
FTM7X  /MTF/FIT.FTN  0  - 
FTN7X  /MTF/GIMPUT. FTM  0  - 
FTM7X  /MTF/IMT.FTM  0  - 
FTM7X  /MTF/LSF.FTN  0  - 
FTM7X  /MTF/LSFS.FTM  0  - 
FTM7X  /MTF/MTFLD . FTM  0  - 
FTM7X  /MTF/MTFLTV.FTN  0  - 
FTM7X  /MTF/MUL.FTN  0  - 
FTM7X  /MTF/OPTICS. FTM  0  - 
FTM7X  /MTF/OUTPUT.FTN  0  - 
FTM7X  /MTF/POUTPUT.FTN  0  - 
FTN7X  /MTF/PREDICT. FTM  0  - 
FTM7X  /MTF/SOURCE.FTM  0  - 
FTM7X  /MTF/TVMTF.FTN  0  - 
FTM7X  /MTF/ VIDEODATA. FTM  0  - 
FTM7X  /MTF/ WANDER . FTM  0  - 
•• 

•*  Marga  tha  ralocatabla  coda  filas  into  fila  MTF. REL. 

••  Fila  MERGE.CMD  contains  a  list  of  tha  filas  to  ba  marged. 

•a 

MERGE  /MTF/ MERGE.CMD  /MTF/ MTF. REL 

aa 

**  Craata  an  indaxad  library  of  tha  ralocatabla  coda  modulas 

aa 

LINDX  /MTF/MTF.REL  /MTF/MTF.LIB 

aa 

••  Purga  tha  ralocatabla  coda  filas 

aa 

PU  /MTF/MTF.REL 
PU  /MTF/ ABEL. REL 
PU  /MTF/ AIRY. REL 
PU  /MTF/ AM TF. REL 
PU  /MTF/BITINV. REL 
PU  /MTF/DIGITS . REL 
PU  /MTF/DIGITIZE. REL 
PU  /MTF/DIMPUT. REL 
PU  /MTF/DIVIDE. REL 
PU  /MTF/FFT. REL 
PU  /MTF/FIT. REL 
PU  /MTF/GIMPUT. REL 


PO  /MTF/INT. REL 
PU  /MTF/LSF. REL  . 

PU  /MTF/LSFS. REL 
PU  /MTF/MUL.REL 
PU  /MTF/OPTICS. REL 
PU  /MTF/OUTPUT. REL 
PU  /MTF/POUTPUT. REL 
PU  /NTF/SOURCE. REL 
PU  /MTF/TVMTF . REL 
PU  /MTF/ VIDEODATA. REL 
PU  /MTF/WANDER . REL 
•• 

**  Link  tha  programs 

•• 

LINK  /MTF/MTFLTV . FTN  /MTF/MTF.LIB 
LINK  /MTF/MTFLD . FTN  /MTF/MTF.LIB 
LINK  /MTF/ PREDICT. FTN  /MTF/MTF.LIB 


(THIS  PAGE  INTEFTIOFA* T "  "TA'T) 
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APPENDIX  E.  TEST  PLAN 


1.0  SUMMARY 

1.1  BACKGROUND 

Previous  work  at  the  Naval  Postgraduate  School  describes  a  system  for 
laser  spot  profile  analysis  using  TV  cameras  to  view  the  spot  of  a  laser 
designator  for  the  evaluation  of  the  laser  performance  (ref  1,  app  B).  This 
report  also  describes  a  proposed  field  test  program  for  the  system. 

1.2  OBJECTIVES  OF  TESTS 

The  first  objective  of  this  test  Is  to  compare  Cn^  data  produced  by 
the  NPS  mechanical  scanner  and  a  video  measurement  system  over  a  coincident 
atmospheric  path.  The  second  objective  of  this  test  Is  to  apply  the  C^ 
data  to  the  evaluation  of  a  ground-based  laser  designator. 


2.0  DETAILS  OF  TEST  PLAN 
2.1  EQUIPMENT  SETUP 

The  general  setup  Is  shown  in  figure  1,  which  is  a  block  diagram  of  the 
bptlcs.  At  one  end  of  the  range  would  be  the  designator,  slit-scan  OTF 
measurement  system  of  the  NPS,  and  a  TV  MTF  measurement  system.  The  target 
screen  would  be  placed  about  1  to  3  kilometers  toward  the  other  end  of  the 
range.  The  screen  is  viewed  by  a  silicon  vldicon.  Finally,  there  Is  a  small 
HeNe  gas  laser  behind  the  screen,  shining  through  a  small  hole  in  the  screen. 
This  is  the  source  for  the  two  OTF  measuring  systems. 

In  the  first  step,  field  measurements  would  be  made  with  a  strapped-down 
laser  designator  with  its  spot  trained  onto  a  uniform  target  screen.  The 
designator  must  be  securely  strapped  down  so  that  any  vibrations  caused  by 
wind  or  personnel  walking  about  would  not  deflect  the  designator  beam  by  more 
than  5  microradians.  (One  microradian  maximum  would  be  desirable.)  This  may 
require  mounting  the  designator  on  a  separate  platform.  These  measurements 
would  serve  to  determine  the  stability  of  the  designator  laser  and  optical 
system  by  separating  the  atmospheric  turbulence  effects  from  the  overall 
effects.  The  beam  deflection,  due  to  atmospheric  effects,  Is  of  the  order  of 
20  to  30  microradians  on  a  turbulent  day.  The  telescopes  with  the  slit-scan 
and  TV  must  also  be  strapped  down  and  mounted  on  the  separate  platform. 

This  separate  platform  can  be  constructed  next  to  the  main  platform,  but 
not  actually  attached  to  It.  Thus,  one  can  reach  the  equipment  on  the  plat¬ 
form  to  work  on  It,  but  vibrations  on  the  main  platform  will  not  affect  the 
instruments  on  the  separate  platform. 

For  the  same  reason,  the  target  screen  should  also  be  tied  down  so  that 
the  wind  will  cause  motion  of  no  more  than  approximately  1  cm  motion  maximum. 
The  size  of  the  target  screen  of  8  by  8  feet  should  be  large  enough  to  contain 
the  laser  spot  under  the  must  turbulent  conditions  at  a  distance  of  1  to  3 
kilometers  from  the  designator.  The  reflectivity  of  the  screen  should  be  uni¬ 
form  for  this  part  of  the  test. 
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Figure  1.  Optical  Configuration  for  Field  Test 


A  low  power  (1  mW)  HeNe  laser  will  be  placed  behind  the  middle  of  the 
target  screen  with  a  small  hole  in  the  screen  to  let  the  laser  beam  through. 
This  HeNe  laser  will  be  directed  toward  the  slit-scan  and  TV  system  located 
near  the  designator.  The  HeNe  laser  may  need  a  small  lens  to  diverge  the  beam 
a  few  milliradians  so  that  the  telescope  will  be  uniformly  illuminated. 

The  TV  MTF  system  consists  of  an  8-inch  Celestron  telescope  with  a  barlow 
lens  of  focal  length  -IS  cm  to  increase  the  effective  focal  length  of  the 
telescope.  A  0.6328  micrometer  narrow  band  Interference  filter  will  be  placed 
just  in  front  of  the  vidlcon  between  the  vidicon  and  the  telescope.  An  adap¬ 
tor  must  be  constructed  to  mount  the  TV  onto  the  telescope.  This  adaptor  con¬ 
tains  the  barlow  lens  and  the  interference  filter.  The  dynamic  range  of  a 
vidlcon  tube  is  not  very  large  and  the  light  cannot  be  changed  by  means  of 
stooping  down  the  aperture,  without  affecting  the  OTF  of  the  system.  Thus, 
there  should  be  a  way  of  inserting  neutral  density  filters  in  front  of  the 
vidlcon.  This  system,  excluding  the  adaptor,  will  be  supplied  by  Fort 
Huachuca.  The  adaptor  will  be  made  by  the  Naval  Postgraduate  School. 

A  second  vidicon  is  needed  to  view  the  target  screen.  It  must  have  a 
silicon  focal  plane  in  order  to  be  sensitive  at  1.06  micrometer,  near  Infrared 
(IR),  wavelength  of  the  designator.  It  will  be  placed  near  the  target  screen 
along  the  line  between  the  designator,  but  slightly  off  that  line  to  avoid 
obscuring  the  screen  from  the  designator  location.  This  vidlcon  will  have  a 
1.06  Interference  filter  mounted  in  front  of  the  lens  to  block  out  the  back¬ 
ground  radiation  from  the  sun.  The  synchronizing  signal  for  this  vidlcon  will 
be  used  to  trigger  the  designator  so  that  it  will  fire  during  the  flyback  time 
between  the  frames  of  the  TV. 

The  output  from  each  of  the  two  vidicons  will  be  recorded  on  separate 
Panasonic  TV  tape  recorders  for  later  processing.  The  output  of  the  NPS 
slit-scan  telescope  can  be  processed  immediately.  It  will  also  be  recorded  on 
a  Precision  Instrument  Co.  FM  analog  tape  recorder  for  later  processing. 

2.2  SUPPORT  REQUIREMENTS 

Each  end  of  the  range  will  require  115  Vac.  The  designator  end  will  re¬ 
quire  approximately  20  amperes  while  the  target  end  will  require  approximate¬ 
ly  10  amperes.  The  power  at  the  target  end  could  be  supplied  by  a  small 
gas- powered  generator  set  of  approximately  1  or  2  kW  power. 

The  test  will  require  a  simple  communication  system  between  the  desig¬ 
nator  end  and  the  target  end.  This  could  be  a  field  telephone  or  walky-talky 
radios.  Also,  communication  to  the  laboratory  or  Greely  Hall  would  be  useful. 

The  setting  up  of  the  equipment  and  the  taking  of  data  should  easily  be 
accomplished  by  four  or  five  individuals.  Two  will  be  from  the  Naval  Post¬ 
graduate  School  and  the  others  from  Fort  Huachuca.  During  the  taking  of  data 
two  or  three  people  will  be  needed  at  the  main  platform  to  take  data  on  the 
NPS  equipment  and  the  TV-OTF  system.  The  other  two  people  will  be  at  the  tar¬ 
get  end  operating  the  TV  that  views  the  spot. 


2.3  SCHEDULE  OF  EVENTS 

The  equipment  will  be  set  up  on  the  first  day  and  tested  for  proper  oper¬ 
ation.  This  includes  alignment  and  focusing  the  telescopes.  The  alignment  of 
the  secondary  mirror  of  each  telescope  is  very  important  for  accurate  measure¬ 
ments  of  the  MTF.  The  silicon  vidicon  will  be  focused  on  the  target  screen 
from  a  distance  that  allows  the  target  screen  to  just  fill  the  field  of  view 
of  the  camera.  A  vertical  scale  (this  could  be  a  meterstlck)  is  temporarily 
attached  to  the  target  and  a  recording  is  made  of  the  video  to  give  a  scale  in 
TV  lines  per  meter  of  target  screen.  This  will  be  required  for  later  process¬ 
ing.  If  the  distance  between  the  TV  camera  and  the  screen  is  changed,  this 
measurement  will  have  to  be  repeated.  The  MTF  of  this  vidicon  should  be 
checked  by  placing  a  video  test  pattern  on  the  target  and  measuring  the  qual¬ 
ity  of  the  TV  image. 

Assuming  everything  is  set  up,  aligned,  and  tested  by  evening,  the  first 
measurements  with  the  OTF  systems  will  be  made  at  sunset  which  is  the  period 
of  least  turbulence.  This  will  give  an  approximate  calibration  for  all  sys¬ 
tems  until  a  laboratory-controlled  calibration  can  be  made. 

The  NPS  slit-scan  OTF  system  will  make  measurements  every  half  hour 
during  the  evening  to  get  a  plot  of  the  optical  turbulence  structure  function, 
Cn,  as  a  function  of  time  during  the  evening.  At  the  same  time,  measure¬ 
ments  will  be  made  with  the  TV  MTF  system  and  recorded.  During  a  period  of 
time  when  the  optical  turbulence  is  low,  the  scale  of  the  TV  system  will  be 
determined  by  using  a  grating  in  front  of  the  TV  telescope.  This  calibration 
should  remain  fixed  unless  the  focal  length  of  the  telescope  system  is 
changed. 

It  may  be  necessary  to  use  the  next  day  to  complete  the  equipment  setup 
and  to  make  the  preliminary  tests  described  above. 

One  day  should  be  devoted  to  taking  data  with  the  three  systems.  This 
will  consist  of  taking  a  1-minute  recording  of  the  video  from  each  TV  system 
and  the  NFS  MTF  system  simultaneously.  These  measurements  should  be  made  at 
half-hour  Intervals  throughout  the  day  and  evening.  This  is  to  demonstrate  . 
that  the  three  systems  are  tracking  each  other;  that  is,  they  are  measuring 
essentially  the  same  optical  turbulence  and  not  some  other  systematic  behav¬ 
ior.  If  the  software  for  the  data  processing  is  in  operation  at  this  time, 
some  of  the  data  should  be  processed  for  numerically  checking  the  results  of 
the  three  systems.  However,  a  qualitative  check  can  be  made  at  the  time  of 
each  recording  by  measuring  the  image  spread  directly  on  the  TV  monitors.  The 
spreading  is  almost  proportional  to  Cn.  (It  is  closer  to  Cnl*2.) 

If  smoke  tests  are  to  be  made,  they  should  be  conducted  the  following 
day.  The  smoke  grenades  should  be  set  off  up-wind  at  such  a  position  that  the 
smoke  will  drift  across  the  optical  path  at  several  different  ranges,  e.g., 
0.1,  0.3,  0.5,  0.7,  and  0.9  of  the  distance  from  the  optics  end  to  the  target 
end.  One  series  of  tests  should  be  conducted  during  the  daylight  hours  and 
one  series  of  tests  after  sunset.  Several  minutes  of  recording  of  each  of  the 
three  systems  should  be  made  during  each  of  the  smoke  tests  as  the  smoke 
passes  through  the  optical  path.  Dugway  Proving  Ground  will  provide  the 
equipment  and  personnel  needed  for  the  smoke  test. 
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The  last  day  Is  to  be  used  for  repeating  any  tests  that  seemed  to  give 
unexpected  or  Inconsistent  results.  This  should  allow  time  to  disassemble  the 
equipment  and  pack  it  for  transportation. 

2.4  DATA  ANALYSIS 

The  NPS  data  can  be  processed  at  the  time  the  data  is  taken  to  obtain  a 
Cq  value.  However,  it  will  be  recorded  so  that  it  can  be  reprocessed  at  a 
later  date  if  there  are  any  questions  about  the  results  of  the  other  systems 
and  how  they  compare  with  the  NPS  system. 

The  processing  of  the  data  from  the  two  TV  systems  will  be  similar  but 
not  Identical.  Ten  seconds  of  data  from  the  recording  of  the  0.6328  TV  data 
will  be  transferred  to  the  Eigen  video  disk  recorder.  This  is  300  frames, 
which  is  the  maximum  it  will  hold.  Each  frame  will  be  transferred  to  the 
Quantex  for  digitizing.  The  digital  data  will  then  be  transferred  to  the 
HP-1000,  where  it  will  be  processed  similar  to  the  processing  used  by  the  NPS 
slit-scan  system  to  obtain  a  Cn  value.  Several  consecutive  10-second 
periods  may  be  needed  for  longer  averages,  especially  for  determining  the 
spectrum  of  the  wander. 

The  processing  of  the  data  from  the  1.06  TV  system  will  also  use  the 
Eigen  video  disk  recorder  to  capture  10  seconds  of  video,  or  300  frames. 
However,  the  Quantex  will  subtract  one  frame  of  background  from  another  frame 
of  the  spot  before  sending  the  data  to  the  HP-1000.  The  processing  here  will 
use  the  subroutines  of  the  prediction  program  PDT  to  calculate  the  Cn> 

Most  of  the  software  has  been  written  for  other  computers  and  needs  to  be 
translated  to  the  FORTRAN  used  by  the  HP-1000.  This  will  be  completed  before 
the  test,  but  may  not  be  debugged  on  the  HP-1000  by  then. 


DEPARTMENT  OF  THE  ARMY 
HEADOUARTERS.  U.  S.  ARMY  TEST  AND  EVALUATION  COMMAND 
ABERDEEN  PROVING  GROUND.  MARYLAND  21005 
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SUBJECT:  RDTE  Methodology  Improvement  Program  Directive,  Field 
Validation  of  MTF  Measurement,  TECOM  Project 
No.  7-C0-RD4-EP 0-006 


Commander 

US  Army  Electronic  Proving  Ground 

ATTN:  STEEP-MT-T 

Ft.  Huachuca,  AZ  85613 


1.  Reference  TECOM  Regulation  70-12,  dated  1  June  1973. 

2.  This  letter  and  attached  STE  Forms  1188  and  1189  (End  1)  constitute  a 
directive  for  the  subject  investigation  under  the  TECOM  Methodology 
Improvement  Program  1W665702D625. 

3.  The  MIP  at  Enclosure  1  is  the  basis  for  headquarters  approval  of  the 
subject  investigation. 

4.  Special  Instructions: 

a.  All  reporting  will  be  in  consonance  with  paragraph  9  of  the  reference. 
The  final  report,  when  applicable,  will  be  submitted  to  this  headquarters, 
ATTN:  DRSTE-AD-M,  in  consonance  with  Test  Event  570/580.  Each  project  shall 
have  a  phase/project  completion  in  FY84  which  is  reflected  in  the  scheduling 
in  attached  STE  Forms  1188  and  1189. 

b.  Recommendations  of  new  TOPs  or  revisions  to  existing  TOPs  will  be 
Included  as  part  of  the  recommendation  section  of  the  final  report.  Final 
decision  on  the  scope  of  the  TOP  effort  will  be  made  by  this  headquarters  as 
part  of  the  report  approval  process. 

c.  The  addreasee  will  determine  whether  any  classified  information  is 
involved,  and  will  assure  that  proper  security  measures  are  .taken  when 
appropriate. 

d.  Upon  receipt  of  this  directive,  test  milestone  schedules  as 
established  in  TRMS  II  data  base  will  be  reviewed  in  light  of  known  other 
workload  and  projected  available  resources.  If  rescheduling  is  necessary,  a 
letter  citing  particulars,  together  with  recommendations,  will  be  forwarded  to 
Commander,  US  Army  Test  and  Evaluation  Command,  ATTN:  DRSTE-AD-M,  with  an 
information  copy  to  DRSTE-TO-O,  within  15  calendar  days. 
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DRSTE-AD-M 

SUBJECT:  RDTE  Methodology  Improvement  Program  Directive,  Field 
Validation  of  MTF  Measurement,  TECOM  Project 
No.  7-C0-RD4-EP0-006 

e.  The  HQ,  TECOM  point  of  contact  is  Mr.  Norman  E.  Pentz,  ATTN: 
DRSTE-AD-M,  AUTOVON  283-2170/2375. 

f.  FY84  RDTE  funds  in  the  amount  of  $50,000  have  been  authorized  for  this 
investigation,  DA&COM  Form  1005  will  be  forwarded  by  the  Comptroller. 

FOR  THE  COMMANDER: 


1  Enel 
as 


GROVER  H.  SHELTON 
C,  Math  Imprv  Div 
Analysis  Directorate 


June  1983 


METHODOLOGY  INVESTIGATION  PROPOSAL 


1.  TITLE.  Field  Validation  of  TECOM' s  Path-Weighted  KTF  Determination  for 

EO  and  DEW  Systems 

2.  CATEGORY.  Modern  Weapons*  Smart  Munitions 

3.  TNSTAr.r.ATTOM .  us  Army  Electronic  Proving  Ground 
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5.  STATEMENT  OF  THE  PRQm.FM.  The  distribution  of  the  refractive  index 
structure  coefficient,  c”s along  an  atmospheric  optical  path  significantly 
affects  laser  beam  quality*  infrared  signatures  and  the  effectiveness  of 
laser  directed  energy  weapons  (LDEW) .  Although  the  feasibility  of  a  TECOM 
conceived  procedure  for  measuring  the  properly  path  weighted  MTF  from  which 
the  path  weighted  C^2  integral  can  be  calculated*  has  been  demonstrated  in 
the  laboratory,  its  validity  under  field  test  conditions  has  not. 


6.  BACKGROUND.  The  degrading  effects  of  atmospheric  turbulence  on  the  per¬ 
formance  of  beam  forming  and  imaging  optical  systems  have  been  well  known  for 
many  years.  Measured  variances  in  laser  beam  intensity*  position*  and  spot 
energy  distribution  have  been  quantitatively  correlated  to  the  degree  of  tur¬ 
bulence  along  the  relevant  optical  path.  During  the  1960's*  work  by  several 
investigators*  notably  Fried*  established  the  theoretical  relationship  between 
the  atmospheric  modulation  tr.ansfer  function  (MTF)  and  the  square  of  the  re¬ 
fractive  index  structure  coefficient  C^-CC^2).  Subsequently*  other  workers 
developed  the  theoretical  relationshipNbetween  C^2  and  laser  beam  degradation/ 
distortion  by  the  atmosphere.  The  MTF  generally  provides  a  measure  of  optical 
image  integrity  remaining  after  passage  over  an  optical  path  while  C^2  is  a 
measure  of  the  atmospheric  turbulence  along  the  path.  Thus*  from  a  measurement 
of  either  C^2  or  the  MTF*  both  image  quality  and  laser  beam  quality  can  be 
derived.  Because  of  the  critical  time  and  specific  path  sensitivity  of  the 
turbulence  characteristics*  existing  methods  for  measurement  of  C2  such  as 
scintillometry  were  deficient  in  establishing  the  true  path  turbulence  structure 
affecting  a  test  item's/system's  performance.  To  resolve  this  deficiency*  the 
Naval  Postgraduate  School  (NFS)  developed  a  mechanical  slit  scanning  technique 
which  directly  provides  the  MTF  from  which  the  properly  path  weighted  C^2  can 
be  derived.  Also*  whereas  the  scintillometer  precision  exceeds  an  order  of 
magnitude  of  the  measurement  itself*  the  NFS  system  is  precise  within  approxi¬ 
mately  10%.  During  the  review  of  the  NFS  work  by  the  TECOM  Laser  Technical 
Committee*  it  was  determined  that  a  video  camera  could  replace  the  slit  scanning 
arrangement  and  the  MTF  could  be  determined  by  the  video  image  analysis  tech¬ 
nique  already  established  at  WSMR  and  YPG.  This  would  provide  a  real  time 


MTF  measurement.  As  a  result  a  methodology  investigation  was  initiated  at 
USAEPG  to  develop  this  video  adaptation.  .The  NPS  was  employed  as  contractor 
to  do  the  actual  work  which  was  completed  as  far  as  demonstrating  that  in  the 
laboratory#  the  video  technique  provides  essentially  the  same  results  as  the 
NPS  mechanical  slit  scanning  equipment.  To  complete  this  work#  a  field  vali¬ 
dation  test  is  required. 

7.  GOAL.  To  validate  a  TECOM  conceived  method  for  measuring  the  MTF  in  real 
time  along  the  optical  path  involved  in  the  field  test  of  any  electro-optical 
(BO)  or  directed  energy  weapon  (DEW)  system  and  thereby  determine  the  effects 
of  the  atmosphere  on  system  performance. 


a.  The  validity  of  TECOM  'a  real  time#  video  based  measurement  of  the  MTF 
along  an  atmospheric  optical  path  will  be  investigated  by  a  field  test/com¬ 
parison  of  results  obtained  with  those  obtained  using  the  NPS  mechanical  slit 
scanning  equipment. 

b.  USAEPG  personnel  will: 

(1)  Develop  a  field  test  configuration  and  identify  all  instrumenta¬ 
tion  not  currently  available  but  required  for  the  field  validation  of  the  TECOM 
real  time  MTF  measurement  procedure.  The  configuration  will  be  that  appropriate 
for  the  testing  of  a  laser  designation  system.  The  standard  will  be  the  static 
MTF  measurement  system  developed  by  the  NPS. 

(2)  Procure  or  obtain  on  loan  all  equipment  required  but  not  available 
at  USAEPG. 

(3)  Measure  the  path  weighted  MTF  a  statistically  significant  number 
of  times  using  concurrently  both  the  TECOM  procedure  and  the  NPS  procedure 
along  as  nearly  identical  optical  paths  as  feasible.  Contractor  assistance 
may  be  required. 

(4)  Adapt  and  utilize  the  image  analysis  software  developed  at  wsmr# 
YPG#  NPS#  or  equivalent  to  reduce  the  data  for  determination  of  the  properly 
path-weighted  C^  integral. 


9.  JUSTIFICATION. 

a.  Association  with  Mission.  TECOM  has  the  responsibility  to  measure  the 
performance  of  varied  types  of  EO  systems  including  laser  DEW.  Optical  tur¬ 
bulence  is  one  of  the  major  causes  of  EO  system  degradation  and  should  be  pro¬ 
perly  measured  during  tests  to  enable  a  valid  performance  assessment  of  the 
EO  test  item. 

b.  Present  Capability.  Limitations.  Improvement,  and  Impact  on  Testing  if 
not  approved.  TECOM  can  currently  use  scintillometry  or  task  another  organi¬ 
zation  to  use  the  NPS  system  for  optical  turbulence  assessment.  Both  of  these 
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methods  are  undesirably  limited.  This  video  capability  will  enable  a  precise# 
convenient#  and  sensitive  determination  of  the  effects  of  the  atmosphere  on 
£0  system  performance.  Zf  the  investigation  is  not  completed#  the  effects  of 
optical  turbulence  during  many  laser  or  other  EO  system  performance  tests  will 
not  be  separable  from  the  data. 

c.  Dollar  Savings.  Dollar  savings  for  testing  per  se  are  not  quantifiable 
since  MTF  data  using  existing  methods  have  largely  been  unusable. 

d.  Workload. 

Test  Program  FY  89  88  87  88  89 

Remotely  Piloted  Vehicle  X  X 

STINGRAY  - -Classified - -  - 

CCLAW  -  -  -  -  -Classified - 

ELOCARS  X 

10.  RESOURCES. 

a.  Financial. 


Dollars  (Thousands) 


In-House 

FY84 

Out-of-House 

Personnel  Compensation 

Travel 

Consultants  and  Other  Services 
Materials  and  Supplies 

Equipment 

9.0 

3.0 

3.0 

5.0 

30.0 

Subtotals 

20.0 

30.0 

Total 

50.0 

b.  Explanation  of  Cost  Categories 

• 

(1)  Personnel  Compensation.  Covers  in-house  labor  costs  for 
principal  investigator  and  other  project  support  personnel. 

(2)  Travel.  Travel  is  required  to  coordinate  the  task  with  other 
TECGM  installations.  Travel  in  FY84  may  be  revised  if  the  field  testing  is 
scheduled  for  W8MR  or  YPG. 

(3)  Consultants  and  Other  Services.  Consultant  support  will  be 
utilized  to  assist  in  performing  the  detailed  technical  feasibility  analysis 
and  to  evaluate  the  field  test  data. 

(4)  Materials  and  Supplies.  Recording  media  and  incidental  supplies 
will  be  required  to  support  the  field  testing. 
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(S)  Equipment.  The  empirical  verification  field  tests  will  require 
a  designator  source.  Current  shortage  of  fielded  items  or  spare  test  items 
indicates  that  a  simulator  unit  will  have  to  be  procured.  FY84  funding  covers 
equipment  contingencies  including  servicing  and  repair. 

c.  Anticipated  Delava.  A  delay  in  EPG's  acquisition  of  a  laser  may 
extend  the  initiation  of  the  field  verification  testing.  However,  an  alternate 
support  capability  may  be  made  available  at  WSMR  or  YPG  to  alleviate  the 
problem. 

d.  Obligation  Plan  (FY84) . 

FO  1 _ 2  3  4  TOTAL 

9.0  50 


Available 


250 

250 


500 

(2)  Resolution  of  Non-avallable  Personnel.  Not  applicable. 
11.  INVESTIGATION  SCHEDOLE. 


0  N  0  J 

In-house  - 

e 

Contract  A  -  - 

Symbols:  -  -  -  Active  investigation  work  (all  categories) 

...  Contract  monitoring  (in-house  only) 

A  Award  of  Contract 

R  Final  report  due  at  HQ.  TECOM. 

12.  ASSOCIATION  WITH  TOP  program,  it  is  expected  that  the  turbulence 
measurement  methodology  initially  developed  for  specific  commodity  items 
such  as  laser  designators  and  rangefinders  will  be  incorporated  into  the 
existing  TOP'S  covering  those  items.  As  more  generalized  procedures  are 
developed  that  apply  to  a  broad  class  of  E-0  test  items,  then  a  new  TOP 
will  be  generated  that  can  be  referenced  in  the  E-0  commodity  item  TOP'S 


Obligation  Rate 
(Thousands) 

In-House  Personnel. 

(1)  Requirements. 


Physicist.  GS-1310 
Elec  Engr.  GS-0855 


3.0 


28.0  .10.0 


Number 

1 

1 


FY84 

Manhours 

Required 

250 

_ 


Total 


500 
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